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AN  IMPROVED  HELIUM-3  NEUTRON  SPECTROMETER 


William  A.    Bair 

Lawrence  Radiation  Laboratory- 
University  of  California 
Berkeley,    California 


ABSTRACT 

3 

This  report  is  a  preliminary  analysis  of  a     Ke  neutron  spectrometer 

being  developed  at  the  Lawrence  Radiation  Laboratory,    Berkeley.     The 

spectrometer  will  be  capable  of  presenting  neutron  spectra  to  an  upper 

3 

limit  of  14.0  MeV.     The  spectrometer  employs  a     He-Kr  gas  mixture  in 

3  3 

a  proportional  counter,    and  the  reactions     He(n,  p)T  and     He(n,  d)D  are 

3 

utilized  to  obtain  the  incident  neutron  energies.     Recoils  from     He,    Kr, 

and  the  stabilizing  gas  (CO2)  are  electronically  removed  from  the  spectra 
by  taking  advantage  of  the  shorter  rise  times  of  the  recoil  pulses.     An 
analytic  analysis  of  the  rise  time  discrimination  pulse  rejection,    the  wall 
effects,    and  the  spectrometer  efficiency  is  presented. 

The  computer  program  for  the  removal  of  the  above  corrections  is 
lengthy,  and  although  an  outline  of  the  methods  of  machine  computation 
is  presented,    this  concept  is  not  carried  to  fruition.     Instead  a  computer 
program  and  an  "inversion  matrix"  are  developed  (based  on  the  analytic 
analysis  presented)  by  means  of  a  "random  probability  analysis."     The 
"inversion  matrix"  should  permit  the  correction  of  a  raw  neutron  spectrum 
output  from  a  multichannel  analyzer  to  a  final  neutron  spectrum  by 
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correcting  in  a  single  operation  for:     (a)  rise  time  discrimination, 
(b)  wall  effects,    and  (c)  spectrometer  efficiency. 

Appendices  include  a  table  of  comparative  physical  and  operational 

3 
characteristics  of  all     He  proportional  counters  and  ionization  chambers 

constructed.,  and    reported  in  the  technical  literature  to  date  (Fall  1965). 

The  complete  electronic  circuitry  necessary  to  construct  the  rise  time 

discriminator  units,    and  the  computer  program  "CATMAN"  developed 

for  the  random  probability  analysis  are  also  presented  as   separate 

appendices. 
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I.     INTRODUCTION 
A.      Background  and  Purpose 

This  report  concerning  an  improved  3  He  neutron  spectrometer  is 
a  preliminary  report  concerning  the  theory,    design,    and  construction 
of  a     He  neutron  spectrometer  by  the  Health  Physics  Department  of  the 
Lawrence  Radiation  Laboratory.     The  actual  operation  and  performance 
of  the  spectrometer  will  be  the  subject  of  a  later  report  by  the  Health 
Physics  Department  and  will  be  authored  by  Mr.    Wai-Kit  Quon. 

The  detection  of  neutrons  and  especially  the  measurement  of  their 
energy  have  from  the  first  attempts  presented  a  more  difficult  problem 
than  the  other  common  nuclear  particles.     This  difficulty  is  because  the 
usual  methods  depend  on  effects  resulting  from  charges  on  the  particles. 
Most  methods  of  neutron  detection  depend  upon  the  scattering  of  particles 
after  a  collision  with  the  neutron,    and  the  measurement  of  the  energy 
deposited  by  the  scattered  charged  particles.     The  drawback  to  any 
method  that  is  designed  to  measure  neutron  spectra  which  is  based  on 
scattering  lies  in  the  conversion  of  the  definite  energy  of  the  incident 
neutron  to  a  continuum  of  energies  of  the  recoiling  particle.     To  avoid 
this  continuum  of  energies,    a  neutron  induced  nuclear  transformation 
was  sought  in  which  the  energy  release  would  correspond  unambiguously 
to  the  neutron  energy.     Investigations  of  several  potential  reactions, 
and  the  elimination  of  those  that  had  sharp  resonances,    low  lying  energy 
states  in  the  residual  nucleus,    and  large  positive- -or  any  magnitude 
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negative  Q  values,    led  to  the  selection  of  the  following  reaction: 
1nQ  +  3He2— ^3HX  +  1H][  +  0.  764  MeV. 

In  addition,   the  %e(n,  p)T  cross  section  is  fairly  large  and  varies 
smoothly  with  energy,    and  the  ability  to  discriminate  against  y  rays 
enhances  this  selection  for  the  task  of  neutron  spectroscopy.     Now,    if 
the  charged  particles  produced  by  the  neutron  interaction  were  to  be 
completely  stopped  in  the  active  volume  of  a  proportional  counter,    the 
ionization  caused  by  the  loss  of  particle  energy  would  be  equal  to  the 
energies  of  the  charged  particles.     This  energy  would  then  appear  as 
free  electrons  and  positive  ions  within  the  active  volume.     If  an  elec- 
trical field  gradient  is  present,   the  electrons  would  migrate  to  the 
anode,    and  the  charge  deposited  by  the  electrons  would  be  proportional 
to  the  energy,    E     of  the  incoming  neutron. 

Several  experimenters  have  examined  the  neutron  interactions  with 

the     He  nucleus,    and  have  constructed  neutron  detectors,    and  to  a 

1-13 
limited  extent,    neutron  spectrometers  utilizing  proportional  counters 

and  ionization  chambers.     The  work  of  the  majority  of  the  previous 

experimenters  has,    in  general,    concerned  itself  with  the  neutron  energy 

spectra  of  only  a  few  MeV,    although  Sayres  reports  success  with  neutron 

energies  up  to  8.  1  MeV   "  and  the  measurement  of  a  cross  section  at 

5,  14 
17.  5  MeV  using  monoenergetic  neutrons. 

In  the  majority  of  the  previous  work  accomplished,    two  problems 
seem  to  inhibit  the  usefulness  of  a  -'He  neutron  spectrometer. 
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First,    an  ambiguity  arises  from  the  fact  that     He  has  a  relatively 

3  3  - 

large  cross   section  for  elastic  scattering.     The  elastic     He(n,  n)    He 

cross  section  is  equal  to  the  nonelastic  cross  section  at  .  072  MeV,    but 
increases  rapidly  with  increasing  energy.     At  4.  2  MeV  the  elastic  cross 
section  is  4.  6  times  as  large  as  the  inelastic  cross   section.     The  ratio 
decreases  slightly  with  further  increasing  energy,    but  the  elastic  cross 
section  is  at  least  double  the  inelastic  cross  section  from  0.  25  MeV  to 
14  MeV.     The  probability  of  anundesired  recoil  pulse  is  thus  at  least 
twice  the  probability  of  a  desired  inelastic  collision.     Unfortunately,    the 
recoil  events,    in  addition  to  being  more  probable,    also  produce  a 
continuum  of  energy  pulses  from  3/4  E     to  zero  depending  upon  the  angle 
that  the  incoming  neutrons  strike  a  -'He  nucleus.     This  prevents  utilizing 
any  simple  subtraction  process  in  the  determination  of  the  initial  neutron 
spectra. 

Second,     the     possibility  of  a  nuclear  event  occurring  in  an  orienta- 
tion such  that  the  charged  particles  produced  by  the  event  intersect  the 
walls  of  the  containing  vessel  before  expending  their  full  energy  was 
recognized.     This  phenomena  gives  rise  to  a  count  that  is  of  less  energy 
than  a  count  that  goes  full  track  in  the  gas  medium.        He  nonelastic 
collisions  with  high  energy  neutrons  produce  tracks  that  are  of  consider- 
able length  in  relation  to  the  dimensions  of  the  proportional  counter,    and 
the  wall  correction  effects  increase   significantly  with  increasing  neutron 
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Solutions  to  the  above  two  problems  were  attempted  in  various 

8,9,  15 
manners.      Brown  constructed  a  "wall-less"  counter  composed 

of  a  peripheral  ring  of  1  6  small  counters  surrounding  a  central  count- 
ing wire  and  an  anticoincidence  circuit  to  remove  the  wall  effects.     He 

then  developed  a  computer  program  to  unfold  the     He  recoil  events. 

16 
Brown's  work  was  extended  by  Wang       who  developed  a  computer  program 

to  unfold  the  "wall  effects"  in  addition  to  the  "recoil  events.  "     Sayres 

10 
and  Coppola        developed  a  "risetime  discriminator"  circuitry  that 

removes  the  recoil  events  (and  a  fraction  of  the  desired  events)  by  taking 

3 
advantage  of  the  shorter  track  length  of  a     He  recoil  in  comparison  with 

the  larger  track  lengths  of  the  inelastic  particles  produced  by  neutron 

3  3,4 

interaction  with  a     He  nucleus.      Batchelor  uses  an  approximate  method 

of  computation  to  remove  the  wall  effects  which  appears  to  be  useful  at 

low  energies,    and  short  range  of  track  in  relationship  to  the  counter 

dimensions. 

The  purpose  of  this  paper  is  to  attempt  to  combine  the  removal  of 

3  3  10 

the     He(n,  n)    He  recoil  effect  developed  by  Sayres       with  the  unfolding 

of  the  wall  effect  developed  by  Brown     and  Wang.  It  is  further  the 

purpose  of  this  paper  to  report  attempts  to  extend  the  energy  range  of 

the     He  proportional  counter  to  14.  0  MeV  (which  appears  to  be  about  a 

practical  upper  limit  for  these  devices)  and  to  develop  a  means  whereby 

the  spectrometer  can  be  utilized  with  some  degree  of  4  jt  efficiency.     The 

spectrometer  will  be  used  by  the  Health  Physics  Department,    Lawrence 
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Radiation  Laboratory,    Berkeley  for  radiation  surveys—hence  the 

requirement  for  4rt   sensitivity. 

B.      Description  of  the  Spectrometer 

3 
The  experimental  apparatus  for  this     He   spectrometer  consists 

of  a  proportional  counter,    electronic  components,    and  a  computer 

program. 

1 .       Proportional  counter 

The  proportional  counter  is  a  cylindrical  tube  with  an  active  length 

of  15  inches,    a  4-inch  outside  diameter,    and  a  wall  thickness  of  0.  035 

inch.     The  diameter  of  the  anode  wire  is  0.003  inch.     The  proportional 

counter  is  designed  with  independently  adjustable  field  tubes.     Gas 

filling  is  10  atmospheres  of  krypton,    2  atmospheres  of  high  purity     He 

(less  than  10         %  tritium),         and  1/2%  CO?  as  a  stabilizing  gas.     The 

tube  was  specially  constructed  by  the  Texas  Nuclear  Corporation  for 

this  spectrometer.     A  line  diagram  showing  the  tube  and  pertinent 

dimensions  are  shown  in  Fig.    1. 

0.035  in.    stainless  steel  tubin 
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Fig.    1..        He  proportional  counter  tube 
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A  discussion  of  the  design  parameters  and  the  design. variables 
is  included  in  Appendix  A. 
2.       Electronic  components 

The  electronic  components  are  designed  to  be  able  to  distinguish 

3  *  3 

between  a  pulse  generated  from  a     He  recoil     and  a     He  inelastic  event. 

To  understand  the  operation  of  the  electronics  (a  detailed  description  | 

of  which  appears  in  Appendix  A),    it  is  first  necessary  to  understand 
what  the  electronics  are  to  accomplish.      Basically,    the  length  of  track 
of  a  JHe  recoil  (or  any  heavy  particle  recoil),    is  dependent  upon  the 
mass  of  the  incident  particle,    the  mass  and  charge  of  the  recoil  particle, 
the  energy  of  the  incident  particle,    the   scattering  angle,    and  the  composi- 
tion of  the  absorbing  medium.      (An  analysis  of  the  collision  dynamics 
involved  are  included  in  Appendix  B.  )    In  general  it  can  be   stated  that 
for  equal  recoil  energies,    a  heavier    particle,    or  a  particle  of  the  same 
mass  but  a  higher  atomic  number,   will  traverse  a  shorter  distance  in 
the  same   stopping  medium  than  a  particle  of  less  mass.     In  a  propor- 
tional counter,    a   track's    radial  component  from  the  center  wire  determine 
the  length  of  time  that  it  takes  to  collect  the  charge,    i.  e.    the  pulse  rise 
time.      Since  the  track  length  of  the  recoils  is  always  less  than  the  track 
lengths  of  the  disintegration  particles  (which  are  the  desired  events),    it 
is  possible  to  select  a  rise  time  that  will  effectively  exclude  all  recoil 


''"The  arrangement  also  serves  to  eliminate  recoils  from  the  krypton  gas 
and  the  CO2. 
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events,    yet  retain  a  certain  portion  of  the  desired  events  which  will 
have  a  longer  rise  time  because  of  their  longer  track  lengths. 
Unfortunately,    the  rise  time  criterion  also  eliminates  many  "good 
pulses"  whose  tracks  are  orientca  with  respect  to  the  center  wire  in 
such  a  way  that  their  rise  times  are  within  the  discrimination  criterion. 
The  electronic  circuitry  developed  in  Appendix  A  is  designed  to 
accomplish  this  pulse  rise  time  discrimination. 
3.      Computer  program 

The  computer  program  is  developed  to  accomplish  the  "unfolding" 
of  the  wall  effects  and  to  restore  those  desired  pulses  which  were 

electronically  eliminated  by  rise  time  discrimination.     The  wall  effects 

1  ft 
computer  program  follows  the  method  outlined  by  Wang,         but  is  not 

identical.     Differences  are  primarily  in  the  computation  of  the  probability 

function  P(E— >  E^)   which  is  the  probability  per  unit  energy  that  the 

reaction  products  of  energy  E  will  deposit  energy'E^  inside  the  sensitive 

volume  due  to  the  wall  effect.     In  addition,    a  new  probability  function  is 

analytically  developed  to  replace  mathematically  the  valid  pulses  which 

were  eliminated    due  to  rise  time  discrimination.     An  inversion  matrix 

is  proposed  which  will  combine  the  two  probability  functions,    and  enable 

a  complete  reconstruction  of  the  neutron  spectra,    above  a  minimum 

energy  that  is  determined  by  our  rise  time  discrimator  setting.     The 

neutron  spectrum  below  this  minimum  energy  cutoff  level  is  unfortunately 

lost  in  the  process  of  recoil  elimination. 


Unfortunately,    the  computer  program  analytically  developed  in 
this  paper  requires  so  much  computer  time  that  it  was  not  considered 
warranted  to  carry  it  to  completion.     Instead,    a  random  probability 
analysis  was  carried  out,    and  the  inversion  matrix  was  developed  on 
the  basis  of  these  random  probability  calculations.     An  outline  of  the 
random  probability  method  is  included  in  Section  II-E  and  Appendix  E. 
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II.     DEVELOPMENT  OF  THE  3He  SPECTROMETER 
A.      Helium- 3  and  Neutron  Interactions 
A  neutron  can  interact  with  a  -*He  nucleus  in  five  possible  ways. 

These  reactions  which  are  energy  dependent  are  as  follows: 

3 

1.  He(n,  p)T  for  neutrons  of  all  energies 

2.  He(n,n)    He  for  neutrons  of  all  energies 

3 

3.  He(n,  d)D  for  neutrons  of  energies  >  4.  36  MeV 

3 

4.  He(n,  pn)D     for  neutrons  of  energies  >  7.  32  MeV 

5.  He(n,  p2n)    H     for  neutrons  of  energies  >  10.  3  MeV 

The  cross   sections  for  these  reactions  are  reasonably  well  known 

18 
in  the  energy  range  0.  001  eV  to  14.  0  MeV       with  the  exception  of  the 

■2 

He(n,  pn)D  for  which  no  data  are  available.     It  is  assumed  that  this 
exceptional  case  has  a  small  cross  section,    even  though  it  is  energetic 

cally  possible.     A  comparative  chart  of  the  other  reaction  cross  sections 

3 

is  included  as  Fig.    2.     In  addition  to  the     He(n,  pn)D     reaction,    the 

3  1 

He(n,  p2n)    H     reaction  will  not  be  considered  in  the  determination  of 

either  wall  effects  or  rise-time  discrimination.     The  reasons  for  this 

omission  are: 

a.  At  the  maximum  energy  expected  to  be  encountered  in  this  proportional 
counter  (i.e.   14.0  MeV),    the  (n,  2n)  cross  section  is  only  15%  of  the  total 

non-elastic  cross  section  . 

3 

b.  The  neutron  induced  disintegration  of  the     He  nucleus  involves  four 

separate  product  particles --only  two  of  which  are  charged.     Four  particle 
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dynamics  are  complicated,    and  the  energy  imparted  to  the  various 
particles  may  span  a  wide  range  of  variables.     Note,    however,    that  the 
threshold  for  this  reaction  is  10.  3  MeV  and  the  kinetic  energy  of  the 
charged  particles  produced  by  even  a  14.  0  MeV  incoming  neutron  cannot 
exceed  3.  7  MeV.      Since   some  of  the  energy  of  the  incoming  neutron  will 
undoubtedly  be  carried  away  as  kinetic  energy  of  the  two  product  neutrons, 
it  is  highly  probable  that  the  sum  of  the  charged-particle  kinetic  energies 
will  be  less  than   3.  7  MeV.     Pulses  from  a  14.  0  MeV  neutron  (n,  2n) 
reaction  would  thus  appear  on  a  multichannel  analyzer  at  less  than  3.  7 
MeV.     Referring  to  Fig.    2,    it  is  to  be  noted  that  the  cross  section  of  the 
-*He(n,  p)T  reaction  is  fairly  constant  between  about  0.  4  and  2.  5  MeV,    and 
the  cross  section  at  these  energies  is  at  least  14  times  as  great  as  the 
(n,  2n)  reaction  at  14.  0  MeV.     Distortions  of  the  neutron  spectra  would 
thus  appear  above  10.  3  MeV  at  essentially  0%  at  threshold,    up  to  a 
maximum  of  15%  at  14.  0  MeV.     This  distortion  would  also  manifest 
itself  as  a  lower  energy  pulse,    (*^3.  7  MeV)  but  the  allocation  of  reaction 
particle  kinetic  energies  would  result  in  an  energy  distribution  that  would 
not  be  unique,    but  rather  span  from  0  to  3.  7  MeV.     For  a  continuous 
neutron  spectra  this  "lower  energy  distortion"  is  not  considered  serious 
due  to  the  more  probable  (n,  p)  reaction  at  these  lower  energies  which 
would  essentially  mask  these  degenerate  pulses  from  higher  energies. 

If  we  were  to  consider  only  single  interactions  within  the  active 
volume  of  our  container  (i.  e.    no  multiple .  scattering  events),    and  to 
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Fig.    2.      Energy  Dependent     He  Neutron  Cross  Sections 
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likewise  consider  no  neutron  disintegration  interactions  with  the  material 
in  the  container  walls,   we  would  find  that  (after  recoils  are  eliminated 
by  rise-time  discrimination  and  wall  effects  mathematically  resolved) 
three  reactions  which  would  still  tend  to  disturb  our  spectrum  (irrespecti\ 
of  all  corrections).     These  reactions  therefore  limit  the  maximum 
efficiency  of  our   spectrometer.      The  three  interactions  are: 

1.  The  disturbance  of  the  energy  spectrum  above  10.  3  MeV  and  below 

3  1 

3.  7  MeV  from  the     He(n,  p2n)    H     reaction. 

2.  The  disturbances  of  the  energy  spectrum  above  7.  32  MeV  and  below 

6.  68  MeV  from  the     He(n,pn)D     reaction  (if  this  reaction  does  in  fact  occu: 

3.  The  disturbance  of  the  energy  spectrum  at  all  energies  due  to  the 
neutron  induced  krypton  disintegrations. 

The  cumulative  error  introduced  by  neglecting  these  reactions  is 
unknown. 

The  remainder  of  this  paper  is  to  be  divided  into  four  parts.     The 
plan  of  "unfolding"  the  raw  spectra  to  develop  a  complete  neutron  spectrun 
is  as  follows : 

Development  of  a  probability  function  to  replace  the  pulses  eliminated 
by  "rise-time  discrimination.  "  (II-B) 

Development  of  a  probability  function  to  remove  the  wall  and  end 


sjs 

This  problem  is  discussed  in  detail  in  Appendix  A,    although  no  attempt 
is  made  to  correct  the  spectrum. 
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effects.     (II-C) 

Combination  of  the  two  correction  factors  to  reconstruct  the  neutron 
spectra.      (II-D) 

Reasons  for  and  development  of  a  random  probability  inversion 
matrix  in  lieu  of  the  above  II-D.      (II-E) 

Finally,    it  should  be  mentioned  here  that  the  work  leading  to  the 
mathematical  expressions  for  wall  and  end  effects  is  primarily  a  restate- 
ment  of  the  work  of  Dr.   W.    R.    Brown     and  Mr.    H.    T.    Wang      ,    although 
some  modifications  have  been  made  to  the  analysis  of  the  later  author. 
A  number  of  their  drawings,    graphs,    and  formulas  are  repeated  in  this 
text  for  continuity  and  to  facilitate  explanations.     They  will  be  designated 
by  a  "dagger"  (f)  if.  attributed  to  Dr.    Brown,    and  by  a  "double  dagger" 
(7* -p) if  developed  by  Mr.    Wang. 

B.     Development  of  a  Probability  Function  to  Replace  Pulses 

Eliminated  by  "Rise-time  Discrimination" 

As  a  result  of  the  arbitrary  orientations  of  the  disintegration  particles 

3  3 

of  the     He(n,  p)T  and  the     He(n,  d)D  reactions,    some  tracks  will  be  posi- 
tioned such  that  the  time  difference  between  the  arrival  of  the  first 
electron  and  the  last  electron  to  the  central  counting  wire  will  be  quite 
short.     This  situation  will  occur  when  the  radial  component  of  the  track 
lengths  is  small,    or  when  the  event  takes  place  near  the  central  portion 
of  the  counter  where  the  electrical  field  gradient  is  high.     In  both  of  these 
orientations  the  count  produced  by  a  "valid  event"  will  be  rejected  by  the 
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electronic  .circuitry  as  though  it  were  a  recoil  pulse.     It  was  considered 
that  if  an  analytical  expression  could  be  developed  for  determining  the 
time  difference  between  the  arrival  of  the  first  and  the  last  electron  to 
the  central  counting  wire,    the  "rise  time"  of  the  pulses  could  be  deter- 
mined mathematically.     A  probability  function  could  then  be  developed 
based  on  the  registered  counts  that  would  be  representative  of  those  good 
events  that  were  rejected  with  the  recoils.     To  develop  this  probability 
function,    it  is  necessary  however,    to  make  four  assumptions: 

(1.)   The  time  required  for  the  charged  particles  to  deposit  their 
energy  in  the  gas  medium  is  instantaneous. 

(2.)   A  nuclear  event  can  occur  with  equal  probability  within  any  unit 
volume  in  the  sensitive  region  of  the  proportional  counter. 

(3.)    The  mobility  of  the  electrons  in  the  gas  .mixture  is  directly' propor- 
tional to  the  field  gradient  over  the  limited  range  of  the  electrical  field 
gradients  in  which  we  are  interested. 

(4.)    The  electron  mobility  in  the  gas  mixture  (2  atmospheres     He, 
10  atmospheres  Kr,    and  0.  5%  CO2)  can  be  computed  by  assuming  that  the 


"As  an  approximation,    it  can  be   shown  that  by  relating  a  specific  particle 
energy  ^i   =    2  rrw^     ;   solving  for  v^;  and  then  relating  an  average  velocity 
to  track  length  (^(E))  and  solving  for  time  ti  =  /  v    _  ()) 

2 

that  the  time  required  for  a  particle  to  deposit  its  energy  is  about  3  nano- 
seconds for  a  proton  (longest  track  length,    highest  velocity)  of  10.  0  MeV 
energy,    and  0.  6  nanoseconds  for  a  10.  0  MeV  %e  recoil  particle  (shortest 
track  length,    smallest  velocity).     Rise  times  of  interest  to  our  counter 
circuitry  are  in  excess  of  3  microseconds. 


-15- 

electrons  behave  as  ions,    and  that  the  mixture  mobility  can  be  computed 

as  can  the  "mixture  mobilities"  of  ions. 

Assumptions  (3)  and  (4)  will  be  di        >ssed  in  greater  detail  in  the 

following  sections. 

1 .      Electron  migration  time  from  any  point   r    in  coaxial  cylindrical 

geometry.       In  a  cylindrical  proportional  counter,    the  electrical  field 

gradient  can  be  determined  at  any  point   r   from  the  center  of  the  counter 

by  the  relationship:  v 

E(r)  =  2 


r&i  —  (1) 


a 


where  E(r)  is  the  field  gradient  at  any  radial  distance    r   from  the 

anode. 

V       =    applied  voltage 

r         =    radius  of  the  cathode 
c 

r         =    radius  of  the  anode 
a 

(This  formula  is  developed  in  Appendix  D.  ) 
For  our  counter  dimensions  we  can  plot  the  field  gradient  at  any^ 
radius  from  the  center  wire  by  using  the  appropriate  applied  voltage  and 
counter  dimensions.     Instead  of  directly  plotting  the  field  gradient  as 
a  function  of  radius  however,    let  us  divide  both  sides  of  Eq.    1  by  the 
total  pressure  of  the  counter  filling  gas  in  cm  of  mercury  and  plot  this 
value  instead.     The  reasons  for  this  will  become  evident  shortly.     This 
plot  can  be  seen'in  Fig.    3  for  representative  counter  voltages. 

An  electron  produced  by  ionization  at  a  radial  position   r   will  drift 
slowly  towards  the  anode  under  the  influence  of  the  electrical  field  gradient. 
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Fig.    3.     Variation  of  Field  Gradient/Filling  Pressure  (E/P) 
for  Representative  Applied  "Voltage 
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As  can  be  seen  in  Fig.    3,    the  electrical  field  gradient  at  regions 
beyond  about  1   cm  from  the  cathode  varies  very  slowly.     After  the 
electron  passes  this  point  in  its  inward  radial  drift,    the  field  gradient 
sharply  increases.     If  we  assume  a  direct  proportionality  between  the 
drift  speed  and  the  E/P  ratio,    we  can  obtain  an  expression  in  cylindrical 
geometry  for  the  time  it  requires  for  an  electron  to  reach  the  cathode,  i.  e. 


t=         ^—  (2) 


19 
where  m,    the  electron  mooility  is  the  slope  of  the  speed  vs  E/P  curve. 

Measurements      '         of  drift  speed  for  ^He  and  Kr  with  0.  5%  CO     indicate 
that  within  the  limited  range  of  E/P  which  is  the  most  important  in  deter- 
mining the  drift  time  of  electrons  for  cylindrical  proportional  counters, 
that  this  assumption  of  linearity  does  not  appear  to  be  unreasonable. 
Substituting  the  expression  for  field  strength  in  the  case  of  cylindrical 
coaxial  geometry,    we  obtain  the  electron  drift  time  from  any  radial 


position   r.  ,r 

I 

t  = 


_dr =   1_     _P_    ln    lc_ 

V0        \  m    V0  r 

r"'        rn/- 


a  [r*Pln  r 


a 


1       P    ft     -rc        r2 

t  -     —    tt-    -<-n     —        T 

m    ^o  ra     [2> 


o; 


a 


t^I     -L     _P    £n    l£   .    r2 
2  2      m      Vo  ra 

Since    _£l    is  very  small,    the  time  necessary  to  travel  from  the  surface  of 

the  anode  to  the  geometrical  center  is  negligible. 
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The  development  of  Eq.    2  and  3  follows  the  logic  developed  by 
Harling.  ^  9 

We  still  must  determine  however,  some  means  of  combining  the 
mobilities  of  the  electrons  in  a  "gas  mixture"  since  data  is  available 
only  for  mixtures  of  individual  noble  gases  and  moderating  gases.     A 

theoretical  development  of  this  "gas  mixture"  mobility  problem  was 

21 
presented  by    Staub,        however  the  necessary  gas  constants  are  not 

available  to  compute  the  final  result.     We  must  therefore  make  use  of 

assumption  (d)  and  assume  that  as  a  first  approximation,    the  electron 

mobilities  behave  as  do  ionic  mobilities. 

While  it  is  realized  that  no   simple  theory  can  account  for  the  electron 

diffusion,    attachment,    recombination,    and  agitation  velocity,    as  the 

electron  moves  through  a  changing  field  gradient,    it  is  considered  that 

this  approximation  will  cause  at  most  an  error  of  a  constant  factor.      The 

electron  mobilities  are  thus  assumed  to  follow  the   same  physical  laws 

as  do  ions --although  it  is  realized  that  the   situation  is  considerably  more 

complex  for  the  fast  moving  electrons. 

21 

Staub1  s  method  of  combining  ion  mobilities        is  as  follows: 

>^12  P         Pl  P         P2 

where  ^2     =    the  ion  mobility  in  the  gas  mixture 

P-l    2  =    ^e  *on  mo1°ility  in  the  component  gases 
P^    2  =     partial  pressures  of  the  component  gases 
P         =     P     +  P2  =  total  gas  pressure. 
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Substituting    m    (the  electron  mobility)  for   u   we  have: 


PHe  1  PKr  1 


+  ^—      *      —  (4) 


mmix         ptot  mHe    ptot  mKr 

Using  the   sources  explained  in  the  footnote,  ''"  we  find  that  for  our  tube 

V  _l 

mHe        =     0.3cm//j-sec-(^-cm        Hg) 


m 


Kr        =1.0  cm/^  sec-(  -^cm"1  Hg) 


Ppje         =     2.  01   atmosphere 
Pr,  =     10.05  atmospheres 

Substitution  of  the  above  values  into  Eq.    4  yields  m     .      =  0.  72 

n  7  mix 

,.-y    -  -i 

cm/u-sec-(cm  •  cm    Hg).  Utilizing  this  value  in  Eq.    3  and  the  appropriate 
counter  dimensions  we  can  plot  migration  time  for  an  electron  from  any- 
radial  point  in  the  counter.      This  has  been  accomplished  in  Fig.    4. for 
representative  counter  operating  voltages. 

2.     Rise  time.      Rise  time  of  the  pulse  has  been  defined  as  the  time 
difference  between  the  arrival  of  the  first  electron  at  the  central  counting 
wire  and  the  arrival  of  the  last  electron.       By  use  of  assumptions  (c)  and  (d) 
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The  mobility  factor  for  Kr  was  taken  from  the  data  of  English  and  Hanna. 
The  mobility  factor  for  the  helium  is  taken  from  the  data  of  Friedes  and 
Chrien,  ^  1  which  is  plotted  in  reference  19  (Harling).     The  reason  for  this 
use  of  Friedes  and  Chrien'  s  data  rather  than  Harling' s  is  that  Harling, 

(1)  made  his  measurements  using  4He  with  only  small  amounts  of  ^He,    and 

(2)  he  used  substantial  percentages  of  CO2.     Although  drift  velocity  of 
electrons  in  ^He  appears  to  be  rather  insensitive  to  moderating  gases,    the 
CO2  content  of  Friedes  and  Chrien  was  considerably  less  than  that  utilized 
by  Harling. 
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Fig.    4.     Variation  of  Electron  Migration  Time  from  Various 
Radial  Positions  for  Representative  Voltages 
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we  have  been,  able  to  determine  migration  time,    and  have  formulated 
an  expression  for  the  time  for  an  electron  to  transit  any  radial  distance. 
Thus  the  rise  time  (t    )  can  be  expressed  as 

t,  =  t2  -  h    -    I   I     f-  in  ^       frf2  -    A  (5, 

r         L        l  Z   m     V  ra       \i  n) 

Note  that  this  is  identical  to  the  time  that  it  takes  an  electron  at 
position   rr       (furthest  distance)  to  migrate  to  rn  (nearest  distance)  if  we 

charge  the  limits  of  Eq.    3  from    r    to    r„    and    r       to    r_.     We  have 

i  f  a  n 

previously  noted  that  the  field  gradient  charges  very  slowly  beyond  about 
1  cm  from  the  anode  wire  (see  Fig.    3).     Therefore,    if  we  select  our  rise 
time  discriminator  in  such  a  manner  that  pulses  which  spend  their  full 
energy  within  the  active  volume  enclosed  by  a    1   cm  radius. are  rejected. 
We  can  see  that  the  assumption  regarding  the  linear  dependence  of 
mobility  and  E/P  is  not  too  serious  a  concept  since  the  pulses  which 
depart  from  this  assumed  linear  dependence  (at  high  E/P  values)  are 
automatically  eliminated.     As  it  turns  out,    the  "dead  volume"  of  our" 
counter  will  extend  to  about  2.  75  cm,    as  will  be  shown  later.     Pulses 
which  originate  in  the  "dead  volume"  and  extend  outwards  will  be  treated 
in  detail  later,    as  will  pulses  which  do  the  reverse.      (Section  II-B.  7) 
3..     Electron  production.     From  the  above,    we  now  have  a  means  of 
calculating  rise  times  for  electrons  that  are  located  at  various  radial 
distances  in  the  counter  active  volume.     We  do  not,    however,    have  any 
indication  of  where  these  electrons  will  be  located—in  fact  there  will  be 
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millions  of  electrons  generated  along  the  path  length  of  a  charged 
particle  that  is   stopped  in  the  gas  medium.      Let  us  now  turn  our  attention 

to  the  formation  of  these  electrons  in  a  gas  medium.     We  have  previously 

3 
discussed  the  various  reactions  that  can  occur  with  a     He  neutron  inter- 
action,   and  have  determined  that  only  three  of  the  five  potential  reactions 
are  of  interest  in  our  counter.     Neglecting  the  energy  distribution  within 
the  three  events  for  the  moment,        let  us  now  determine  the  ranges  of 

the  various  particles  in  the  gas  medium. 

22 

Friedlander  and  Kennedy        present  an  empirical  relationship  for  the 

determination  of  the  ranges  of  charged  particles  in  a  gas  mixture.     This 

relationship,    adapted  for  our  specific  usage,    can  be  used  to  determine 

3 

the  ranges  of  protons,    tritons,    deuterons,    and     He  recoils  as  follows: 


Range  formulas  for  the  individual  gases  (mg/cm    ) 
Protons 


Rp  in  Kr  (E)      =    Rp  in  air  (E)    [l .  89  -  0.  25  log10  (E)]   -  0.  36' 
Rp  in  He  (E)      =    Rp  in  air  (E)    [o.  82  +  0.  043  log10  (E)] 
-'He  Recoils 


(6 


{E]  [1.89  -0.25  1og10(!]]-0.36J 
R3He  in  He  (E)=    |    Rp  in  air    (^j[o.82  +  0.043  log10  (•£)] 


Ro         in  Kr  (E)  =    A    /Rtd  in  air 
°He  4     I    c 

'( 


'The  energy  distribution  is  discussed  in  Appendix  C  and  in  Section  II-'B.  6, 
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Tritons 


RT  in  Kr  (E)      =     3  }Rp  in  air  (£)   [l .  89  -  0.  25  log10  (  j]  ]  -  0.  3 6> 

J(6c) 


Rrp  in  He  (E)      =     3  R-^  in  air 

Deutorons 


(f)    [0.82  +  0.  043  logoff)] 


Rd  in  Kr  (E)       =     2  |rp  in  air  ( ^  j   [l .  89  -  0.  25  log10  [  —  )   ]-  0.  36> 

/E\    r  (Ml  '(6d) 

Rd  in  Ke  (E)       =     2  Rp  in  air    ( ^-     [0.  82  +  0.  043  log10  [4f-    J 

The  particle  ranges  can  be  converted  into  centimeters  by  use  of  the 

relationship: 

„  Rtot    (rag/cm2) 

I(cm)=-^- — —  (7) 

fHe    PHe+    PKr    PKr 

\14.7/  \14.  7/ 

where  2.  is  the  particle  range  in  centimeters  in  the  gas  mixture, 
P    =  partial  pressure  of  the  individual  gases  (psi) 
P    =  gas  density  at  STP 
R.    .    =  total  particle  range  in  the  gas  mixture     (see  below) 

1  _  WHe  WKr 


*WE)  RHe(E)  RKr<E) 

W    =  weight  fraction,    W^je  +  W^;r  =  1 

PHe     PHe 


(8) 


W 


He"       PHe     PHe    +    PKr     PXr 


wKr=      Kr    pKr —      =    1  -  WH( 

PHe    PHe    +    PKr    ?Kr 
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where 

P    =    partial  pressure  of  gas  (psi) 

P     =     gas  density  at  STP 

•°He  =    0.  16629  g/liter  (f) 

PKr  =    3.455  g/liter  (7) 

The  above  relationships  enable  us  to  conveniently  compute  the 

I 
ranges  of  the  charged  particle  in  the  gas  mixtures,    and  are  amenable 

to  machine  computation.     The  ranges  computed  as  above  are  plotted  in 

Fig.    5  for  comparison.     Note  the   substantial  difference  in  the  length  of 

the  track  for  various  particles  of  the  same  energy--particularly  above 

about  1.  0  MeV.     It  is  this  difference  in  track  length  that  will  enable 

us  to  remove  the  recoil  spectrum  from  our  spectrometer. 

4,     Rise  time  discriminator  values.      The  rise  time  discriminator   setting 
for  the  pulses  can  now  be  determined  since  we  now  have  a  means  of 
determining  the  pulse  rise  times  and  the  length  of  track  of  various 
particles  in  the  gas  medium.     If  we  recognize  that  it  is  only  the  radial 
component  of  the  track  length  that  is  important  in  our  determination  of 
rise  time  descrimination,    we  can  select  the  most  unfavorable  recoil 
situation  that  may  be  encountered.     Setting  the  discriminator  at  the  rise 
time  appropriate  to  this  most  unfavorable  case,    we  will  remove  from 
consideration  all  pulses  that  can  be  attributed  to  recoils  (He,    C,    O,    and 
Kr).     This  worst  case  can  be  considered  to  be  composed  of  two  factors; 

(1)  The  maximum  energy  that  can  be  imparted  to  a  recoil  nucleus,    and 

(2)  the  longest  rise  time  that  this  maximum  energy  recoil  can  cause,    due 
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Fig.    5.     Ranges  of  Various  Charged  Particles  in  Gas 
Mixture  of  10  Atmos.    Kr,    2  Atmos.    3He,    and 
0.  5%  CO. 
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to  its  geometrical  position  in  the  counter. 

Appendix  B  (Case  1)   summarizes  the  maximum  energy  that  can  be 
imparted  to  the  various  recoil  nuclei.     It  can  be  readily  seen  that  the 
3 He  nucleus  can  absorb  up  to  a  maximum  of  3/4  the  energy  of  the  incidei- 
neutron.      The  other  possibilities  for  energy  absorption  are  of  consider- 
ably smaller  magnitude.     The  ~>He  recoil  is  thus  selected  as  the  recoil 

of  concern.      The  ranges  of  the  other  recoil  particles  in  the  gas  mixture 

3 
is  likewise  less  for  equal  energies  when  compared  with     He.     The  worst 

geometrical  case  that  can  occur  is  if  the  track  length  of  the  recoil  ^He 

nucleus  is  oriented  in  a  directly  radial  position  with  one  end  of  the 

track  at  an  infinitesimal  distance  from  the  wall  of  the  cathode.     This 

positioning  is  illustrated  in  Fig.    6. 

Rise  time  discriminator  setting  can  then  be  determined  by: 

(1)  calculating  the  maximum  recoil  energy  that  can  be  imparted 

Eo  (max)  =  0.75  E.,  (max) 

JHe  in 

(2)  determining  the  path  length  in  the  gas  medium, XJ^->         (max))from 

JHe 

the  equations  presented  in  Part  II-B-j.2. 

(3)  determine  maximum  rise  time  by  substituting  into  Eq.    5  the  value 
b,    the  radius  of  the  cathode,   for  r^,    and      [b  -i(Eo         (max))]    for  r    . 

£"16 

This  rise  time  discriminator  setting  (t„  )  will  then  exclude  all  recoils 

°  v  r  ' 

from  the  resultant  spectra. 
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Fig.    6.        He  recoil  track  position  to  give  maximum  pulse 
rise  time. 


5,       Determination  of  radial  positions  of  the  nearest  (rn)  and  furthest  (rf) 

3  -  3 

electrons  from  the  anode  for      Ke(n,  p)T  and     He(n,  d)D  reactions. 

Consider  the  fact  that  to  define  the  positioning  of  any  event  in  space,    it  is 

necessary  to  define  six  degrees  of  freedom.  -     Three  of  these  degrees  of 

freedom  are  necessary  to  define  the  positioning  of  the  center  of  gravity, 

and  three  are, required  to  explain  the  orientation  of  the  event  with  respect 

to  the  center  of  gravity.     For  the  purposes  of  discussion,    let  us  assume 

that  we  can  define  the  "center  of  gravity"  as  the  position  where  the  neutron 

strikes  a     He  atom.     Since  the  reaction  occurs  within  the  boundaries  of 

the  counting  tube,    we  can  define  the  positioning  of  the  event  (center  of 

gravity)  in  cylindrical  coordinates  of  radius,    length,    and  angular  dependence. 

The  operation  of  the  counter,    however,    is  independent  of  the  positioning 

of  the  event  in  terms  of  length  (except  for  end  effects,   which  will  be 

discussed  later)  and  completely  independent  of  any  angular  relationship, 
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with  respect  to  the  anode  wire.      This  means  that  the  positioning  of  the 
center  of  gravity  can  be  defined  in  terms  of  a  single  parameter    r. 

To  define  the  orientation  of  the  resultant  ionization  tracks  produced 
by     He  disintegration  let  us  first  turn  our  attention  to  the     He(n,  p)T 
reaction.     For  a  neutron  of  given  energy  En,    and  given  resultant  proton 
energy  E0,    we  can  determine  the  energy  K<-n  of  the  triton  particle  and  the 
laboratory  angle  aT    between  the  proton  vector  and  the  triton  vector  by 
use  of  the  relationships  developed  in  Appendix  B. 


ET        "    EN  +  Qn,  p  "  EP  ( 1 


and 

(2ET  +  Qn    p) 
cosaT   = —  (1C 

Z*1'  3  Em  Ep 

Note  thato;  is  always  obtuse  in  the  laboratory  system  of  coordinates. 

Let  us  now  define  a  coordinate  system.     The  coordinate  system  for 
our  calculations  is  selected  as  follows: 

(1)  The  event  occurs  at  a  perpendicular  radial  distance  r~  from  the 
center  line  of  the  anode.  The  end  point  of  Tq  defines  the  zero  point 
(0,    0,    0)  of  the  coordinate  axis. 

(2)  The    z    axis  is   selected  to 'include  the  vector    r    and  extensions  thereof 

(3)  The   x  axis  extends  parallel  to  the  anode  wire  and  passes  through 
the  point  (0,    0,    0). 
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(4)     The    y   axis  passes  through  the  point  (0,  0,  0)  defined  as  above, 
and  is  mutually  perpendicular  to  the    x   and    z  axis. 


(0,0,0) 


■  x 


anode  wire 


Fig.    7.    .Selection  of  coordinate   system. 


With  this  selection  of  coordinate  axis,'  the  event  is  angularly  independent 
of  the  center  anode  wire.     Calculations  of  the  path  orientations  of  the 
ionizing  tracks  will  be  made  in  terms  of  polar  coordinates  using  these 
axes  as  reference  lines.     We  must  now  determine  the  distances  of  the 
nearest  and  farthest  electrons  from  the  anode  wire.     For  the  purpose 
of  illustration  let  us  visualize  the  Ep  and  Erp  track  lengths  as  solid  lines. 
The  orientation  of  these  two  vectors  can  have  three  degrees  of  freedom. 
Let  us   select  the  three  rotational  axes  as  shown  in  Fig.    8. 

Neglecting  for  the  moment  the  ET  vector,    let  us  turn  our  attention 
to  the  Ep  vector. 

We  are  able  to  define  a  position   r^  which  we  call  the  tip  of  the  Ep 
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Rotation  about 
the    x   axi  s 


Asmuthal  rotation 
around  the    z    axis 


-x 


wire 


Rotation  about  Ei 


vector 


Fig.    8.     Rotational  axes. 

vector.  .    (We  will  refer  to  the  Ep  vector,    although  in  reality  we  are 
speaking  of  a  definite  length  of  path  of  a  proton  of  set  energy  in  the  . 
counter  gas  medium.  )    We  can  completely  describe  this  point   r_    with 
respect  to  our  coordinate  axis  by  the  polar  angle  92,  the  asmuthal  angle 
0~  and  the  length  of  the.  vector  Ep.     We  can  thus  determine  r ^  (Ep,  62,  $?) 


92     /      Ep    r2(Ep'92'02) 


-0. 


x 


Fig.    9.     r2(Ep,e2,02)  . 

We  need  an  expression  of  r-?  in  terms  of  Tq   since  in  determining  migra- 
tion times  we  are  concerned  with  the  radial  distance  from  the  central 
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Ep[l-sin2  52  cos2-^] 
Ep  sin92sin(-^>2) 
Ep  sin  8Z  cos  (-#_) 


1/2 


Plone   conjoining  x-oxis  and  vector   Ep 

Vectors  lie  ina  plone  parallel  to  yz  plane1 


mu. 37200 


Fig.    10.      Determination  of  ^(r      E    ,  0_,  Q^) 


\    - 
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counting  wire  to  the  point   r2    (Ep,  0,,  02),    i.  e.   we  need  an  expression 
of   r?(rn,  Ep,  92,  02).     To  determine  this  expression  let  us  pass  a. plane 
which  contains  the   x    axis  and  the  vector  Ep  (a  plane  can  be  defined  by- 
two  intersecting  lines).     This  plane  will  intersect  the    z    axis  at  (0,  0,  0) 
and  will  make  an  angle  {3  with  the    z    axis  in  the  minus    y   direction. 
Referring  to  Fig.    10  we  can  derive  an  expression  for     r2(rQ,  Ep,  Q~,  02). 
3y  the  law  of  cosines: 

r     =  r^  +  E^|l  -  sin2  02  cos2(- 02)]  -  2rQEp[l  -  sin2  02  cos2(- 02)]         .cos(18C 

but,    cos  (180  -    )  =  -  cos   |3 

Ep  cos  92 

and  from  Fig.    10,    -cos    £= -r#2    , 

Ep  [l   -   sin2  92  cos^  (-  02)J 

substituting  this  value  for  cos  £  into  the  above  equation,   we  find  that 

2         2  2  r  7  2  -, 

r     =  r     +  E     i  1  -  sm     92  cos     (-  02)J  +  2  rQ  Ep  cos  Q^ 

or  (I! 

r2(r0,Ep,  02,  02)  =  /r2  +  Ep[l  -  sin2  02  co  s2  (- 02)j  +  2  rQ  Ep  cos  92 


7ro+E^[1- 


Let  us  now  consider  the  possibility  that  since  we  have  not  placed  a 
restriction  on  where  our  event  can  occur,    that  a  skew  line  may  develop 
and  that  the  closest  point  to  the  anode  wire  may  in  fact  lie   somewhere 
along  the  length  of  vector  Ep,    rather  than  at  the  end  point.     Since  Ep 
is  a  definite  length  in  the  gas  medium  for  a  given  energy  of  the  proton, 
let  us  call  this  length ii(Ep).     An  axiom  of  solid  geometry  is  that  the 
closest  point  of  approach  for  skew  lines  is  at  a  point  where   a  line 
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connecting  the   skew  lines  is  perpendicular  to  each.      This  means  that 

a  minimum  distance  can  be  determined.     We  have  an  expression  for  the 

length    x 2^     Therefore,    if  we  take  a  partial  derivative  of   r2   with  respect 

to  Ep   (note  that  this  is  actually  J^(E-o))and  set  it  equal  to  zero,    we  can 

determine  a  distance  A_p  along  the  Jc(Ep)  vector  whereby  the   skew  lines 

are  closest.     We  will  call  the  closest  distance  between  the  vector  Ep 

and  the  anode  wire  r  .. 

4 

rl/2 
„  r,^ — •    =    —    :  r     +,v(Ep)      [l  -  sin     0?  cos     (-  97)J  +  2£(E-q)  rn  cos  9?  I 

\    ')  I  7  •>  \ 

•  j2^(Ep)  J  1  -  sin^  02  cos^  (-  02)  j    +2rQ  cos  9- 

Set  the  above  expression  equal  to  zero  and  solve  for  IL-p\   note  that  the 
first  bracketed  term  can  never  equal  zero.   For  real  values  of  SL-p  we 


r0  cos  y2 
XP  =     "    [1  -  sinz  02  cosz  (-02)j  (12) 

where    0  <  i p  <i(Ep). 


The  negative   sign  indicates  that  i2p  cannot  be  positive  unless  92  is 
greater  than  90°.      This  is  in  accordance  with  the  physical  situation 
envisioned  and  the   selection  of  the  coordinate  axes.     At  92  =     re,  JCp  must 
equal  Tq  and  the  Ep  vector  is  perpendicular  to  the  anode  wire.    Xp  cannot 
exceed  the  physical  length  of  Ep  which  has  been  assigned. 
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Substitution  of  this  value  of  JL-p  into  Eq.    11  for  Ep  determines 


that: 


r2  cos2  02  •     [l  -  sin2  02  cos2  (-  02)]    -  2  rg  cos2  92 

'r     +. 


0     t-  sin202  cos2 


i    v2 
'2'J 


[l  -  sin     02  cos2  (-02)] 


or 


r4=        r0' 


cos      02 


[I  -  sLi2  02  cos2  (-  02)j 


where    —    S.Q-    <     ^ 
2  2 


and         0      <ip    <    &(EP) 

3 

Let  us  now  turn  our  attention  to  the  other  portion  of  the     He  (n,  p)T 

event.      Let  us  consider  the  Erp  vector. 

In  a  manner  similar  to  the  determination  of  the  end  of  the  proton 
vector  Ep,we  can  position  the  end  of  the  triton  vector  E-p.     We   shall 
call  this  position  r    (Erp,  0-,  ,  0,  ) 


r. 


(180  -  0i)i, 


£j-n. 


r^ExOi^x)^      / 


.x 


—  wire 


Fig.    11.     r1(ET,01,01) 


In  the  case  of  the  triton  however,   we  cannot  arbitrarily  describe  the 
position  of  the  end  of  the  E-p  vector  with  respect  to  the  anode  wire,    since 
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there  exists  a  unique  relationship  for  the  positioning  of  the  triton 
vector  for  a  given  neutron  and  proton  energy.    (Ej\f  and  Ep  respec- 
tively. )     This  relationship  is  characterized  by  the  angle  OJ  t  ,    which  is 
the  obtuse  angle  between  the  proton  and  triton.     We  must  therefore 
consider  this  constraint  upon  the  direction  that  the  triton  can  travel. 
We  can  visualize  this  restriction  by  the  following: 


ANODE 


Fig.    12.     Rotation  of  Erp  'vector  about  Ep  axis. 


If  we  were  to  rotate  the  Ep  vector  in  Fig.    12  about  its  own  axis 
while  imagining  that  the  E      vector  is  rigidly  attached  to  the  Ep  vector 
at  an  angle  aT  .  from  the  axis  of  the  Ep  vector,    we  can  see  that  the  tip 
of  the  ET  vector  would  trace  a  circle  in  a  plane  that  is  perpendicular 
to  the  Ep  vector.     The  determination  of  the  radial  distance  to  the  center 
wire     r-i  (r«,  Erp,  9-.,  0  •.)  must  take  this  degree  of  freedom  into  account. 

Let  us  now  examine  Fig.    13.     The  interior  angle  of  the  cone  described 
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by  the  rotation  of  ET  is  (180  -  O'.-r  ).      Let  us  now  pass  a  plane  through  the 
z  axis  and  containing  the  Ep  vector.     Since  this  plane  contains  the    z    axi 
it  will  by  definition  be  perpendicular  to  the  plane  containing  the   x   and   y 
axis.     We  will  use  this  plane  to  relate  9.    to  9?,  and  0^  to  $?'      ^et  us 
consider  the  rotation  of  the  E^  vector  from  a  "zero  position"  where  the 
Erp  vector  is  at  a  maximum  position  away  from  the  center  wire  (i.  e. 
top  position)  and  is  contained  in  the  (z,  Ep)  plane.      The  Ep  vector  is 
then  rotated  in  a  counter  clockwise  manner  (viewed  from  the  tip  of  the 
Ep  vector).     The  angle  of  rotation  we  call  7S '.     By  projecting  the 
Erp  vector  back  to  the  plane  containing  Ep  and  .z   we  can  determine  the 
incremental    angle    between  0  •,    and  (180  -9?)  this  incremental  angle  is 


sin 


I     |  E      sin   Q^  cos  y 


[_ET  (1 


•in2  «L  sinz  7)1/2i 


or 


.    -1 
sin 


sin^L    cos    y 
_(1  -  sinZ  aL  sin2  7)1/z 


therefore  Q-,    is  related  to  Q?  by 

r 


_^  sin  O^  cos  7 


h    =  (180  -  97)  J_  sin"        g * TT7 

1  2  !(1  -  sin2  aL  sin2/)1/2 


or 


x  =     180  -  02  ±  sin*1 


sine;  l  cos  7 


,(1  -  sin^a 


L 


•   2    a/2 

sm   7  )        . 


(3 


where   a.    as  determined  by  E-^T  and  Ep  is  a  constant,    7  varies  from 


0to2jr.      The   sign  of  the  incremental  angle' is  positive  for   -  J£  <     7   <^L    and 

2  2 


negative  for   2L  <    7    <".£.   # .      92  is  arbitrarily  selected. 
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Troce  of  circle  formed  by 
tip  of  Ej  vector  after 
roration   about  the  Ep  vector 


Ej  (l-sin2  aL 


Anode 


Note:  Vector  ET  shown   exaggerated  in  size  with 

respect  to  Ep  for   clarity  of  angular   relationships. 


mu.  37201 


Fig.    13.     Determining  Relationship  Between  0£  and  9,,    and 
02  and  0, 
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If  we  likewise  were  to  consider  the  perpendicular  projection  of  the 
tip  of  the  E^,  vector  to  the  (Ep,  z)  plane  we  can  readily  determine  that 
since  this  projection  is  perpendicular  to  the  (Ep,  z)  plane,    it  must  of 
necessity  be  parallel  to  x  -  y  plane.     We  may  therefore  compute  the 
incremental  angle  between  0i    and  0_,    and  determine  a  relationship 
between  them.     The  incremental  angle  is 

1            ET  (1  -sm^c;L  sin      y) 
cos  — 

Therefore,    the  relationship  between  0-t    and  0?  can  be  expressed  as 

(180  -0X)  =  {180  -02)  -  cos"1  jjl  -  sin2aL  sin27)j1/2 

or 

+  ™"1    f/i  _.,-«2^_    miJL„  n1/2 


0:  =  02  Icos         [(1  -  sin2aL  sin27  )1/2J 


where  the  negative   sign  of  the   square  root  is  taken  for  0\  y   <  tf    and  the 
positive   square  root  is  taken  for      rt     <y  <2^ 

CC-r    is  a  constant  determined  by  fixed  Ep  and  E-nt 
0o    is  arbitrary 
Having  now  established  a  relationship  between  Q±  and  O2,    0]_  and  02 
through  the  constant  angle  c^and  the  rotational,  angle    y,    let  us  now 
proceed  to  determine    r,  (rQ,  ET,  9i  ,  0,  ). 

Referring  to  Fig.    14,    the   same  procedure  can  be  followed  as  was 
used  for  the  computation  of  r2(r      Ep,  9^,  02)  in  Fig.    10.      We  arrive  at 
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ET  sin  fi|  sin(l80- 


Ej  sin  5 1 
—  Ej  cos  0| 
ET[l-sin2a,  cos2  (180- 


S  r|(ET,©|,^|) 


Anode 


mu. 37202 


Fig.    14.     Determination  of  r    (r      E       9^  0    ) 


•  40. 


the  relationship 
2_„2   ,  ^2  f, 


.2  n 


r     =r„  +E       1  -  sin^  9,    cos     vlSO  -  0,).' 
i        0  T  L  i  J-.J 

Er  cos 


2rQE     [l  -  sin2 


91   cos2  (180  -  0J 


but  from  Fig.    14,    cos  (o  = 


)„  [1  -  sin2  0      CosZ  (180  -  0i)j 
r  [  1  \J 


7TT1/2 


substituting  this  value  into  the  above  relationship  we  find  that 


r.    =     /r0  +  E2    fl  -  sin2  0,    cos2  (180  -  0,  )1    -  2  rn  ET  cos  0, 
1       J     0  t  L  1  v  1  .j  0      T  1 


( 


The  possibility  of  a  "skew"  condition  must  likewise  be  considered 
in  case  of  the  E^  vector.      The  length  [y-ry)  must  be  determined  and 
substituted  into  Eq..l6.     The  partial  derivative  of  r.   with  respect  to 
.X.(Erp)  must  be  taken  and  set  equal  to  zero, 
thu  s : 
■^r^  1      f    2 


6JL(ET)        2  0 


f(ET)    [l  -sin2  ©-,    cos     (180  -  0  )]  -2r0i(ET)cos 


"    j  2J:(ET)   [l-  sin     01  cos     (1 
Solving  forjtrp  we  find  that 


2  M     1 

"0  -01)J    -  2  r0  cos  0 


=     0 


JLr*     = 


rn    COS 


[l  -  sin2  0X  cos2  (ISO  -0i)j 


for        0<^T<i(ET) 

substituting  this  value  back  into  Eq.    16  and  calling  this  closest  point 

of  approach  r,,    we  find  that 
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2       i-q  cos2  01->  f  1  -  sin29-i    cos2  (180  -  0-,)]         2  rj  cos     91 


3         °    [i  -  sin2  Ql   cos2  (180  -  01/;'J  [l  -  sin2  01   cos2  (180  -  0 x )J 


or 


(li 


/-  CO  s 

r  _    =    r  _  /i 


3     -     o  n  _-^2  Q     „     2 


[l  -  sin"  9X  cos2  (180  -  01)1 


Summarizing  the  calculations  for  radial  distances  from  the  center 

wire,    given  E~T.    E„,    and  r„  we  must  calculate  four  additional  radial 
to  N         p  0 

distances  for  each  Ep,    E_  vector  pair. 


ET       =    EN  +  Q  -  Ep  (9) 

r0         =    r0 


r2  =    A/r2  +  Ep  [l  -  sin2  92  cos2  (- 02)]+2  rQEp  cos  02  (11) 


2  o 
cos     92 

r4         =    r0,/    l'-_ -f" 2  ,     r/    ,-. (13) 

i  -  sin     92  cos      (-  02) 


/  ;  i  -  sin     y?  cos     \-  w?r 

A/  L  ^  ^ 


(only  necessary  to  compute    r^_  where    —  <_92^ 
and  0<^p  <  E-o) 


rx  =    /r2  +  Ei-[l  -  sin2  01   cos2  (180  -  0  J\     -Z  tqEt  cos  Q1  (16) 


cos     9-1 
r         =    r„    /  1 


3  0  /  '  r,  2  „  2   , 


[l  -  sin2  91   cosZ  (180  -  0X)J  (18) 


anci : 

9^      -     iou- «2 


,         j  i        sin<Aj_,  cos  2T                                                (2EttQ) 
=     lSO-O^lsin1] _T __  wiierecosd=-  i (14) 

"  L  ' 


II  -  sin    <3*T    sin  <T;  2/3  x^rp  E- 
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The  sign  of  the  incremental  angle  is  as  follows:     plus  for  -  ±ir<"^<7l, 
minus  if    —(V^-LIL 


0     =     02    -    cos"1      ("l-sin2   6LL  sinVj  1/2  (15. 

The  sign  of  the  incremental  angle  is  taken  as  negative  for  0<  V  <77' 
and.  positive  for   7T   <  <?  <i  2    . 

It  has  thus  been  shown  that  all  possible  combinations  of  radial 
distances  from  the  center  wire  can  be  calculated  in  terms  of  known 
E-yT,    Ep,    and  Tq  in  terms  of  O2,    0?,  & ,    and  cU  . 

If  r    ,    r-,,    r->,    r,,    and  r      are  calculated  for  a  given  event,    and  the 
maximum  and  the  minimum  radial  distances  taken  as  the  positions  of 
the  most  distant  electron,    Xr,    and  the  nearest   electron  il   to  reach  the 
central  counting  wire,    we  now  have  the  distances  to   substitute  into  the 
equation  for  determining  rise  time. 
Deuteron  radial  distances 


For  the  reaction  involving  deuterons,    there  is  an  analogous  r'elation- 
ship  to  the  proceeding  if  we  replace  E^  with  Eqoj    Ern  with  E-.-,,    andd^ 
with  the  relationship:     (see  Appendix  B) 

(Qn.d+  ED2  +  EDl) 
4/ZDl      ED2 


coso?       =-  iif _  (19) 

JLi 
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thus,    given  E-T,    E-rs?'    anc^  ri 


EDl=EN-ED7^Qn,d  .(2°) 


r0       =    r0 


r         =      /r     +  E         [l-sin292cos     (- 92)]+ 2  rQ  ED2  cos  92      (21) 

V  °     D2 


cos2q2 

4  r0  /  r -T- ,    ,     H     , (22) 


->2  (-02)] 
(only  necessary  to  compute  r     when    —  -<  9 2 <£  7/ and 

o<^D2<i(ED2; 


1/  i  1  -  sin^  92  cos' 

ir 
2 

0<X_  <i(En?) 


rl      7r0+EDl     L1  -  si^2  Gl  cos2  (18O-01)]-2roED1cosO1       (2J) 


cosz  §i 

r^      =      rn  /  l  "  r 9 0 1      0<xDi  <^(En1)  (24) 

[l-sin2  9,    cos2  (180-0!)]  D1 

sino^      cosTT 

and         Q1       =     180  -  92  ±  sin-1     L _ I  (25) 

■|[l  -  sin2  olL  sin2£r]1/2j 

_  3^7-  _, 

The  sign  of  the  incremental  angle  is  as  follows:     plus  for     " — <  tf<!L  , 

Z  ^ 

minus  if     f^^T     and  0^   =  02    +    cos"       |"(1  -  sin2c£L  sin2^),1'2  (26) 

where  the  incremental  angle  is  positive  for  -rr     ¥  <  2<  and  negative  for 

o<  $ <ir. 
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6.       Development  of  the  probability  function  (P-n)  for  the  replacement 
of  pulses  removed  by  rise  time  discrimination. 

We  have  now  developed  a  means  whereby  the  rise  time  of  any  valid 
pulse  (i.  e.    not  a  recoil)  can  be  determined  for  any  possible  energy  distrib 
tion  and  spatial  configuration  of  the  ^He  (n,  p)T  and  the     He(n,  d)D  reaction 
If  we  were  to  vary  Ep,    r„,    cos  9,    0,    and  <f  over  all  possible  ranges, 
a  tabulation  of  rise  times    tr  can  be  made  for  each  incoming  neutron  energ- 
Before  this  tabulation  can  be  accomplished  however,    three  additional 
conditions  must  be  considered. 

(1)     The  volume  probability  of  an  event  occurring  at  a  distance    r    is  great 
at  points  that  are  distant  from  the  anode.     In  Fig.    15  below,    an  end  view 
of  the  counter  tube  is  shown.     The  volume  probability  g(r)  of  a  reaction 
occurrence  at    r    within  a  differential  element   dr    is: 


g(r)dr  =      2t1t  dr      -     2r  d] 


7Tb2  b2 


where   bis  the  radius  of  the  cathode. 


Fig.    15.'    Volume  probability  in  the  radial  direction. 


Ep  will  vary  from  Ep  max  to  Ep  min  for  a  set  neutron  energy.     Any  ene; 
within  these  limits  is  equally  possible.     Cos  9  is   selected  as  a  variable 
since  the  selection  of  9  would  result  in  undue  weighting  of  the  rise  times 
near  the  polar  angles.    A  discussion  of  these  factors  is  included  in 
Appendix  C.     It  is  assumed  that  the  applied  voltage  is  held  constant  in  the 
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Therefore,    events  occurring  at  various  distances  r   from  the  centerline 
of  the  counter  must  be  weighted  in  accordance  with  their  volume  prob- 
ability. 

(2)     Thus  far,  the  finite   size  of  our  tube  has  not  been  considered  in  the 
determination  of  rise  times.     It  becomes  obvious  that  some  of  the  energy- 
vectors  will  intersect  the  wails  of  the  cylinder,    and  that  the  ionizing 
tracks  will  be   stopped  short  of  their  full  length.     In  this  case  a  pulse  of 
lower  energy  and  rise  time  will  result.     An  additional  constraint  must 
therefore  be  placed  upon  the  calculations  of  the  various  r's  (i.e.    r.,    r?, 
r-j,    r  /)  in  the  preceding  section,    viz,    if  any  of  the  r's  exceeds  b,    then 
the  pulse  must  be  considered  as  a  "short  pulse.  "     Short  pulses  will  be 
divided  into  two  possibilities: 

(a).     The  radial  component  of  the  track  lengths  is  still  sufficient  to  cause 

t 
the  pulse  rise  time  to  be  above  the  rise  time  discriminator  level   t    ..     This 

will  cause  a  count  to  be  registered  in  a  lower  energy  channel.      The  rise  time 
is  computed  using   b   as  rr.     This  type  of  pulse  will  be  designated  as  a 
"wall  effect  pulse"  (W.  E.  P.  )  and  will  be  totaled  separately, 
(b)      The  radial  component  of  the  track  lengths  is  insufficient  to  cause  a 
pulse  to  be  registered  in  any  energy  channel.     The  rise  time  will  be 
computed  however,    using  b  as  the  maximum  of  the  r's  and  catalogued  in 
the  appropriate  rise  time  channel. 

For  "bookkeeping"  purposes  it  is  necessary  that  the  above  two  possi- 
bilities for  "short  pulses"  be  catalogued  in  the   same  energy  channel 

appropriate  to  the  (E     +  Q)  considered  for  a  pulse  that  goes  full  length 

N 
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in  the  track  medium.     They  are  not  to  be  catalogued  in  the  energy  channel 
corresponding  to  their  reduced  length.     Wall  effects  'are   separately  removl 
and  if  these  degenerate  pulses  were  removed  twice,    it  would  distort  the 
spectrum. 

Each  time  the  rise  time  discriminator   setting  is  changed  the  above 

1 
two  possibilities  ior  "snort  puxses"  must  oe  recomputed. 

(3)     To  compare  the  rise  time  spectrum  from  the  (n,  p)  and  (n,  d)  reaction: 

it  is  necessary  that  the  two  reactions  have  a  common  reference  conditions 

This  common  reference  condition  can  be  met  if  it  is   stipulated  that  the 

number  of  rise  time  situations  considered  for  each  incoming  neutron  ener  \ 

be  the  same.    The  variables  for  the  rise  time  spectrum  (for  a  set  neutron 

energy  E-^j,   fixed  rise  time  discriminator  setting,    and  fixed  applied  voltagi 

are:    three  angles  (cos  0,    0,    and  ff),    distance  Tq,    and  an  energy  distributic 

of  Ep.    The  angle  variations  are  identical  for  each  case  considered,    the 

variation  of  rQ  will  be  identical,    but  the  E  can  vary  from  an  Ep  max  to 

Ep  v-p^-f-,*    the  values  of  which  are  dependent  on  the  energy  of  the  incoming 

neutron.     We  must  therefore  consider  in  our  tabulations  that  an  identical 

number  of  energy  increments  (not  identical  energy  increments)  be  taken 

between  Ep  max  and  Ep       .      regardless  of  the  energy  of  the  incoming 

neutron.      Note  that  s^nce  only  a  ratio  is  eventually  desired,    this  condition 

merely  increases  the  numerator  and  denominator  of  the  ratio  by  a  constan 


All  three  of  the  above  conditions  are  likewise  applicable'  to  the  deuter 
reactions,    as   shown  in  the  following;  discussion. 
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A  tabulation  of  the  number  of  occurrences  of  pulses  having  the 

same  rise  time  for  each  (E-^-  +  Q)  channel  can  now  be  formed.     This  can 

be  accomplished  in  the  following  manner   (see  Table  I):     For  each 

(E._  +  Q         )  channel,    i.  e.   the  energy  of  the  incoming  neutron  plus  the 
J\  n,  p 

Q  of  the  reaction,  a  reaction  site    r      is  selected.      The  energy  distribu- 

n,  p  a 

tion  of  the  particles  produced  by  the  reaction  is  varied  in  a  number  of 

energy  increments  from  E-,  _„__-  to  E-o  w,:^,.     Note  that  the  total  energy 

available  is  constant  at  (E-^  +  Q)  since  we  have  not  varied  the  energy  of 

the  incoming;  neutron.      The  three  angular  relations  are  also  varied  over 

all  possible  configurations.     For  each  energy  increment,    and  for  each 

angular  configuration,'  the  pulse  rise  time  is  calculated  and  tabulated.     The 

number  of  pulses  that  have  rise  times  at  1  jJ.sec,    2  usee,    etc.    are 

summarized  and  multiplied  by  the  weighting  function    — r.     The  reaction 

site  is  moved  a  small  radial  distance  outward  and  the  process  is  repeated. 

This  tabulation  continues  until  the  reaction  site  is  an  infinitesimal  distance 

from  the  cathode.     The  volume  weighted  number  of  pulses  of  the  same 

rise  time  is  then  totaled; 

r=b  r=b  r=b 

)ail  counts  •    dr  = 

r=b  1      - 


)N(WEP)  •   £.    +  }N(l^isec)I.    +)N(2^sec)^    +  .  .  .    + 
r=0  r=D 


AT(30^jisec)  — 2 

>  b 

r  =  0 


dr 


j 


where  dr  can  be  any  constant  and  is  added  merely  to  remove  the  dimension: 
dependence  of  cm      .     For  finite  differences  in  r,    the  interval  between  the 
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r     and  r,    could  be  selected. 
a  b 


The  above  tabulation  is  repeated  for  all  values  of  (E^T  +  Q)  from 
E.T=  0,    to  E-^  =  14.  0  MeV.     Note  that  because  of  the  previously  imposed 
conditions,    that  the  summation  of  all  the  counts  for  each  (E^  +  Q)  tabula- 
tion are  equal. 

Tabulating  the     Ke(n,  p)T  and  the     Ke(n,  d)D  reactions   separately, 
we  can  construct  a  matrix  for  each  reaction  which  is  a  tabulation  of  the 
number  of  events  having  identical  rise  times  in  each  energy  channel. 
Such  a  matrix  is   shown  in  Table  II. 

An  identical  tabulation  to  Table  II  would  be  constructed  for  the 
(n,  d)  reaction. 

The  (n,  p)  and  (n,  d)  tabulations  could  be  directly  combined,    but  the 
reactions  are  not  equally  probable.     However,    if  we  multiply  all  the 
entries  of. the  (n,  d)  tabulation  by  a  factor  that  is  proportional  to  their 
relative  probabilities,    then  we  could  directly  add  the  corresponding 
elements  in  the  two  matrices  to  obtain  a  combined  matrix  of  rise  time 

vs.    (E     -f-Q).     Such  a  factor  is  available  to  us  in  the  form  of  the  cross 

N 

sections  of  the  two  reactions,     crn   p(Ej^)  and   o~n   ^(E.-).      (Note  that  these 
are  the  cross   sections  of    o-(E„)  and  not    o~(EN  +  Q).)     The  combined  matrix 
would  appear  as   shown  in  Table  III. 

Therefore  the  percentage  of  counts  removed  per  energy  channel  is 
equal  to  the  sum  of  the  counts  to  the  left  of  the  rise  time  cutoff  (see 
Table  III)  divided  by  the  total  number  of  counts  in  that  energy  channel. 
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If  the  percentage  of  counts  removed  is  designated  as  "Pr>  "  then  the 

percentage  of  counts  that  were  counted  is  (1  -Pr).      Thererore,    consicerix 

the  raw  spectrum  that  is  received  as  output  from  the  multichannel  analyze 

and  designating  the  counts  from  any  individual  (E-^-f  Q)  channel  as  C-^,    we 

can  say    (1  -  Pr)Crn  =  CN-  or    CRN  = — :  ,   where  CRN  is  the  valia 

**• 
counts  we  would  have  received  without  rise  time  discrimination. 

A"p  "  can  be  computed  for  each  (E-xttQ)  energy  channel  of  the  multi-  { 
channel  analyzer,  and  after  correction  we  will  obtain  the  corrected  neutro; 
spectrum  with  the  rise  time  discrimination  removed. 

7.     Discussion  of  error  introduced  by  assumption  of  linear  dependence 
of  electron  velocity  and  (E/P). 

It  has  previously  been  assumed  (Section  II- 3.-  1)  that  the  mobility  of 
the  electrons  in  the  gas  medium  is  a  constant  value  and  that  the  velocity 
of  the  electrons  in  their  inward  drift  towards  the  anode  varies  directly 
as  the  (E/P)  value,    i.e.    v  =  m  •    E/P.     This  relationship  appears  to  be  a 
reasonable  approximation  in  the  regions  where  the  field  gradient   E   varies 
slowly.     It  can  be  seen  from  Fig.    3  that  from  the  cathode^o  about  a 
distance  of  1..  0  cm  from  the  anode,    this  condition  of  a  slowly  varying 
(E/P)  holds  quite  true.     At  distances  closer  than  1.  0  cm  from  the  anode, 
the  field  gradient  increases  rapidly,   becoming  higher  the  closer  the  anode 
is  approached.     At  regions  very  close  to  the  anode,    gas  multiplication 
takes  place,    and  an  "electron  avalanche"  is  produced. 

The  region  of  interest  for  electron  migration  time  calculations, 
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however,    is  the  region  between  the  radius  where  gas  multiplication  begins 

23 

(i.  e.    a  few  mean  free  paths  of  the  anodc--or  as  described  by  Wilkinson, 

"a  distance  from  the  wire  surface  of  the  order  of  its  diameter")  and  the 
region  where  the  field  gradient  begins  its  rapid  assent  (i.  e.    about  1.  0  cm 
radius). 

First,    it  is  known  that  the  mobility  of  electrons  is  not  a  constant 
value  over  unlimited  variations  of  E/P.     At  some  value  for  each  component 
gas,    a  saturation  condition  is  reached  whereby  further  increase  in  E/P 
does  not  result  in  an  increase  in  electron  drift  velocity- -and  in  fact,    the 
velocity  may  even  decrease  with  increasing  K/P.     These  saturation 
effects  are  probably  the  result  of  electron  collisions  with  gas  molecules 
that  are  no  longer  purely  elastic,    and  a  portion  of  the  electron's  kinetic 
energy  is  imparted  to  the  gas  molecule  in  an  inelastic  collision.      The 
probability  of  an  inelastic  collision  apparently  increases  with  increasing 
electron  kinetic  energy,    but  the  increased  energy  losses  due  to  the  inelastic 
collisions  cause-the  drift  and  agitational  velocities  to  remain  at  a  constant 

(or  nearly  contant)  value.     In  the  case  of  krypton,    this  terminal  velocity  is 

1  20 

reached  at  an  E/P  value  of  1.  5  volts/ cm  •  cm      Hg        when  the  krypton  has 

1/2%  CC>2  intermixed.      English  and  Hanna   °  attribute  this  peak  velocity 

as  the  result  of  the  "Ramsauer  Effect.  "     On  the  other  hand,   no  such 

24 
terminal  velocity  peak  has  been  reported  either  in  pure  He,         or  m 

3  ll>  i9  , 

Ke-CO?  mixtures.  In  addition,    the  velocity  vs.    E/P  curves  that 

are  available  from  various  experimenters  do  not  go  beyond  an  E/P  value 
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-1  -  20 

of  more  than  6.0  volts /cm  •    cm    '    Hg  for  krypton,         and  10.  0 

volts/cm  •    cra~x  Hg      '        for  ^He,   whereas  the  E/P  value  varies 

from  O..75to  over  190  volts/cm  •    cm"     Hg  in  our  counter,    in  the  region 

under  consideration.     It  is  thus  impossible   (at  present)  to  even  estin 

the  electron  drift  velocity  in  the  inner  1..0  cm  radius  of  the  counter  tube, 

Quantitatively,    we  can  say  that  the  electron  velocity  in  this  region 

will  be  considerably  less  than  the  value  we  have  assumed.      This  phcnoml 

of  smaller  and  slowly  changing  velocity  is  due  to  the   saturation  effects  oil 

3 

the  Xr-CO_  mixture -'-and  probably  an  identical  effect  from  the     He-C09 

mixture  that  occurs  at  an  E/P  value  higher  than  those  reported. 

Let  us  now  examine  the  effects  of  this  uncertainty  of  electron  velocij 
in  this  region  in  our  counter  tube.     Equation  5   states  that  the  rise  time  all 
a  pulse  is  equal  (for  a  specified  applied  voltage)  to  a  constant  times 
differences  in  the  squares  of  two  radii,    i.e. 

t_  =  const.   •   1  r     -  r  \dj 

r  V  *        n/ 

We  have  likewise  pointed  out  that  for  tracks  that  occur  entirely  within 
the  region  where  E/P  varies   slowly,    that  we  may  consider  the  "rise 
as  the  time  that  it  takes  the  furthest  electron  to  migrate  to  a  radius  equ; 
to  that  of  the  origina'   position  of  the  nearest  electron. 

Let  us  now  go  back  to  the  determination  of  our  original  setting  of 
the  rise  time  discriminator.     The  value  of  the  rise  time  discrimination 
t     was  determined  by  solving  Eq.    5a  with  the  values    rr  =  b  and 
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u      0 


n 


3         (— E-^t).     Once  this  value  has  been  determined,   we  can 


rie 


interpret  Eq.    5a  as  the  area  between  two  concentric  circles  of  radii    r 


and  r 


n 


I  ig.    16.     Interpretation  of  rise  time  values  as  a  constant  area. 


The  total  area  between  the  circles  remains  a  constant,    i.e. 


,2         2s 
TJf(r     -  r    )  -rjrt 

f         n  r 


2  m  Jr 


V.  £n  rc 


when    r      =    0,    this  expression  reduces  to: 


rf  (min) 


2  m  P 


N 


V    in.  rc 

0        — 

a 


(5b) 


Where  rf(min)  *s  ^ne  0-uter  radius  of  the  "dead  volume,  "  and  any  tracks 
that  lie  completely,  within  the  cylinder  bounded  by  rf/min)>    they  are 
rejected  by  rise  time  discrimination.     From  the  above  discussion  we 
realize  that  ^f/-nin)  is  n°t  as  large  as  calculated  in  Eq.    5b,    due  to  the 
uncertainty  of  the  electron  migration  velocity  in  the  region    ra<r  <1.  0  cm, 
since  the  velocity  tends  to  be  considerably  less  than  the  linear  relation- 
ship originally  assumed  in  a  high  E/P  region.     For  a  14.  0  MeV  neutron, 

3 
which  produces  a  iO.  5  2v4eV  maximum  .  He  recoil,    we  find  that 
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t     =  6.1  usee  and   rf(min)    =  2.7  era  (for  Vq  =  5500  V).     Since  rf(min) 

(actual)  is  less  than  the  value  calculated  for  a  fixed  t'  ,    the  "dead  volume" 
v         .      '  r 

of  our  counter  is  less  than  calculated.      The  question  of  "how  much  less" 
is  one  that  cannot  at  present  be  analytically  solved.     It  has  been 
previously  stated  that  it  was  considered  that  the  linear  dependence  of 
velocity  with  the  (E/P)  factor  was  valid  at  regions    r>  1.  0  cm.     If  the 
rf(min)  (actuad)  is  not  less  than  1.  0  cm,    then  our  assumptions  are   still 
valid  for  those  events  whose  tracks  lie  entirely  in  the  active  volume 
beyond  I,.  0  cm  radius. 

We  have  assumed  that  the  probability  of  a  reaction  occurring  in  any 
unit  volume  of  the  gas'  medium  is  identical.      The  volume  enclosed  by  the 
1.  0  cm  radius  comprises  only  4%.  of  the  total  sensitive  volume.     The  track 
orientations  can  likewise  be  completely  at  random,    and  so  the  probability 
of  a  track  having  any  portion  of  its  length  in  the  volume  enclosed  by  a 
1.  0  cm  radius  is  dependent  upon  the  disintegration  particle  energies,    the 
reaction  site,    and  the  orientation  of  the  reaction  particles.     While  the 
probability  of  passing  into  this  "inner  wall"  region  can  be  computed  in  a 
manner   similar  to  the  probability  of  a  particle   striking  the  outer  wall 
(see  Section  II-C),    there  is  no  way  to  relate  this  probability  to  the  rise 
time  spectrum,    since  electron  velocities  are  unknown  in  the  region-- 
also  every  accepted  count  will  be  registered  as  a  "full  energy  count.  " 
In  addition,    even  if  the  tracks  of  the  charged  particles  do  enter  into  the 
region  of  the  "dead  volume,  "  their  total  radial  orientation  may  be   such 
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that  the  pulse  rise  time  is  greater  than  t   ,    and  the  counts  will  be 

r 

registered. 

It  can  thus  be  concluded  that: 

(a)  Pulses  that  are  the  result  of  tracks  lying  completely  inside  the  "dead 
volume"  will  be  excluded  due  to  rise  time  discrimination,    regardless  of 
energy  and  orientation.      (4%  of  the  active  volume.  ) 

(b)  Tracks  that  lie  completely  outside  the  "dead  volume"  may,    or  may 
not,  be  counted  due  to  rise  time  discrimination.      Those  pulses  above  a 
minimum  energy  that  were  removed  due  to  track  orientation  can  be 
mathematically  replaced. 

(c)  Those  tracks  that 'have  some  portion  of  the  length  within  the  "dead 
volume"  may,    or  may  not  be  counted.     The  actual  rise  time  of  these 
pulses  will  be  larger  than  the  calculated  rise  time,    but  the  rise  time 
discriminator   setting  will  be  the   same  as  in  (a)  and  (b)  above.      This  means 
that  more  counts  of  each  energy  wouic  be  registered  than  would  be  indicated 
by  our  theory  (i.e.    more  counts  because  more  pulses  would  have  longer 
rise  times  and  thus  more  would  escape  rejection  due  to  rise  time  discrimina- 
tion. ) 

The  magnitude  of  the  error  introduced  by  (c)  above  is  unknown,    and- its 
effects  must  await  the  comparison  of  the  final  spectrum  as  developed  from 
our  counter  with  a  spectrum  obtained  by  other  means  from  the  same  source. 
It  is  not  anticipated  that  this  error  will  be  large. 
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II-6-  DEVELOPMENT  OF  A  PROBABILITY  FUNCTION  FOR 
THE  UNFOLDING  OF  THE  WALL  EFFECTS. 

The  development  of  a  probability  function  to  remove  the  vail  effects  w? 
primarily  developed  by  Browtijand  then  extended  by  Wang.  The  overall      I 
development  is  quite  long,  being  in  general  the  substance  of  two  theses, 
and  only  certain  aspects  will  be  presented  herein  along  with  a  correction 
to  Wangs'  "Wall. Effects  Function".  For  a  complete  derivation  of  the  entire 
series  of  probability  functions,  readers  are  referred  to  Ref  v  9  aa&  l6. 
The  notation  of  Brown  and  Wang  will  be  continued  in  this  paper  to  aid 
persons  who  are  interested  in  development  of  these  functions,  but  in 
general  the' ■  j  functions  will  only  be  presented  and  identified  in  their 
final  form. 

1.  .Wall  Effects 

The  sensitive  volume  of  the  proportional  counter  tube  was  previously 
identified  in  Fig.  1.  Within  this  sensitive  volume  of  length  L  and 
radius  b,  let  us  select  a  reaction  site  at  r,  .  and  assume  that  the 
track  length  of  charged  particles  is  -<-  .  For  the  analysis,  the  whole 
gas  filled  volume  (not  just  the  sensitive  volume)  may  be  considered  to  be 
divided  into  four  regions  which  will  be  identified  in  our  subsequent 
discussion, 
a.  Region  A 

Region  A  is  the  core  of  the  sensitive  volume  defined  by  radius  of 
b  --£and  length  L  -  2-£as  shown  in  Fig.  17 


t..jJ 
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t 

L L-2A   _ 
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It  is  clear  that  only  tracks  of  full  size  can  be  found  in  this 
region.  Therefore,  the  probability  PA(P)  Per  unit  track  length  that 
a  track  -will  have  a  length  between  £  and  f  +  d(>  in  this  region  is  :: 

pa(P)#  =  S  (P  -X)dP 
For  a  single  reaction,  happening  in', region'. A '  .<•' the  probability  is 

OQ 

pa  (P)a?  =  (   -(<?  -X)dP  =  1  (+) 

b*  Region  B 

Region  B  is  the  hollov  cylindrical  volume  with  an  inner  radius  of 
b  -X,     an  outer  radius  of  b,  and  a  length  L  -  2  A--,   as  in  Fig.  18 


c-n 


t 

Figure  18   ■-  Region  B 
For  region  B,  Wang  classifies  the  track  length  into  two  categories? 
Case  1 

For  a  single  reaction  in  this  region,  those  tracks  not  intersecting 

•  '  ff 
the  wall  will  have  a  length-*4*,  the  full  sized  track.  The  probability, 

P-oi  (\  )}   Per  unit  track  length  that  a  track  will  have  a  length  between 

\°  and  (v  +  d  P)  in  this  region  is  ; 

P3L  (/°)d,.p  =  FB1  (P)ff(f  -4)  (.T) 

where  the  factor  F   (P)  represents  the  fractional  area  of  a  sphere  with 
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radius  JC   centered  at  (0,0,0)  that  falls  inside  the  sensitive  volume. 

Referring  to  Fig.  19,  it  is  seen  that  the  fractional  surface  area 
of  the  sphere  inside  the  cylinder  is 


sin  0  d  9  d  ft 


FB1(f  >- 


u«r 


1       r°2 
1  -     |       f         sinedQ 

Ql 


vhere  S  =  spherical  surface  area  outside  of  the  cylinder. 


Brora  makes  an  approximation  that  the  transition -from  9     to  9     along 
the  line  of  intersection  "between  the  sphere  and  cylinder  can  he  represented 


by  the  sin     0,  and  he  .concludes  that 


counting  vire 


Figure  19     Intersection  of  a  Sphere  and  a  Cylinder 


Bl 


<e>-^[(^>i<r2)  -ft)- 


PB1  (p  )*e   -1  +  1 


h2  -l2> 


M!  -fg-ij^-^V  (+ 


Case  2. 


Those  tracks  that  do  interest  the  detector  vail  will  lose  a  portion 
of  their  length.  The  probability  of  hitting  the  vail  is  then  the  spherical 
surface  outside  the  cylinder  i.e. 


FB2^   "   "  \ 


fba-P2 


L\    ae 


)  i  4 


-*  )-J. 


*     r 


Then 


;he  variation  over  r>    gives     3?       =   (p  ,r) 


-61- 


"B2 


-   (P,   r)=FB2   (p,.*) 


^(FB2(e  ,  r))  =  PB2  (p  ,  r)df> 


b2  -e* 


i 

8" 


2 

t-  2p    r 

i.e2 


2p 
2  p  r 

b. 
-^2  + 

5r 

2b 

e2 

3r 

}de 

dp 


c.     Region  C 


Regions  C  are  the  portions  of  the  sensitive  volume  of  distance  jl 
from  the  ends  of  the  sensitive  volume,  as  in  Fig.  20 


TT/\ 
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Figure  20  -  Region  C 
The  track  can  be  classified  into  two  different  cases. 

Case  1 

For  a  simple  reaction  in  this  region,  the  probability  of  having  a 

full  track  length  is 

pci  W  d<?   =  Fci  (e )  s  (^-f )  de 

***£   -  j\an   sin  cc  da  dp  =.  1  -  -  J  sin   a  du 


Fci  <?>  - 


■  l-l 


cos  a 


0 
-f  cos' 


1  f# 
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=  1  - 


1 
2 


Figure  21  -  Region  G 

Intersection  of  Sphere 
vith  Active  Volume 


Case  2 


The  probability  of  having  a  curtailed  track  length  is 

P_  (P)  d?  =  (1  sin  a  da)  :x  2 
°Z  2 

since     cos  a  =  z  ,   sin  a  da  =   (n2)  dP.  therefore 

PC2  k°>  dP=    ^       if  z<e 


=  0 


if   z><° 


cL-.     Region  D 

Region  D  are  the  portions  of  the  gas  filled  volume  beyond  the 
sensitive  volume,  as  inEig.,..22 '.  • 


(+) 


Figure  '£d    -   Region  D 

Nuclear  reactions  occurring  at  sites  outside  the  sensitive  volume 
can  also  cause  energy  to  be  deposited  inside  the  sensitive  volume  of  the 
counter. 

If  the  probability  of  depositing  energy  from  a  track  length  inside 
the  sensitive  volume  per  unit  track  length  of  (°  is  P  (f )  then:  ■ 
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P^  (<°)d?     =     -    sin  a  da    x  2 
D  z 


but  cos  a  = 


p^l 


therefore  sin  a  da  =   -  7-3 — n-y  dv   for   ^>  \   ,  and 


dP 


?D   ^  =  (1^2  fQr   (    l_  z)>f>?  0 


ftt) 


:ome 


Figure, -23-Region  D;   Track  Lengths  Inside  Sensitive  Volume 

e..     Summary  of  the  General  Wall  Analysis 

For  an  arbitrary  track  length  £    the  gas -filled  volume  can  he 
divided,  into     four  regions,  as  in  Fig.   'dh 


(+t) 


Figure  'dh:   Summary  of  Counter  Regions 

The  probabilities  per  unit  track  length  that  a  track  will  have  a 
length  between  p  and  (p  +  dp  )  in  each  region  are  listed  in  Table  IV. 
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\ 

f.     The  probability,    P(E  -aE-),    per  unit   energy  that  reaction  products  of 
energy  E  will  deposit  energy  E.;   inside  the   sensitive  volume  due  to  the 

wall  effect 

According  to  the  foregoing  analysis,  now,  if  we  assume  that  the  pulse 
height  is  proportional  to  the  energy  deposited  inside  the  sensitive  volume, 
then  it  is  reasonable  to  assume,   for  the  first-order  approximation,    that 

E      =    K, 
Ei    =    Kp, 
Eb   =    Kb, 
EL  =    KL, 
where   K   is  some  proportionality  constant  that  will  be  cancelled  out  shortly. 

Therefore 

P(E  -*Ei)  =    K'p(X^p), 

where   K     is  another  proportionality  constant.     For  the  whole  counter 


p(X-»p)=  pB2(p,r)gl(r)dr  +        pC2(  p,  z)g2(z)dz 

;  d  -p  Jo 


+  |  PD2^P'  z)g3(z)dz< 

0 


From  Subsec.   II-D.  lb. 


>\  2b       3rl 

"2  *  TZ}>    pC2  -  "2'     PD2 


.p  P      J  p  (-^-P) 

and  the  volume  fractions  are,    respectively, 

,    .  ,  277'rdr       2rdr 

gi(r)dr  =  — - 2~  = — T~     ' 
1  -filoc  b 


dz 
g2(z)dz  =  g3(z)dz  =  — 
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Therefore: 
P 


c*  f )  -  a  M 


b     r„      2 


+  1 


IjFe) 


2" 


CLZ 


_J 


\ 

1  2b 

;  +  r~ 


'd?  dr 

n-2 


//J     ^x  1  1  1  1  1 

+   "L  +  2L  =  2b     +  L 


20 


and  since  p  (i!— >^  )do     =  p  (E->E..  )dE 


then  P(E  ->  E.  )•  .*     ~- 

1  ^ 


E, 


(27) 


"Where  E,    is  the  energy  width  of  the  particle  in  question,  and  ET    is 

D  Li 

the  energy  length.  Note  that  this  is  the  probability  or  fraction  of 
giving  a  track  of  energy  E.  due  to  the  vail  effect.  Surprisingly,  this 
probability  turns  out  to  be  constant  independent  of  the  energy  of  the 
incoming  neutron.* 


2.  Probability  of  a  Track  Not  Kitting  the  Container  Envelope 

Brown  and  Wang  derive  expressions  for  the  average  laboratory  energy 
of  the  various  particles  in  a  plane  perpendicular  to  the  counting  wire. 
The  expressions  cannot  be  analytically  resolved,  so  numerical  integration 
is  performed  on  these  expressions  and  the  resultant  energies  are  then 
approximated  by  use  of  the  following  formula: 


Ep  =  0.505  EN+  0.450 


t+r 


*  Wangs  more  complicated  probability  function  is  due  to  the  improper 

selection  of  limits  of  integration  for  region  C  and  the  approximation 

utilized,  in  determination  of  P  (p>)  d*  . 

B2  '    V> 
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ET       =  0.^0  E^  +   0.150 


Jf 


=  0.42     EAT  -   1.27 


(+) 


Where   the  designation     E     indicates  the  average  laboratory  energy 
in  a  plane  perpendicular  to  the  anode  resulting  from  a  neutron  of  energy 
E   .       All  units  are   in  MeV. 

From  these  expressions  of  particle  energy,,   the  probability  that  a 
simultaneously  emitted  proton  and  triton,   or  a  pair  of  deuterons,   will 
not  hit  the  vails  of  the  containing  vessel  are  derived.     Using  the 
notation  of  Wang, 

f^(E)     =     probability  that  a  proton-triton  will  not  hit  the  wall 
fc(E)     =     probability  that  the  deuterons  will  not  hit  the  wall 


fl(E)      =  flp(E)      '      flT(E) 


d-*T)l/2 


(23)  (f\ 


where  f,   (E)  =.  the  probability  that  a  proton  track  will  not  hit 


the  wall* 


'IT 


-  the  probability  that  a  triton  track  will  not  hit  the  wall 
a  factor  added  to  account  for  the  forward  track 


(l-:t>1/2  = 


component  of  the  triton. 

3E 
'  K   =    ( 


JN 


Tryr —  72     where  Q         =0.764   (MeV) 
**N  +   4  Vp  ^'P 


'IP 


1 
(E) 


it  I  2b 


1   - 


f1T(E)     - 


1-  *   (M 

«      V  .2WJ 


(29)  (+) 

(30)  (t) 


*  Both  Brown  and  Wang  calculate  two  different  values  for  f-  (E)  and 
equate  them  for  energy  boundary  conditions  of  their  tubes.  *      In  our 
case,  the  energy  boundary  conditions"  (because  of  a  larger  size  tube 
and  greater  filling  pressure)  are  such  that  we  do  not  ever  reach  the 
conditions  that  Brown  and  Wang  encountered,  (i.  e.  the  E  . ,  for  our 
tube  is  at  17.O  MeV  which  is  above  the  designed  maximum  of  14.0  MeV) 
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vhere  X.     -   length  of  proton  track  in  the  gas  mixture  for  a  proton 


of  energy  E 


L- 


and  .$   =  length  of  triton  track  in  the  gas  mixture  for  a  triton 


or  energy  E^   . 


vhere 


f5(E)   - 


^ 


i  -  -  IS 

Tt     \  2d 


Tt  2^    1/2 

(1+  V 


1/2 

r      -2 


•V       J 


51)    li 


and 


/  EN 


x°=  \ 5  sn  +  *  v ; 


\l/2 


,      %A  =  5-2^   (MeV) 


J2_  =  length  of  a  deuteron  track  in  the  gas  mixture  for  a  deuteron 


■of  energy  E  . 


The  above  nomenclature  and  designation  "will  he.  used  in  subsequent 
sections. 


-69- 


D.   Combination  of  the  Two  Probability  Factors  to  Reconstruct  the 
Original  Neutron  Spectrum. 

The  original  purpose  of  our  neutron  spectrometer  was  to  accurately 
portray  a  neutron  spectrum  in  the  energy  range  from  essentially  0.1  MeV 
to  l^-.O  MeV.  The  spectrometer  operation  was  such  however,  that  in  order 

to  remove  the  more  numerous  recoil  pulses  which  were  smeared  across  the 

3  E 
energy  spectrum  from  (j-       n)   '  to  zero,  advantage  was  taken  of  the 

differences  in  pulse  rise  times  between  valid  pulses  and  recoil  events. 

The  removal  of  the  recoil  pulses,  however,  was  not  without  its  disadvantages. 

It  was  discovered  that  to  obtain  a  neutron  spectra  without  recoils,  it  is 

necessary  to  forgo  all  information  at  the  lower  energies  since  the  track 

lengths  of  lower  energy  pulses  were  shorter  than  the  maximum  energy  recoil 

pulse.   These  lower  energy  pulses  were  thus  removed  with  the  recoils.   In 

addition,  the  orientation  of  the  tracks  of  a  valid  event  may  "be  such  that 

the  radial  components  are  less  than  the  minimum  required  to  produce  a 

pulse  that  exceeds  the  rise  time  minimum  cut  off. 

The  loss  of  the  lower  energies  of  the  neutron  spectrum  is  unfortunate, 
and  a  means  to  overcome  this  deficiency  has  yet  to  be  devised.  A  means 
to  effect  the  replacement  of  the  higher  energy  pulses  that  were  removed 
due  to  their  specific  orientations  within  the  counting  tube  has  been 
developed  and  presented,  earlier  in  this  report.  The  probability  function 
(P  )  can  be  directly  applied  to  the  raw  spectrum  from  the  multichannel 
analyzer  to  replace  "rise  time  removed"  pulses  in  this  situation. 

The  other  portion  of  the  problem  in  the  determination  of  a  reconstructed 
neutron  spectrum  has  to  do  with  the  removal  of  the  wall  and  end  effects. 
As  mentioned  earlier,  this  problem  was  essentially  solved  by  Brown  and 
extended  by  Wang.  Modification  to  Wang's  probability  function  have  been 
made,  and  the  revised  probability  function  P(E-»E.  )  has  been  previously 
presented. 

What  remains  to  be  accomplished  is  to  present  a  computer  program 
outline  of  the  means  to  unfold  the  various  effects  and  to  arrive  at  a 
complete  neutron  spectrum. 

For  the  purposes  of  this  section  and  to  add  clarity  to  the  computer 
program,  it  is  considered  that  the  following  subroutines  have  already 
been  accomplished; 
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(a)  Assume  that  the  raw  spectrum  has  been  corrected  for  rise  time 
discrimination.   The  rav  counts  in  each  (E  +  Q)  channel  have  been 
divided  by  the  probability  function  (l-P  )  and  the  revised  channel  counts 
arc  now  designated  C   . 

(b)  The  counter  efficiency  has  been  calculated,  and  an  energy  dependent 
value  of  the  efficiency  can  be  shown  in  the  various  reaction  cross - 
sections.   The  efficiency  of  the  counter  was  developed  by  Wang  to  be  a 
function  of  (a)  the  nuclear  cross  sections,  (b )  the  counter  dimensions, 
and  (c)  the  gas  pressures. 

where 

KL   =  number  of  He  atoms/cm  of  sensitive  volume  at  filling 
^He 

pressure 

c—  (E,T)  =  He  (n,p)T  cross  section  (see  Fig.  2) 
n,p^  It       v  '  '  v         ' 

D  =  Effective  thickness  of  the  detector  for  4jt  geometry 

=  Total  sensitive  volume 
effective  area 

where  sensitive  volume  of  sphere  =  sensitive  volume  of  cylinder 

■■  '      k         "b  2 

~  rt  IT  =  «  o     L 

b  -  (|  -o2l:)V3 

thus  effective  area  =  %   R  =  '*(£  b2L)2' 


D  =  7TT2-,?/-  *  (t)   *    l 

rt(-^  o  L) 


•l&&h       ^ 


f . (E  )  =  Probability  the  proton  or  triton  will  not  hit  the  wall 

J.  i\ 

A  similar  expression  can  be  written  for  the  efficiency  of  the  He 
(n,d)D  reaction 

1   n,d<V  "  \     °a,d  (EP)D  f5  <V  I 

ne 
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vhere  the  notation  is  similar  to  the  above 

It  can  be  seen  that  the  efficiency  can  he  reflected  in  a  new  cross- 
section  notation  such  that: 

Y)  {■£    ) 

o-   (E  )*  =  N   Do-   (E  )  =     n>v      ,N 
n,p  "V  -Vn^NJ      fl  ^V 

vhere  o"   (E .  )-*  is  the  cross  section  that  has  been  corrected  for  the 

efficiency  of  the  counter. 

and 

c    ,(Ej*  =  Bl   DC  M   (eJ  =   ^A^n 

vhere,  o"  .(E,_)-*  is  the  cross  section  of  the  5TT  (n,d)D  reaction  that  has 
n,dv  W  ^Kev  ;  ; 

been  counted  for  the  physical  efficiency  of  the  counter. 

'Assume  nov  that  the  higher  energy  rise  time  discrimination  removal 

counts  have  been  replaced,  and  the  cross  sections  have  been  corrected 

for  the  counter  efficiency  as  indicated  above.  A  computer  program  to 

unfold  the  vail  effects  and  to  calculate  the  resultant  true  neutron  ■ 

spectrum  could  be  as  follovs : 


(1)  Let  E.t  -B  -^ 


P 


vhere  E . ,  =  Actual  neutron  energy 
J 

E .  = "Registered  neutron  energy 

.  J 

(2)  Let  A.t  =  Number'  of  neutrons  of  energy  E.t 

E.t  ^  E. 
0  r     J 


(3)  Let  C--  =  Number  of  counts  registered  in  top  energy  channel  of 
the  multichannel  analyzer  (after  correction  for  rise- 
discrimination) 

(5)  Solve  for  A.N, 

V   .  . c™- 
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^  CR(»-D  -  V-D'  <rn,p<IW*0!H-l)  +  AN'  WV*  * 

•where  P_  (E  t-^E/.t  -  x )  is  the  probability  that  an   (n,p)  event  from 
P       I\  \.A-iJ  -i    _     ]_ 

channel  E„  will  be  registered  in  channel  E/..  ,  \         P(E..-- »*E.  )  =  2=r 


*L 


(7)  Solve  for  A^^, 

(8)  Repeat  this  procedure  until  ve  reach  channel  A,      ,.  x, 

where: 

1*3-1 

°"n,p  CVjW*fl  (E(N-jp 

(9)  Let  E,T  • .  -   =  E„T  -   (Q         -Q     , )  and  the  deuterons  begin  to  appear. 


i.  e.  E       .    .        .  -  E„  -  (0.76^  +  3.27)  =  EATt   -  3.27 
registered         N       v  '         N 


i=J 


(10)  VM)=A(H-5-l)'°i/(»-M)',Ml(VM))+      2j    A(N"°' 

i=0 

'     ^(u-ip*  t"fi  (E(N-i)3'pp  (E(i^T  Vj-1)5     +  ak'  "i/V)* 

fj.  (E/-,  »  i  \)-  The  last  tern  in  this  expression  indicates  the 

deuteron  effects  from  ErT,  and  f,_(E/,T  .  .  \)  =  the  probability  that  the 

N'      5  (N-j-1)7 

deuterons  of  energy  E/„  .  ,  \  vill  not  hit  the  wall. 

(11)  Solve  for  A(K_._l)t 
(N-J-l)'«  =  CR(N-o-l)  "       A(N-i)'  ^p^Oi)''*  L^l^N-i)^ 


pp  (E(¥^TE(n-j-i)) '-  V  °"n,d<V*  f5  (Vj-i)' 

divided  "oy   V(Vj-l)')*fl(E(K-j-l)) 
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(lwd)   Then  in  general: 


i  =  g  +  k 

E 

i  =  0 


(N-i)'  n,p^  (N-i)w 


^i  (Vi))J,rP  {E(i^TE(;.-j-i-k))  -  A(N-k)'  °-n,d  <V*)->" 

i  =  k  -  1 


f5  (E(N-^-l-k))  ~ 


V'i)'  ^d^N-i)'^ 


f_6   (E(N-j-l-k)^l 


i  =  0 


Pd  (E(N-j-l-i)  E(N-o-l-k)) 
divided  by    ^^■(E(K.j.1_k),)^  V^Wl-k^ 

uhere  P,   (E/.T  .  - — rt"E/.T  .  .  ,  0   is  the  probability  of  a  deuteron 
d  v  (N-j-l-i)  (N-j-l-k)'  •        *  J 

having  scattering  from  a  higher  energy  channel  into  channel  E/„  '.  -  ,  \ 

and  f^  (E/,T  .,,<,)  =  the  probability  that  the  deuterons  having  energy 

E/„  .  ,  .  \  vill  intersect  the  vail  =  l-f_  (E/.T  . '  .  ,  \) 
(Bi-j.-1-k)  •  [_   5   (N-j-l-k)^ 


(15)  This  procedure  is  followed  until  the  spectrum  is  unfolded  to  the 
love st  energy  channel  that  has  escaped  rise  time  discrimination. 

(Ik)     The  unfolded  spectrum 

...A, 


-AN"  A(N-1)"  A(N-2)W  A(N-3)' 


•(N-j-l-k) 


is  then  the  neutron  spectrum  corrected  for  (a)  rise  time  discrimination, 
(b)  vail  effects,  and  (c)  counter  efficiency. 


-74. 


E..    Random  Probability  Inversion  Matrix 
The  foregoing  sections  have  considered  in  detail  the  computation's 
necessary  to  change  the  raw  energy  spectra  to  a  neutron  spectra  by  the 
removal  of  the  rise  time  discrimination  effects,    the  wall  and  end  effects, 
and  corrections  for  the  spectrometer  efficiency. 

Unfortunately,    as  can  also  be  seen  from  the  foregoing  sections,    the 
number  of  computations  necessary  to  remove  these  various  effects  and 
corrections  is  exceedingly  large.     There  are  four  position  variables,    two 
independent  energy  variables  (one  reaction  particle  variable,    one  incoming 
neutron  variable),    risetime  discriminator  setting,    and  an  applied  voltage 
condition- -a  total  of  8  variables.     The  number  of  variables  and  constraints 
on  individual  computations  make  their   solution  by  analytic  means  virually 

s 

an  impossibility,    and  even  when  considered  for  high  speed  computing 
machines,    the  number  of  computations  involved  would  require  several 
hours  of  machine  time.     It  was  therefore  decided  that  in  order  to  compute 
an  "inversion  matrix'-'  which  could  be  applied  to  the  ."raw  spectra"  output 
of  the  multichannel  analyzer  recourse  must  be  made  to  a  "random  prob- 
ability method  being  that  it  will  require  only  a  comparatively  small 
amount  of  electronic  computing  machine  time,    and  in  addition  will  permit 
computation  of  all  corrections  to  the  raw  spectra  simultaneously,    as  well 
as  elimination  of  several  assumptions  that  previously  were  necessary  to 
apply  to  the  more  general  case. 

The  random  probability  analysis  for  our  counter  has  been  developed 
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by  means  of  a  FORTRAN  PROGRAM  having  the  code  name  "CATMAN.  " 
An  outline  of  the  random  probability  analysis  is  presented  below: 

(1)  A  neutron  of  a  specific  energy  E-^-  (selected  as  the  maximum  neutron 
energy  that  our  counter  is  expected  to  detect)  comes  into  the  counting 
tube.      The  direction  of  the  incoming  neutron  is   selected  at  random,    and 
since  the  detector  is  to  have  a  4  7f  response,    this  direction  may  be 
completely  arbitrary. 

(2)  The  neutron  passes  into  the  volume  of  the  tube  (not  necessary  to  limit 
this  to  the  "active  volume"  between  the  field  tubes)  and  transits  a  distance 
into  the  gas  medium.     Since  each  unit  volume  is  presumed  to  have  an  equal 
probability  of  a  reaction  occurring  in  that  volume,    the  position  of  the  event 
is   selected  at  random  by  a  "random  number  generator"  which  determines 
both  longitudinal  distance  along  the  anode  wire,    and  radial  position  Tq. 

(3)  An  event  occurs  at  a  radius  Yq  from  thepcounter  center  line  and  at 

some  length  X  along  the  counter  length.     The  event  can  be  either  a  "^rIe(n,p)T 

3 
reaction,    or  a     He(n,.d)D  reaction.     The  probability  of  the  type  of  reaction 

is  based  on  the  ratio  of  the  cross  sections  at  the  energy  of  the  incoming 

neutron  to  the  sum  of  their  cross  sections,    i.  e.   the  probability  of  a 

He(n,  p)T  reaction  is: 


crn.D(EN) 


n,  p 


^p^N)  +  °~n,d(EN) 
and  the  probability  of^a.    He(n,  d)D  reaction  is: 

o-n,  p(%)  +  o-n>  d<EN>       ' 
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(4)    The  reaction  particles  are  assumed  to  have  isotropic  distribution  in 

the  center  of  mass  system  of  coordinates,    and  the  direction  of  one  of  the 

reaction  particles  is  selected  at  random  (three  dimensional  selection 

based  on  0        ,    cos  0        ,    and  equal  but  opposite  momentum  vectors).     The 
cm  cm 

energy  of  the  particles  in  the  CM  system  is  divided  among  the  particles 
as  follows: 


EP=m      =f<iEN  +  Qn,p> 


1    /3 


■3 

/      He(n,  p)  reaction 


T  cm  4    v4      in  n,  p' 


TT 


4     4 

1    ,  3 


Dlcm    =ED2  cm=  j(4EN  +  Qn,d)    (      He (n,  d)D  reaction 


(5)  The  reaction  is  converted  into  the  Laboratory  system  of  coordinates, 
and  the  energy  of  the  various  particles  is  computed  in  the  Laboratory 
system.     Note  that  any  energy  between  a  maximum  and  a  minimum 
Laboratory  energy  is  equally  probable  and  can  occur  for  the  proton  and 

Q  I 

triton  particles.     Likewise,    if  —  E-!Sj>    Qn  M   (threshold  energy),    the 
two  deuteron  particles  will  vary  between  a  maximum  and  minimum 
Laboratory  energy  with  equal  probabilities.     (From  Appendix  C.  )    The 
Laboratory  energies  are  determined  by  the  random  selection  of  0Cm  an<^ 

cos  ecm- 

(6)  The  length  of  track  of  the  various  Laboratory  particle  energies  are 

computed  by  the  relationships  previously  presented. 

(7)  The  track  lengths  and  orientations  are  then  compared  to  the  reaction 
site  to  determine  if  any  portion  of  the  tracks  has  intercepted  the  walls 
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of  the  tube,    and/or  has  a  portion  of  its  length  in  the  gas  volume  determined 
by  the  perpendicular  radial  projections  of  the  interior  ends  of  .the  field 
tubes  and  the  end  walls  (end  effects). 

(8)  If  the  track  of  one  of  the  particles  intersects  the  wall  or  the  end 
portion  of  the  counter,    the  portion  of  the  energy  represented  by  the  track 
that  is  physically  within  the  active  volume  is  calculated.     Note  that  this 
energy  is  not  directly  proportional  to  the  track  length  since  the  energy 
deposition  is  not  linear  (see  Fig.    5).      The  energy  deposited  in  the  active 
volume  is  equal  to  the  total  energy  of  the  particle  minus  the  energy 
represented  by  the  length  of  track  that  passes  outside  the  active  volume 
(if  the  same  gas  mixture  were  present  outside  the  container).     This 
conclusion  can  be  drawn  from  the  concept  that  after  a  charged  particle 

has  lost  a  certain  amount  of  energy,    it  still  has  the  same  length  of  path 

f 

to  travel  as  an  identical  particle  that  begins  its  life  with  an  energy  equal 

to  the  reduced  energy  of  the  particle  in  question. 

(9)  The  various    rn    2    3    4)  va^ues  (radial  components  of  track  lengths) 
are  calculated  to  determine  rise  times,   with  the  condition  that  an  appro- 
priate constraint  be  placed  on  those  tracks  that  intersect  the  walls  of  the 
tube,    or  pass  into  the  end  regions. 

(10)    The  rise  time  of  the  pulse  is  then  determined,    based  on  the  time  that 
it  takes  for  the  nearest  (rn)  and  most  distant  electron  (r.f)  to  reach  the  , 
anode. 


1  P  o 


2        2 

r      -  r 
f         n 
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(11)  A  weighting  frunction  is  applied  to  the  result  to  account  for  the 

greater  volume  probability  at  points  more  distant  from  the  anode. 

g(r)dr  =-|   dr 
b^ 

where    r    is  the  reaction  site   selected  in  (3)  above. 

(12)  The  weighted  pulse  is  catalogued  by  energy  (E-^  +  Q),    rise  time, 
and  reduced  energy  (if  applicable). 

(13)  Steps  1  through  12  are  repeated  several  thousand  times  for  the  same 
incoming  neutron  energy.     The  probability  of  the  various  particles  not 
hitting  the  wall  is  computed.     £-,  (E),!  the  probability  of  the  proton-triton 
tracks  not  hitting  the  wall,    and  fc(E),    the  probability  of  the  deuteron 
combinations  not  hitting  the  wall,    are  calculated  after  "weighting"  for 
volume  probability  considerations.     The  probability  of  causing  a  pulse  in 
a  lower  energy  channel  P(E  •— *  E- )  is  likewise  calculated  from  the  track 
length-energy  relationship  described  in  step  8,    after  the  volume  prob- 
ability function  has  been  applied.     All  three  of  these  probability  functions 
are  combined  into  a  new  probability  function   P    (E  — ^-E^)  for  the   spectrum 
unfolding  procedures. 

(14)  The  energy  of  the  incoming  neutron  is  reduced  by  one  energy  unit, 
and  steps  1  through  13  repeated. 

(15)  The  energy  of  the  incoming  neutron  is  reduced  in  successive   steps' 
until  the  lowest. value  of  interest  is  obtained,    or  until  the  length  of  track 
of  the  charged  particles  are  less  than  the  length  required  to  escape  rise 
time  discrimination. 
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(16)     The  weighted  number  of  counts  in  their  respective  energy  channels 
are  totaled  and  accounted  for,    as  explained  in  the  preceding  section.     A 
matrix  arrangement  can  then  be  compiled  and  corrected  for  the  efficiency 
of  the  counter.      The  matrix  arrangement  could  appear  as: 


n 


•11        a12        a13' 


■21        ^22       ^23' 


lln 


'2n 


"nl 


an2        an3' 


nn 


*1 


x 


A, 


A 


A 

n 

Total  neutrons 
sent  into  each 
energy  range. 


Counts  regis-  a--  =  fraction  of  neutrons  of  energy 

tered  in  each  j  that  are  counted  at  energy   i. 

energy  incre- 
ment for  random 
probability 
analysis. 


(17)     The  problem  now  resolves  itself  into  determining  an  inversion  matrix 
whereby  we  can  calculate  the  original  neutron  spectra  from  the  counts 
registered  in  each  energy  channel  of  the  multichannel  analyzer,    i.  e.   to 
•compute  an  inversion  matrix  such  that: 
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A, 


A. 


A 


N 


Neutron 

Spectra 


bll        b12       b13' 


b21        b22        b23' 


bnl        bn2        bn3' 


Inversion 

Matrix 


1 


'In 


2n 


b 
nn 


C. 


C(N-1) 

CN'| 


Raw  Counts 
from 

Multichannel 
Analyzer 


(18)  Note. that  in  steps  16  and  17  we  have  referred  to  a  matrix  arrange- 
ment and  an  inversion  matrix  for  our  illustration.     The  actual  matrices, 
if  they  were  to  be  constructed,    would  be  an  array  of  500x500  or  more, 
and  would  be  difficult,    if  not  impossible,    to  invert.     Instead  of  actually 
constructing  these  matrices,   the  computer  would  solve  this  problem  by 
a  subtraction-division  process  similar  to  that  explained  in  the  preceding 
section  II  (D). 

(19)  The  complete  random  probability  analysis  in  FORTRAN  language 
(with  appropriate  comment  cards)  is  included  in  Appendix  E. 
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III.  CONCLUSIONS  AND  PRESENT  STATUS  OF  PROJECT 

3 
At  the  present  time,    the     He  neutron  spectrometer  described  herein 

as  under  development,    has  not  been  tested.      The  electronic  circuitry  for 

the  rise  time  discrimination  however  has  been  constructed  and  operated, 

using  the  small  JHe  tube  used  by  Wang.     This  tube  is  filled  with  10  atmoS' 

3 
pheres  of     He  (no  heavy  stopping  gas),    has  only  a  15/16  in.    inside 

diameter,    and  an  energy  width  E-^  of  only  0.  5  MeV  and  an  energy  length 

Et    of  2.  5  MeV       for  protons.  '     The  intrinsic  tube  efficiency  (due  to  the 

small  size)  however  is  extremely  low,    even  at  these  moderate  energies, 

and  it  was  not  considered  that  the  expense  of  developing  an  inversion 

matrix  was  warranted. 

The  "inversion  matrix"  for  the  larger  tube  that  is  on  order  from  the 
Texas  Nuclear  Corporation,    has  essentially  been  completed,    and  the 
"random  probability  program"  (CATMAN)  which  is  used  to  generate  the 
"inversion  matrix"  has  been  completed. 

Future  work  in  the  development  of  this     He  spectrometer  could  lie 

in  some  means  of  electronically  comparing  the  rise  time  of  pulses  to  the 

pulse  magnitude  (energy)  to -determine  if  the  pulse   should  be  "gated"  into 

the  multichannel  analyzer.     This  comparison  feature  would  be  in  lieu  of 

a  set  value  (t    )  for  the  rise  time  discriminator  circuitry,    and  would 
r 

enable  the  registering  of  those  valid  low  energy  pulses  which  are  now 


'"These  figures  have  been  recalculated  in  this  report,    since  Wang  made 
an  error  in  the  density  of  JHe  on  pages  87  et  seq.    of  Reference  16. 
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automatically  eliminated  because  their  short  track  lengths  cannot  over- 
come the  rise  time  discriminator  setting.     The  rise  time   setting  would 
then  be  some  function  of  energy,     i.e.    t     (E).     This  improvement  would 
enable  the   spectrometer -to  still  differentiate  between  valid  pulses  and 
recoils  without  sacrificing  all  low  energy  valid  pulses  in  the  process.     It 
seems  somewhat  paradoxical  that  efforts  should  be  made  to  extend  the 
range  to  energies  in  the  neighborhood  of  1.  0  MeV  from  the  upper  energies 
when  the  original  work  of  Batchelor^  was  aimed  at  extending  the  range  of 
his  spectrometer  beyond  1.0  MeV  from  the  lower  energies  by  eliminating 
the     He  recoils. 

A  subsequent  report  ^Part-Id^  will  be  written  by  the  Health  Physics 
Group  of  the  Lawrence  Radiation  Laboratory,    Berkeley  concerning  the 
actual  functioning  of  the  ^He  neutron  spectrometer,    and  departures  (if 
any)  from  the  predicted  operational  characteristics. 
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Appendix  A-  DESIGN  OR  EXPERIMENTAL  APPARATUS 

3 

1.   Design  of  the  Pie  proportional  counter  tube 

a.   Previously  constructed  He  proportional  counters. 

The  design  of  the  helium-3  proportional  counter  tube  was  .completed 

after  an  extensive  search  of  the  technical  literature.   To  aid  in  the 

selection  of  the  counter  tube  parameters,  a  summary  of  comparative  sizes , 

3 
gas  compositions,  filling  pressures,  and  other  pertinent  data  on  all  He 

proportional  counters  and  ionization  chambers  constructed  to  date  was 

tabulated  (Table  Al).   The  data  necessary  to  complete  Table  I  is  not 

reported,  but  it  is  apparent  that  a.  variety  of  heavy  monotonic  gases  and 

stabilizing  gases  have  been  employed  by  experimenters.   Counter  dimensions 

and  shapes  have  likewise  been  somewhat  varied. 

"b.   Lawrence  Radiation  Laboratory  usage  requirements 

3 

inis  Ke  proportional  counter  will  ultimately  be  used  for  radiation 

surveys  by  the  Health  Physics  Department  at  the  Lawrence  Radiation  Lab- 
oratory. The  objective  of  this  counter  therefore,  is  to  be  capable  of 
detecting  and  accurately  evaluating  neutron  spectra  of  energies  to  about 
14.0  MeV  regardless  of  the  neutron  source  direction.  It  is  likewise 
necessary  that  the  apparatus  be  capable  of  being  readily  transported  and 
placed  in  relatively  inaccessible  locations.  The  usage  requirements  thus 
consist  essentially  of:  (l)  hit   resolution,  (2)  relatively  small  size,  (3) 
the  high  voltage  requirements  remain  within  the  capacities  of  existing 
Laboratory  equipment,  i.e.  less  than  6000  volts  (an  additional  considera- 
tion that  was  imposed). 

c.  Design  variables 

In  the  design  of  a  proportional  counter,  the  designer  has  essentially 
five  variables  that  are  available  to  him,  viz:  shape  (geometry),  dimensions 
of  the  container  (cathode),  dimensions  of  the  anode,  gas  composition  and 
pressure,  and  applied  voltage.  It  is  to  be  noted  .that  once  a  design  has 
been  fixed,  the  only  variable  available  to  the  experimenter  is  the  applied 
voltage.   Unfortunately,  the  five  variables  each  effect  the  final  prop- 
erties of  the  counter  and  are  essentially  inter-related. 

d.  Design  parameters 


TABLE  A-l 
Physical  Characteristic*  of  3He  Proportional  Counters  and  Ionization  Chamber*  Reported  In  the.  Technical   Literatur 


Type 

Counter 

counter 

ahape 

Prop. 

Cyl. 

Prop. 

Cyl. 

Prop. 

Cyl. 

Prop. 

Cyl. 

Prop. 

Cyl. 

r     Cathode  Anodo  Effcc-  Opera-  Gas  Energy 

radius  radiua       tive          ting  multi-  range 

(in.)  (in.)  length  voltage  plica-  (MeV) 

(in.)  (volt.)  tion  


l-Mhn;   Caa    (a 


Stabilising  gaa 


CO,  Cll. 


0-1.0  .071 

1.0071) 


1.0625      .002  4.70       2600      ~10  0-1.2       0.355 

(0.355) 


1.965        .0025  5.51        6000- 

8000 


0-5.0      (l.Occ)  6-10 

(SSTP 


Prop. 

Prop. 


Cyl. 

Cyl. 


1.0625      .002  4.79       4400       -10  0-2.5        1.355 

(1.355) 


0.9375      .002  9.00 


0.980        .001  6.00       2450  4-5  7.0  (1.09) 

0.480        .001  4-12       ~1200  (1-10) 


Prop. 

Cyl. 

0.  480 

.001 

i.O      -1900 

Prop. 

Cyl. 

0.480 

.003 

2100 

0.480 

.003 

920 

0.480 

.001 

1050 

0. 3125 

.006 

5000 

0. 3125 

.005 

5000 

Prop.        Cyl.  2.00  .004  15.0 


(2) 


(10) 

(2) 
(6) 

(9.7) 
(9.7) 


.0066 


-17.  5   (0.976) 

4.0 

.0418 

(0.013)5) 

4.25 

.  0445 

(0.992) 

4.0 

.0418 

(0.0501) 

4.25 

.  0445 

0.3 
0.3 


Expcr.mcn 

Date. 
Total  and 

Preaaure       Reference 


R.    Batchclor 
(1952)  Rel.  1 


Datchclor,    Avca.  and 
S'*yrme(1955)Ref.  3 


Bloom,   Reilly,  and 
Topple  (1955)  Ret.  2 


Batchclor  and 
Morriaon(1960)Ret.  4 


5.02 

A.  R.  Sayre ■ 

4.  32 

(1960)  Ret.   S 

5.03 

4.  34 

5.60 

W.K.    Brown 

(1962)  Re/.  8,9 

Mille,  Caldwell,  and 
Morgan(1962)  Ret.  25. 
Wang(1963)  Ret.    16 

Tcxaa  Nuclear  Corp. 
(19641  Ret.   26 

J.  L.  Fricdea  and  R.  E. 
Chrien  (1964)  Ret.  11 


Sayre •  and  Coppola 
(1964)  Ret.    10 
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TABLE  A-l   (com.  ) 


Type       Counter     Cathode     Anode      Effcc-    Opera-  Gas 

counter      shape  radiu.     radiu.        tive         ting  multi- 

(in.)  (in.)       length   voltago  plica- 

(>n.  )      (volt.)  tion 


Prop. 

Cyl. 

0.4375 

.0025 

10.0 

4000 

4  50 

0.4375 

.0025 

10.0 

2500 

11 

0.4375 

.0025 

10.0 

5000 

93 

0.4375 

.0025 

10.0 

4000 

42 

0.4375 

.002  5 

10.0 

2400 

32 

0.4375 

.0025 

10.0 

750 

1.0 

0.4375 

.0025 

10.0 

1000 

1.2 

0.4375 

.0025 

10.0 

2000 

2.4 

0.4375 

.0025 

10.0 

3000 

8 

0.4375 

.0025 

10.0 

4200 

15 

0.4375 

.0025 

10.0 

4500 

390 

0.4375 

.001 

6.0 

1000 

100 

Ion 

Cyl.      ~  1.375  ~ 

).  14 

9.25 

6000 

Ion 

Sphere 

3.94 

1.295 

N.A. 

Ion 

Sphere 

0.  511 

.0197 

N.A. 

>1600 

Energy 

r 

llinc  Gas  (at 

mospne 

c 

1 .», 

STP 

Experimenter, 

range 

He 

Stc 

PP 

»e  B« 

Slabil 

«a. 

Date, 

and 

(3He) 

A 

Kr          Xe 

co2 

0.61 

CH4 

N2 

Pre. sure 

Reference 

5.  5 
(7) 

6.  11 

O.    K.    Harllng 

(1965)  Ret.    19 

5.  5 

(7) 

0.611 

6.  11 

9.02 

(7) 

1.00 

10.02 

1  1 .  64 
(?) 

0.61 

12.25 

5.805 

(7) 

0.  305 

6.  11 

5.93 

(7) 

£&' 

6.59 

5.93 
(7) 

(DF») 
0.  06 

6.59 

5.93 
(7) 

(OF,) 
0.66 

6.59 

5.93 

IBFj) 

6.  59 

(7) 

0.66 

5.93 

(7) 

(BF,) 
0.66 

6.59 

8.71 
(7) 

0 

97 

9.68 

5.92 

(7) 

7 

7 

7 

5.92 

2.7 
(2.7) 

3.8 

0 

07 

6.  57 

Freeman  and  Weit 

(1962)  Re*.   7 

3.  1-1.0 

0.46 
(0.46) 

2 

5 

0.  118 

J. 08 

N.  P.    Olailtov 

(1957)  Rev.    13 

J.  2-0.  8 

2.5-4 

2.  S 

-4 

0.1-0. 16 

5.1-8+ 

A.    I.    Abramov 
(1958)  Rev.    12 
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In  addition  to  the  three  usage  requirements,  there  exist  several 
mechanical  and  electrical  requirements  that  the  counter  must  fulfill. 
These  requirements  are  : 

(1-)  To  fulfill  the  Ifcc  requirement,  the  filling  gas  must  he  capable  of 
completely  stopping  maximum  energy  product  particles  within  the  active 
volume  of  the  container. 

(2-)  The  wall  thickness  of  the  container  must  "be  sufficient  to  contain  the 
filling  gases  under  pressure. 

(3)  The  gas  fillings  must  not  provide  competing  reactions  that  can  be 
mistaken  for  bona -fide  events.   Recombination  should  be  minimal- 

(1+.)  Gas  multiplication  to  some  degree  is  necessary  to  overcome  electrical 
noise,  and  to  obtain  independence  from  the  reaction  site  within  the  counted 

(5.)  "Jitter-time"  (the  time  it  takes  for  an  electron  at  a  point  on  the 
cathode  to  migrate  to  the  anode)  must  be  relatively  small,  and  within  the 
resolving  capabilities  of  the  associated  electronics. 

(6-)  Electric  field  distortion  must  be  minimized,  and  field  gradient  at 
equal  distances  from  the  anode  should  be  equal. 

(•7.)  Electrical  breakdown  and  "Corona  discharge"  must  be  eliminated. 

e.   Effect  of  design  variables  on  design  parameters 

The  changing  of -any  design  variable  has  a  profound  effect  on  the 
design  parameters.  The  interdependence  of  the  parameters  to  several 
variables  makes  any  proportional  counter  a  compromise.  The  changing  of 
the  variables  has  an  effect  upon  the  parameters  as  follows; 

(1-)  Shape  (geometry) — effects:  kit   response;  field  gradient  distribution 
(thus  "Jitter  time");  space  charge  distribution. 

(2-)  Dimensions  of  the  cathode — effects:  amount,  pressure,  and  composition 
of  the  filling  gas;  field  gradients  (thus  "jitter  time"). 

(3.)  Dimension  of  the  anode — effects:  field  gradient  ("jitter  time"  and 
gas  amplification);  corona  discharge. 

(k)  Gas  composition  and  pressure — effects:  stopping  power;  jitter  time; 
gas  amplification;  competing  nuclear  events;  relative  cross  sections  of 
the  composite  gas. 
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(5)  Applied  voltage --effects:  field  gradient  (thus  Jitter  tine  and  gas 
multiplication);  corona  effect. 

f •   Selection  of  the  appropriate  design  variables 

(1)  The  geometry  selected  was  a  cylindrical  proportional  counter.  A 
cylinder  was  selected  instead  of  a  sphere,  due  to  the  adverse  effect  of 

a  concentration  of  positive  space  charges  near  the  anode  after  the  occurr- 
ence of  gas  multiplication.  .While  this  condition  likewise  exits  in  a 
cylindrical  proportional  counter,  the  space  charge  effects  are  distributed 
over  the  length  of  the  anode  wire  and  are  not  concentrated  at  a  single 
point  as  in  a  spherical  configuration.   It  was  also  considered  that  the 
resolving  tine  lag  would  be  less  (thus  the  counting  rate  higher)  and 
multiplication  more  uniform  in  a  cylindrical  chamber.  This  selection, 
however,  necessitated  some  sacrifice  of  the  4it  response  of  the  proportional 
counter--a  deficiency  that  was  partially  overcome  by  increasing  the  gas 
pressure.  The  field  gradient  distribution  is  radial  along  a  central  axis 
in  a  cylindrical  tube,  except  near  the  ends,  instead  of  a  direct  radial 
distribution  as  in  a  spherical  chamber.  The  end  distortion  of  the  field 

gradient  of  the  cylinder  can  be  eliminated,  however,  through  the  use  of 

27 

appropriately  designed  field  tubes. 

(2)  The  dimensions  of  the  cathode  were  arbitrarily  selected  to  be  a  maxi- 
mum outside  radius  of  2.0  in.  and  an  active  length  of  15-0  in.  These 
dimensions  were  considered  to  be  about  the  maximum  that  would  meet  both 
the  portability  requirement  and  still  maintain  a  reasonable  field  gradient 
(thus  jitter  time)  near  the  outer  portions  of  the  tube's  radius.  It  is  to 
be  noted  that  the  tube  size,  which  has  been  dictated  by  practical  considera- 
tion likewise  limits  selections  of  the  amounts,  pressures,  and  composition 
of  the  filling  gases.  The  thickness  of  the  cathode  walls  is  determined 

by  the  filling  pressures. 

(3)  Gas  composition  and  filling  pressures 

(3« a.)  The  gas  composition  must  be, such  that: 

(3-a.l)  The  reaction  particles  from  a  lij-.O  MeV  incident  neutron  will  be 
stopped  in  a  diameter  of  the  cylinder.  This  will  to  a  certain  measure 
restore  a  portion  of  the  k-x   geometry  resolution. 

(3- a. 2)  The  electron  migration  through  the  gas  mixture  should  be  sufficiently 
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rapid  that  extremely  long  jitter-times  are  not  encountered. 

(3- a. 3)  Competing  nuclear  events  are  minimized. 

(3«a.^)  A  stabilizing  gas  be  employed  to  absorb  photon  produced  from 
bremsstrahlung  during  gas  multiplication  and  to  prevent  spurious  pulses. 

(3. a. 5)  The  gas  multiplication  (which  is  a  function  of  the  pressure)  be 
sufficient  to  provide  a  pulse  that  is  both  proportional  to  the  energy  of 
the  incident  neutron  and  sufficient  to  mask  electrical  noise.  A  gas 
multiplication  of  at  least  10  is  considered  desirable.  U}  ] 

(3-b.)  The  selection  of  the  gas  composition  was  dictated  by  several 
factors . 

(3«b.i)  The  cross  section  of  the  Ke(n,p)T  reaction  decreases  by  a  factor 

-1  o 

of  10  for  neutrons  at  0.1  MeV  to  11  MeV.    It  vas  decided  that  some 
compensation  for  this  decrease  in  cross  section  should  be  made  if  neutrons 
of  higher  energies  are  to  be  detected  in  meaningful  numbers.  This 
reduction  in  efficiency  at  higher  energies  is  likewise  compounded  by  the 
fact  that  since  the  higher  energy  product  particles  will  travel  further 
in  the  gas  medium,  their  chances  of  intersecting  a  wall  are  increased 

and  thus  the  probability  of  the  production  of  a  full  energy  count  is 

3 
decreased,  it  was  considered  that  two  atmospheres  of  ,He  would  compensate 

somewhat  for  this  loss  of  efficiency  at  higher  energies. 

(3-b. 2)  Eelium-3  is  a  gas  with  relatively  low  stopping  power  for  charged 
particles.  A  heavy  monotonic  gas  must  be  added  to  insure  that  the  charged 
particles  stop  within  the  diameter  of  the  cylinder.  Previous  experimenters 
have  used  argon,  krypton,  and  xenon  for  this  purpose.  Argon  was  rejected 
in  this  design  for  two  reasons:  First,  the  pressures  of  argon  necessary 

to  stop  a  charged  particle  are  greater  than  those  of  either  Kr  or  Xe  for 

33 
the  same  energy.  Secondly,  the  A(n,a)S   reaction  becomes  significant 

2  'SO 

above  about  2.2  MeV.  '    Pulses  due  to  this  reaction  can  be  removed  by 

rise  time  discrimination,  but  the  reaction  data  is  incomplete  and  at 
higher  energies  an  argon  (n,p)  reaction  occurs.   The  disintegration  cross 
section  at  a  single  energy  of  l^-.O  MeV  has  been  measured  for  argon, 
krypton,  and  xenon  and  found  to  greater  for  argon  than  either  Kr  and  Xe. 
(The  disintegration  particles  produced  are  not  however  identified  in 
reference  31- ) 

Xenon  was  likewise  rejected  due  to  the  uncertainty  of  the  type  of 
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dis integration  particles  produced,  and  absence  of  data  concerning  the 
disintegration  cross  section  at  intermediate  energies.  The  rejection 
of  xenon  was  difficult  since  the  heavier  of  the  noble  gases  would  have 
required  smaller  filling  pressures  than  either  argon  or  krypton. 

Krypton  was  thus  selected  to  be  the  "stopping  gas".  Unfortunately, 

krypton  likewise  has  a  disintegration  cross  section  (tentatively  identified 

32 
as  an  (n,p)  reaction  by  Shamu   )  which  produces  a  scattering  of  disintegra- 
tion particle  energies  over  the  entire  energy  distribution,  the  main  portion 

31  32 

lying  in  the  h-Q   MeV  region.   ;    No  attempt  will  be  made  to  remove  the 

protons  produced  from  krypton  disintegration  from  the  resultant  spectrum 
due  to  the  uncertainty  of  their  energy  distribution. 

Having  selected  a  helium -3,  krypton  mixture  it  was  necessary  to 
determine  the  appropriate  filling  pressures.  This  was  accomplished  by  a 
trial  and  error  method  of  computation  based  upon  determining  the  length 
of  track  of  a  proton  having  the  energy  of  the  average  cosine  of  scattering 
angle  (Lab  System )from  a  1^.0  MeV  incident  neutron.  When  the  length  of 
track  was  equal  to  the  radius  of  the  cylinder,  it  was  considered  that  the 
pressure  requirement  had  been  met,  and  that  a  configuration  that  could 

approach  the  desired  kx   resolution  had. been  accomplished.  It  was  found 

3 

.that  10  atmospheres  of  krypton  plus  2  atmospheres  of  He  met  tnis 

condition . 

(3..b.3)  A  stabilizing  gas  was  considered  essential  and  0.5>>  CO  was 

selected  for  this  purpose.  Carbon  dioxide  was  selected  in  lieu  of  methane 

3  7 

which  has  been  successfully  utilized  by  Batchelor,   Freeman  and  west,. 

and  Sayres   in  their  most  recent  counters.  The  reason  for  this  change 

is  that  our  counter  is  designed  to  evaluate  neutron  spectra  to  lt.0  MeV 

(considerably  above  the  energies  evaluated  by  previous  experimenters), 

and  it  was  considered  that  the  introduction  of  a  source  of  additional 

recoil  protons  (from  the  methane)  into  the  tube  would  provide  unnecessary 

complications --especially  in  view  of  the  additional  proton  background  due 

to  the  krypton  disintegrations  (which  was  not  a  problem  at  lower  neutron 

energies). 
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(3«bA)  For  comparison  purposes,  the  cross -sections  of  the  various  gas 
reactions  are  plotted  in  Fig.  Al  as  a  function  of  neutron  energy.  The 
cross  section  for  the  C0o  was  obtained  by  adding  the  energy  dependent 
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Fig.     Al  .      Neutron  Reaction  Cross  Sections 


Source  of  Data 

(1)  Calculated  from  6~rp(E)  carbon  +  cr_(E)  oxygen  (Ref.    33' 

(2)  Ref.    34  L  L  3o 
(3-7)     Ref.    18 

(8)  Ref.    32 

(9)  Ref.    18 
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total  cross  section  of  carbon  to  twice  the  total  cross  section  of  oxygen. 
In  addition,  a  relative  plot  of  the  cross- sections  is  presented  to  displ 
the  collision  probabilities  ,!seen"  by  a  neutron  as  i_o  enters  the  tube 
(Fig.  A2).   It  should  be  noted  that  the  disintegration  cross  section  of 
the  krypton  at  10  atmospheres  pressure  becomes  a  factor  of  major  con- 
sideration at  higher  energies,  and  probably  will  determine  the  upper 
energy  limit  of  further  refinement  on  this  type  counter. 

(3-b.p)  The  dimensions  of  the  anode  'were  now  considered  and  several  applie 
voltage  and  wire  sise  combinations  "ere  attempted.  Jitter  time  "was       i 
calcul3.ted  by  use  of  a  "weighted"  electron  mobility  factor  (presented  on 
page  19  )  and  found  to  vary  from  18-22  ^.seconds  for  applied  voltages  of 
6000-5000  volts.   The  gas  multiplication  was  calculated  from  extrapolations 

q 

of  Brown's  data  and  found  to  be  approximately  10-15  for  the  above  range 
of  voltages,  using  a  3  mil  diameter  wire.  The  jitter  time  is  considerably 
longer  than  desired,  but  increasing  the  wire  size  (for  fixed  cathode 
dimension  and  filling  pressure)  decreases  the  gas  multiplication  and  dees 
not  appreciably  decrease  the  jitter  time  in  the  voltage  ranges  considered. 
The  jitter  time  is  relatively  insensitive  to  changes  in  applied  voltages, 
therefore  if  we  seek  to  reduce  the  jitter  time,  the  only  feasible  method 
•is  to  reduce  the  diameter  of  the  counter  and  increase  the  gas  filling 
pressure.  This  alternative  is  undesirable  since  an  increase  in  gas 
pressure  would  increase  the  relative  disintegration  cross  section  of  the 
krypton — a  problem  that  is  already  larger  than  negligible  proportions. 

(5-b.6)  \le   must  therefore  accept  the  fact  that  the  counter  will  be 
comparatively  slow,  and  design  the  electronics  accordingly. 
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2.  Design  of  the  Associated  Electronic  Components 


a.   Design  Parameters 


^ 


The  electronic  components  of  the  lie  proportional  counter  must  he 
designed  to  meet  several  requirements.  The  majority  of  these  require- 
ments are  common  to  proportional  counters  in  general,  however  there  exist 
several  unique  problems  associated  with  the  broad  range  of  rise  times 
of  the  pulses  (essentially  3-22  j^sec  is  our  region  of  interest).   The 
design  parameters  include: 

(1)  The  pulse  produced  by  an  ionizing  event  must  be  linearly  amplified 
and  be  capable  of  accurately  portraying  the  energy  of  the  ionizing  event. 

(2)  The  electronics  must  have  a  low  signal  to  noise  ratio. 

(3)  The  high  voltage  supply  must  be  well  stabilized,  since  variations 
in  the  high  voltage  will  cause  a  direct  variation  in  the  pulse  height. 

(k)   The  electronics  must  be  capable  of  distinguishing  between  a  desired 
event  and  an  undesired  event  (recoil). 

(5)  The  electronics  must  be  sensitive  to  the  relatively  large  range  of 
rise  times,  and  to  accomodate  pulses  that  vary  in  both  energy  and  rise 
time . 

(6)  Provision  must  be  made  to  permit  variation  in:  the  high  voltage 
supply,  the  counter  field  tube  settings,  the  rise,  time  discriminator 
setting,  and  the  time  constants  of  the  double  differentiating  network. 

b.  Electronic  design 

(1)  A  block  diagram  of  the  electronics  arrangement  constructed  to  meet 
the  above  parameters  is  indicated  below.  The  basic  arrangement  of 
components  is  similar  to  that  described  by  Sayres  and  Coppola   with  the 
exception  of  a  new  Peak  Sampler  Unit  and  a  Time  to  Height  Converter  Unit 
which  were  designed  and  constructed  by  Mr.  R.  M.  Brown  of  the  LPL  Physics 
Technical  Support  Group.  The  purpose  of  the  Peak  Sampler  Unit  is  to 
provide  a  gating  to  the  Multichannel  Analyzer  at  a  time  when  the  pulse 
has  reached  its  maximum  value . 

(2)  To  explain  the  operation  of  the  electronics,  consider  the  formation  of  a 
pulse  in.,  the  proportional  counter  as  the  result  of  a  collision  of  a 
neutron  with  one  of  the  contained  gas  molecules  (-)x:e,  Kr,  COp)«  The 
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collision  may  result  in  either  a  desired" event ,   i.   e.   disintegration  of 

3 
the     he  nucleus ,   or  an  undesired  evcnt--the  production  of  a  recoil. 

(2. a.  )  The  electrons  that  were   liberated  as   the  charged  particle   (s) 
transit  the   gas  are  accelerated  "by  the  electrical  field  gradient  and  a 
charge  is  deposited  on  the  anode.    .The  .size  of  the   charge  is   dependent 
upon  the  energy  of  the  particles  that  are  stepped  within  the   gas  medium, 
and  the  amount  of  gas  multiplication  in  the  area  surrounding  the  anode. 

(For  gas  multiplication  above  about  10,   the  energy  spectrum  is  undistorted 

29 

by  the  position  of  -che  ionizing  event  within  the  tube   providing,  that 

saturation  of  positive  ions  in  the  sheath  surrounding  the  anode  has  not 
been  achieved. )  The  electrical  pulse  begins  to  accumulate  when  the  first 
electron  reaches  the  anode,  and  reaches  a  maximum  when  the  last  electron 
finally  arrives  at  the  center  wire.   Depending  upon  the  type  of  event, 
the  initial  energy  of  the  incoming  neutron,  the  radial  distance  from, 
and  the  resulting  track  orientations  with  respect  to  the  anode,  this 
charge  build-up  period  may  vary  from  a  fraction  of  a  micro  second  to  as 
much  as  22  useconds  for  a  given  voltage  gradient  distribution  within  the 
counter.   The  actions  of  the  following  component  units  begin  when  the 
first  charge  begins  to  appear  on  the  anode,  and  continue  until  ail  the 
charge  has  been  collected. 

(2.b.)  The  incoming  pulse  is  amplified  through  a  charge  sensitive 
Preamplifier  and  an  Amplifier. 

(2.c.)  The  pulse  is  sent  into  two  parallel  channels  after  leaving  the 
Amplifier.   One  portion  of  'che  still  rising  pulse  is  sent  directly  to- 
the  Peak  Sampler  Unit  where  it  accumulates  to  a  maximum—at  which  time 
it  is  either  "gated"  into  the  Multichannel  Analyzer,  or  is  discarded. 

(2.d.  )  The  other  portion  of  the  pulse  is  sen'c  to  a  Double  Differentiating 
Circuit  (figure  A-k-)-  where  an  output  pulse  is  obtained  whose  duration  is 
dependent  upon  the  rise  time  of  the  input  pulse.   (As  an  approximation, 
the  output  pulse  from  this  unit  represent  the  "rate  of  rise  time"  since 
simple  PC  circuits  are  not  true  differentiators.)  The  time  between  base 
line  crossings  of  the  doubly  differentiated  pulse  is  thus  related  to  -che 
rise  time.    The  exact  relationship  between  the  differentiating  time 
constants  and  the  rise  time  is  not  critical  since  we  are  only  interested 
in  setting  a  discriminator  value  which  eliminates -all  pulses  having  a 
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rise  time  less  than  a  predetermined  value ,   and  retaining  those  that  have 
a  rise  time  greater  than  this  value.  The  calibration  of  the  rise  time 
circuitry  is  accomplished  oy   utilizing  a  signal  generator  with  pulses  of 
known  rise  times.   The  shape  of  the  input  and' output  pulses  are  shown  on 
Pig.  A-k.      The  differentiating  time  constants  are  equal  to  each  other  a_c 
each  of  the  three  settings  (l,  2,   and  5  usee),  and  can  be  varied  depending 
upon  the  desired  output.  The  counter  normally  operates  with  the  differen- 
tiating time  constants  set  at  2.0  usee.  For  an  analysis  of  low  energy 
spectra _,  the  differentiating  time  constants  are  set  at  1.0  usee,  and  at 
the  higher  5.0  [j.sec  time  constant  when  we  are  primarily  interested  in 
the  higher  end  of  the  energy  spectra. 

(2.e.)  The  output  of  the  Double  Differentiating  Circuit  is  directed  to  a 
unit  called  a  "Pulse  Height  Discriminator".   The  purpose  of  this  unit  is 
to  convert  the  input  pulse  to  a  square  wave,  and  from  this  square  wave  to 
generate  a  "start  pulse"  and  a  "stop  pulse"  as  the  input  double  differen- 
tiated pulse  crosses  the  zero  potential  axis.  The  electronic  arrangement 
for  this  component  is  shown  in  Fig.  A-5  along  with  the  circuitry  necessary 
to  produce  the  5  volt  bias  source. 

(2.f.)  The  "start  pulse"  and  the  "stop -pulse"  are  then  directed  to  a 
.Time -to -Height -Converter  unit.  A  block  discription  and  a  pulse  formation 
diagram  is  displayed  in  Fig.  A -6.  The  "start  pulse"  sets  the  Bistable 
Multivibrator  thus  initiating  the  ramp  generator  (consisting  of  a  constant 
■current  source  charging  a  capacitor).  The  charging  process  continues 
until  it  is  stopped  by  the  action  of  the  "stop  pulse".  The  "stop  pulse" 
is  fed  into  a  Monostable  Multivibrator  which  produces  a  2.5  u-sec  output 
pulse.  This  pulse  turns  off  the  constant  current  source,  and  thus 
terminates  the  charging  process.  The  height  of  the  ramp  at  the  point 
where  it  is  terminated  determines  the  height  of  the  2.5  u-sec  wide  output 
pulse.   The  capacitor  is  rapidly  discharged,  and  the  Start  Multivibrator 
is  reset  for  the  next  pulse.  The  output  of  the  Time -to -Height  Converter 
is  thus  a  pulse  of  magnitude  such  tat  it  is  directly  related  to  the  zero, 
cross  of  the  double  differentiating  circuit  output--thus  directly 
related  to  the  pulse  rise  time.   The  pulse  width  of  2.5  usee  is  of  no 
significance  other  than  a  convenient  width  for  the  following  components. 
A  circuit  diagram  of  the  Time -to -height  Converter  is  shown  in  Fig.  A-7* 
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(2.g.)  The  output  of  the  Time -to -Height  Converter  is  then  sent  to  a 
Single  Channel  Analyzer  (Fig.  A-8).   The  Single  Channel  Analyzer  can  be 
set  to  accept  ail  pulses  above  a  certain  height  (time)  or  to  accept 
pulses  on  a  window  basis.  The  analyzer  settings  are  adjustable.  Numbers 
on  Fig.  A-8  refer  to  standard  LKL  drawings  .   The  Single  Channel 
Analyzer  produces  an  output  pulse  only  if  the  input  amplitude  is  within 
the  window  setting. 

(2.h.)  The  output  pulses  of  the  Single  Channel  Analyzer  (rise  times 
within  desired  ranges)  are  then  passed  to  the  trigger  input  of  the  Peak 
Sampler  Unit  (Fig.  A-9).  The  original  pulse  has  been  building  up  to  a 
maximum  and  when  the  "trigger  signal"  arrives  from  the  Single  Channel 
Analyzer^  the  pulse  is  gated  into  the  Multi-Channel  Analyzer.  If  the  . 
"trigger  signal"  does  not  enter  the  Peak  Sampler  Unit;  the  pulse  is  led 
to  ground. 


*  Data  on  the  units  of  Fig.  A-8  can  be  found  in  Reference  35, 
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Consider  the  general  case  of  a  particle  of  mass  im,  and  v  colliding 

j-      1 

with  a  particle  of  mass  m  assumed  to  he  at  rest.  Particle  m  absorbs 

the  collision.  After  the  reaction,  particles  m  and  m,  with  velocities 

v„  and  v,  leave  the  reaction  site  at  angles  S  2-nd  Y   respectively. 
.?      4  L       L 

Utilizing  non-relativistic  dynamics,  the  laboratory  system  of  coordinates, 
and  mass  numbers  of  masses,  we  can  write  the  three  general  equations  for 
the  conservation  of  momentum  and  energy. 

a.  Conservation  of  Momentum: 

ml  vl  =  n3  Y3  cos  ^  +  mk  vk   cos  'V/L  (B-l) 

0  =  m,  v7   sin  SL  -  r\vj,  sin  ^'l  (3-2) 

"b.  For  the  Conservation  of  Energy: 

1  m  v  +  Q  =  1  m7  v  +  1  mk  v.  (3-3) 

2  2  -*  ^   2   '   r 

where  Q,  is  the  energy  equivalent  corresponding  to  the  mass  difference  (in 
non-relativistic  particles)  between  the  initial  masses  (before  collision) 
and  the  resulted  masses  (after  collision),  i.e. 

Q  =  (mx  +  m2)  -  (m^  +  rn,, )  (B->0 

Note  that  Q,  may  be  positive  or  negative.   Q  values  are  tabulated 
and  are  readily  available,  or  may  be  calculated  by  using  Eq..  (B-^) 

Equations  (B-l),  (3-2),  and  (3-3)  may  also  he  written  in  terms  of 

1   2       jr 

energy  by  making  use  of  the  relationship  E  =  ^  mv  of  mv  -Kj2   Em 
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Thus  equations   (3-1),    (3-2),  and  (B-3)    (ai^cr  clearing  of  the 
constant     $2     factor  in  Eos.    (f-l)  and  (3-2))  becomes: 

Jr\  Ei    =  K  S,   cos^,     +JraL   3'>      cos   S  (3-1. a.) 


r~r  _.    vv 


0  =     Jrn^  E^   sinb,     -^  E^  sin    ^  (B-2.a.) 

E1+   Q  =  S3-r  E4  (3-3. a.) 


In  the  following  analysis,  we  shall  work  with  energies  rather  than 
velocities.  ATote  that  by  the  simultaneous  solutions  of  the  above  three 
equations,  we  may  eliminate  any  two  parameters. 

Case  I  Simple  Scattering 


Let  us  first  consider  the  case  of  simple  scattering  wherein  m  =  m_ 
and  m  =  m.  .  It  is  readily  seen  from  Eg..  .(B-4)  that  in  this  case  Q  ='  0. 
It  is  desirable  to  derive  an  expression  for  the  energy  imparted  to  the 
stationary  target  mp  in  terms  of  the  incident  particle  energy.  Making 
the  above  mass  substitutions  into  Sqs..  .(.3-1. a. )  and  (B-2.a.)  transposing 
the  last  term  of  Ecs.  (B-i.a.)  and  either  term  in  (3-2-a. ),  squaring 
both  equations  and  adding  we  can  eliminate  the  angle  ou  to  obtain 


mi  El  "  2  Jral  ^2  El  \   C°S  X    +  ra2  \   =  ml  E5 


Solving  Eq.  (3-3. a.)  for  E,  and  substituting  this  value  into  the 


above  relationship,  one  obtains : 


or 


E„    -  2  /m.    m.  En    2,.      cos  V.  ■        +  ™~  s>,   =  rni    ^    "  EL  ) 
i.  V     1     2     1     4  t  2      4  l'u.  4 


■2  /m,    nu  E,    S,"      cos  t. 


y^  ,u2  ^  *4     —    .L  -r  -2  -^  -        ^  ^ 


\   0\  +  m2)   =  2./mLm2  3]_  E^ 


cos   V, 
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squaring  both  sides  and  solving  for  E.  yields 


4  '"1  m2 


C!  +  ^f 


E,  co 


.2  w 


(3-5) 


Therefore ,  the  energy  of  the  recoil  particle  is  dependent  upon  (a) 
the  masses  the  incident  particle  and  target  particles,  (b)  the  energy  of 
the  incident  particle,  and  (c)  the  scattering  angle  the  target  particle 
makes  to  the  direction  of  the  incident  particle.  ?or  a  maximum  recoil 
it  is  obvious  that  a  head-on  collision  (V,  =  6) "Will  impart  the  maximum 
energy  to  the  recoil  particle. 

Our  counter  is  designed  to  analyze  neutrons  of  energies  approaching 
I^-.O  MeV.   It  is  of  interest  to  determine  the  maximum  energy  that  can  be 
imparted  to  a  potential  recoil  particle.   In  this  case,  let  us  assume 
that  Vj_  =  0.  Using  the  atomic  masses  of  the  neutron  and  the  Eelium-3 
nucleus,  and  the  chemical  atomic  weights  for  krypton,  carbon,  and  oxygen 
respectively,  we  can  substitute  these  values  into  Eg,.  (3-5)  to  obtain:' 


Maximum  Recoil  Energy 

General         for  Ih-.O  MeV  neutron 


3 
1.  E     He  recoil 
max 


2.  E    carbon  recoil 
max 


3«  E    oxygen  recoil 

max 


k.     E    krypton  recoil 
max 


0.75 


j:k 


0.222  EL 


0.0466 


10.5  MeV 
3-95  MeV 
3.11  MeV 
0.65  MeV 


Case  II  -  He  (n,p)T  reaction 


In  the  case  of  the  He  (n,p)T  reaction  we  can  write  the  equations 
for  conservation  of  momentum  and  energy  as  follows :  Let  us  identify  in, 
and  E,  as  a  neutron  having  mass  of  unity  and  energy  ETT;  mass  m0  will 

be  the  He  molecule  having  a  mass  of  3j  w,2  shall  designate  m_  as  the 
triton  having  a  mass  of  3  and  energy  E^;  the  proton  having  a  mass  of 
unity  will  replace  m,  and  will  have  an  energy  S  .  Making  these  substitu- 
tions into  Eos.  (B-i.a.),  (3-2. a.),  and  (3-3. a.)  yields: 


\ 


=  p  Em  cos  6.  +  7  E_,  cos 
4    1      i-   v  ? 


k/5 


(3-1. b.) 


Ill 


0  =  J 3  £„  sinS.     '   +/T  sir 


+,/  Ep  sinY  (2-2-b.) 


*         *>P  T    '    -?  (3->b.) 

It  will  be  of  interest  in  the  determination  of  rise  times  for  us  to 
have  an  expression  for  the  angle  a, between  the  proton  vector  3^  and  the 
tariton  vector  2„.      This   expression   can  be  obtained  by  squaring  3cs.    (3-1. b.) 
and     (3-2.b.)  and  adding  the  two  equations.      This  yields: 


EAT  =  5  ET  +  2J  3  ST  s?     (cos    ^u      C0G^    "  sinfu  sin£   )  +  -E 

Recognizing  that  the  bracketed  tern  is  equal  to:  cos    (  V'     +    5T    )  = 

cos  aL       _,  we  can  solve  for  cos  a. 

e,t     "  ^  S_  -  E. 


■  cosa.   = 


2    (3  Em  3J1/2 


iji  *"p. 


The  solution  of  equation   (3-3 .b.)  for  S     and  the  substitution  into 


the  expression  for  cos  a.      yields; 

■  Q  *"  -;- 
.  n,p 

2/51 


(°'n  r)  ":"  2ETl 
cos  a,  =  -  l  ";p     xi—  (B-6) 


Note  that  a. is  always  obtuse  since  Q    is  a  positive  quantity. 
u  n_,p 

Q    =0.7o4  MeV 
n,p  • 

The  sin  a,  will  also  be  of  interest  in  our  computations.  This  can'  be 
easily  derived  using  the  Pythagorean  Theorem  to  produce 

fiz     (2  Ep  -  E    )    -Q2 
sin  a     =  '      — =*?  (2-1) 

2/5ETEP 


Case  III  He(n,d)D  Reaction 

In  a  dueteron  producing  reaction,  let  us  identify  the  general  case 
as  follows :   E  and  unity  are  again  substituted  for  3.  and  m  ;  m  is  the 
He  molecule  of  mass  3;  E_n  and  mass  2  are  assigned  as  E,  and  mv ;  and  3 
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and  mass  2  are  utilized  for  E,  and  m,  .  Substitution  of  those  values 
into  Eqs.  (3-1. a.),  (3-2. a);  and  (B-3«a.)  yields: 


u; 


^   %?2  ED1    C0S^  -  .  2  2D2  °03  V-  --I-  =  -> 


0  =    (2  Z 
4        Dl 

sin  £T      -;- 

.  "  ^D2 

E       +     Q         ,     = 

N         n ,  a 

^Dl     +     "J2 

sin  F,  (3-2. c.) 


(3-3. c.) 


We  are  again  interested  in  the  angle  c:;   cev.   .  bhe  two  deuteron 
vectors.  Proceeding  as  in  the  previous  case  we  find  that 

N       Dl      "D2 

cos  aL  =  -^ 

N/^Dl..  3D2 

Solving  (B-J.c )  for  2  and  substitution  into  the  value  for  cos  a, 

cos  a  =   — ;      — -  (3-o) 


kJ^ 


D2 
Note  that  if  E^.  +  E_0  =  Q,  .  (which  is  negative),,  then  a    =  x.     If 

ED1  +  ED2  "^  Si  d  '   aL     iS  acute>  and  if  ED1  *  SD2^  %   d'   ^L  iS  obtuse 
Q  ,  =  •■-  3.27  MeV  and  is  threshold  at  4.36  MeV 

At  threshold  a,  =  0 

The  sin  of  a.   can  be  determined  as  in  case  II 


i.     — (3-9; 

VE3I  S32 
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distribution  of  energy  and  angular  dependence  in  spherical  coordinates 

A.   Reason  for  the  Variation  of  Cos  0  rather  than  Q  in  Determination 
of  Rise -Time  Probability  Function. 


:iement  01  sunace  area  is  determined 


is  portrayed  in  the  attached  sketch. 


Assume  a  radius  of  unity 


■>  x 


dA  =  sin  0  d  0  d  0  =  -d  0  d  (cos  0) 


( 


c-i; 


If  we  wish  to  cover  every  possible  element  of  the  area  of  the  sphere, 
we  can  see  that  the  differential  area  is  not  composed  of  merely  d©  and 
d0,  but  d0  and  sin  0d0.  Varying  just  0  and  0  would  place  undue  emphasis 
on  the  polar  regions  where  the  sine  of  0  varies  slowly  with  0.  Therefore, 
we  must  vary  0  and  cos  0  to  obtain  uniform  distribution  of  the  tips  of 
our  energy  vectors  over  all  possible  configurations. 

B.  Probability  of  a  Lab  Energy  Particle  having  an  Energy  EL  after 
Interaction. 

Given:  A  particle  of  mass  rm,  with  energy  E1T  strikes  a  particle 


1 


1L 


of  mass  m   at  rest.   Product  particles  are  xn_  and  m,      with  added  energy  Q. 

To  Find:  What  is  the  energy  of  ny   in  the  Lab  system  of  coordinates. 

This  problem  is  similar  to  the  relationships  developed  in  Appendix  "3':, 
however  it  will  be  shown  that  any  distribution  of  energy  of  the  product 
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particles  is  equally  probable  within  a  maximum  and  minimum  energy  rang 

Notation:        m   =  mass 

E    =  energy    =    —    mv 

P    =  momentum    =    mv 

v      =  velocity 

0     =  azimuthal  angle 

0     =  ooiar  ancrie 

subscript  x,  y,  z  =  components  of  vector  quantities  in  x,  y,  or  z  directio 


subscript  1,2,3,4=     particle  identification 


subscript  L,  cm   =     coordinate   system; 


L=  Laboratory  system 
cm=  center  of  mass   system 


Select  the  coordinates  of  the  cm  system  so  that  the    z    axis  is  along  the 
direction  of    jv,  T 

i  JL  JLi     J 


SL, 


1L      -     2    '"1    [V1L 


m,       v1T 


or 


v 


1LI 


2E 


1L 


If     ml 


P    '    •  =    m,   j  v      I  =      2  m,     E 

M-M  1        1L  1        1L 


m 


1 
O— > 


v 


1L 


m2 
O 

|v9^  !  =  0 
!    2L  i 


Laboratory  system 
before  collision. 


m. 


o- 


m- 


v0  cm 


lvlcmi       |v2  CXI1i 


Center  of  mass   system 
before  collision. 


cm 


=      P 


!      T     ;P0  ;    T    |  ?- 

1  cm!        i     c  cm :       j     $  cmi 


V 


cm 


^1 
mj_   +  m.2 


v 


1LI 
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V,  j  =       'V  i  T     j     -         V 

1  cm|        |     1L  S         I    cm 


but  jv,  T    i!    and    I  v  are  parallel 

t    lLi  !  1     cm  1  r 


->-, 


m- 


V,  =      Vt  t  -      V  =      V,  T         1    -  r- 

lcmj  J   1L  cm  |      j    1L  I        m-i  +m- 


m- 


=   i  v 


lLi      m.  t  m? 
L    l  A 


-±»  . 


!vOT         -      v 


-£>    ! 


m-i 


=  -  v. 


and  |  v_         I    = 

I    2  cm!  I    2L  cm  |  .    lEt|m  +m 


.     Efr0t     before  cm)  ~     Total  energy  available  in  the  cm  system 

before  collision 


^  2 


'(tot.  before  cm) 


~     m.       v,   „      I       +    -r    m0      v->  „^ 
2  11  crni  2        2         2  cm 


1_    rill  m2       j  _v>,2 
2   m^_  +  m2     |  V1L.I 


=    E     * 


m2 


1L   V  mi    +  m^y 


and  E 


(tot.    after  cm) 


E  .    .     v    ,  +    Q  =E 

tot.   before  cm 


m. 


E    =  E 


IE      rn,  +  m~ 


+    Q 


(C-2) 


but,        E      = 


3  cm? 
v  2  m-j 


+ 


| -^4  cmj 

2  mA 


and  since  in  the  cm  system,    JP3  cmj 


j-^4  cm] 


>.2     r 


ip3  c^r  r  1 


1 

E       =    — *    j    —    +    

2  L  m3         m4j 


^> 


solving  for      j  P3  cm| 


2  mo  mA  E 


3  cm!     ,|      m     +  m4 
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A 


/ 


/ 


* " ^K  %, 


~ 


r 


->-x 


3x  cm 


S?_  "'     '  cos    (90-G)        COS  0 
!  3  cm .  v  cm  '  en 


'3   cm 


.n  9       ccc 
cm  •  e 


'3y  era     = 


!?.,        I  :os    (90-9)        sin  0 
3  cm  s  cm  rcm 


cm 


sin  9       sin   $ 


cm 


cm 


'3z  cm     = 


'3   cm 


cos   9 


cm 


OUo 


■3  xL 


;>x  cm 


3.  yL 


3y  cm 


'3   zL       =     m. 


p  v 

3z  cm       +        z  cm 

r.u 


"but  the  last  term, in 
brackets  =     v 


c  •  m  •  1 


>L  m 


',  ~  m.    1  v,  _ 

3  cm  I     cos  9Cm  +  1  I    11. 


^\+   V./ 


cos   9         2  inv     m, 
cm    i        .>> 


m_ 


'\J    (21- 


+■ 


2sl    E,T 

3  <\        1       1L, 


m,, )        (mL       rn2  ) 


IpH 

!  31! 


=    "p 


'3xL 


"3yi 


^Zij 


--*  id 


-    IT.,    '  f     (sin 
3  c4     } 


2   x  -29  .2 

)         cos     0         +  sin         cm     sin 

cm  r   cm  '  cm 


a2  0™i  + 


2 

cos      9 


cm 


2m-,  m,    E 


'3  "^ 


^\ 


™3  +  \J       (mi 


Pm     in        ^ 
L      3        1'-" 


4        cos   9 

r~jj  cm 


+  m    v~ 


1„  111—       111  1  j_i      _ 

J-  3      ^  ;.   IX 


+   m.  N 


n2)  Al  ft  -  *%) 
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■  5L 


o 


m,  m,  E '  +  ,  m,  m 
3   4 


^  1  3  "1L 


3 


cm 


t*\  +  ^n2;    /\|^  +  *\ 


m_  nv  m,  E._ 
L  3  4  1L 


I  3L 


m,  E 


J3L 


2  ^  '   "  (fflj  +  ^_) 


n  n  " 

1  3  1L  2  cos  0 

•7 r2  +       c: 

[ml  -  m  r         7~. 7-rr- 


m-m,  m,  3n_ 
1  3  4  1L 


j  /\j  (m  +  n,J 


Substituting  for  E:  from  equation  C-2  end  regrouping  yields: 


C   =  E 
3L    I 


LI   Uli 


ml  ^3 


\     -m,   Q     2  co 

2  +  ^~     +-r- 


i  +  —  2 1  +  -f ' 

^(m1  +  Dig)^  +  m^)      (m1  +  m^)  j  (m^  -;-  m^)   (r.^  +  n^) 


N 


,m?  Eit 
m.  m7  m,,  E-T(  — ^ — ^  +Q) 
i  3  4  H  hi.  +  m    ' 


(m3  +  m^) 


(c-3) 


Note  that  E„T  is  composed  of  two  parts : 
3L 

(a)  A  constant  value  dependent  only  on  the  energy  of  the  incoming 
particle. 

(b)  A  variable  that  is  dependent  on  the  cosine  of  the  scattering  angle 
0  in  the  center  of  mass  system. 

In  part  A  of  this  appendix,  we  discussed  the  fact  that  in  spherical 

coordinates,  the  parameter  to  vary  to  insure  complete  coverage  of  a 

spherical  surface  was  «he  cosine,  of  the  polar  angle  and  not  the  angle.  T/.7e 

can  therefore  conclude  that  since  the  cosine  of  0    can  vary,  any  value 

cm 

is  equally  probable  between  -1  and  +1.   Thus  the  Lab  energy  of  the  product 
particle  is  equally  likely  to  have  any  energy  between  a  maximum  determined 

by  cos  0   =1  and  a  minimum  of  cos  0   =  -I. 

cm  cm 


m. 


If  ir^  =  1,  m2  =  3,  ~* 

equation  C-3  becomes  for  the  proton  particle 


771    --"^  o  c*  nri  -f-V^i  TT/a »   [  yi  -o  r— i  _-,  — ^  --•- ^  p.^ 
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EPL" 


I  Ew   +  I  Q 


COS     Q  ,y 

n     kl  — r A   4       KL  ICL     rvp  (C-4) 


If  si  =  1,  m    =  3,  m     =  3,  m,    =  i     for  the  triton 


JTL 


o  EML  + 


cos   0 


% 


j   :    *  3IX  %.,y 


K1^     >  J 


The  total  energy  in  the  Lao  system  of  coordinates  EAr.  +  Q. 


:"!_   "n ,  p    pi 


3 


"EL  '  "^i,p   j  £■  K_  +  f  Q    +  cos  9 

;1    8  EL     k   ?i,p       cm  ik 

i fi A* 


Wl 


3  K„  Q 


3 

3 


EL 


Q    -  cos  © 


k   n>? 


V;9  2 


Since  the  scattering  angles  in  the  cm.  system  are  equal  and  oppc;;ite 


J3  -r     0  =    2 

EL       %p         EL 


."o 


If  m.,    =  1,  3?     =  3  j  ku  =  iQi    =2     as  in  the 

-L  ^  ^  LT 

reaction,,  equation  C-3  hecar.es  for  deuteron  I 


case  for  the  He  (n,d)"D 


-  if, 


r. 


!y>  2 

TZ  =     =■   (E._  +  Q  .)  +  cos  0      .'-  EA_  +  E._>  Q 

D  IjL.    2  V  NL    n^°y  cir.  /\  ■    NL     EL  h,d 


2 


E  \'k   - 


(C 


Deuteron  '  2  is  equal  to  deuteron  '  1,  except  that  when  addin; 
to  determine  total  energy,  the  sign  of  one  of  the  scattering  angles 
must  he  reversed. 
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Appendix  D.  Derivation  of  Field  Strength  and  Voltage  Souipotentiai 
Equations  for  Coaxial  Cylindrical  Geometry 


E(r)  -  voltage  gradient  at  any 
radial  position  r. 


V(r) 


V 


=  voltage  at  any  radial 
position  r 

=  radius  of  the  cathode 

=  radius  of  the  anode 

>pli« 

anode 

-  applied  potential  to  the 
cathode 


lie     r?    V    r. 


;hus    K?~Y  =   0 

.Boundary  conditions 


(1)  V(r  )  =-V 


(2)  V  (r>  0 


in  cylindrical  geometry 


r^V 


=  C. 


^V 


\ 


V  =  C1  inr  +  C2 


Applying  the  3.  C 


V(ra)  =  CL  Xn  ra  +  C2  =  Vq 


V(rc)  -  Cx  J?n  rc  -,  C2  ^ 


(jB .  C .  no .  1  } 
(^B.  C  no.  2  J 


Solving  the -above  two  equations  to  determine  C  and  C~  we  ohtaii 

°  "  i      2 
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V(r)  = 


.in: 


V, 


or 


V(r)  = 


A 

"a 


but         E(r)  =   -grad  V(r 


v.N 


E(r)  =  Vc 


ria 


a 
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Appendix  E 
Random  Probability  Analysis 
CATMAN 
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c 
__r„ 

DM  IS  THE  CEUTERCN  MASS.    EN'H  IS  NEUTRON  MASS.  FE3M  IS  HELIUM- 
■>.    y/i^.  fi  vni  t  rs  tmf  fi  Fr.TRnN-vni  t  tm  frgs.  fni  im  is  the  max. 

c 

ENERGY  OF  THE  CHANNELS.  GNC  IS  THE  CEUTERON  REACTION  ENERGY. 
ON?  r<^  THF  PftnTr.N-TRITnN  RFAGTIUN  FNFRGY.  PIFMTN  IS  THE  MINIMUM 

c 
c 

DIFFERENCE  BETWEEN  THE  SQUARE  OF  THE  MAX.  RADIUS  AND  THE  SQUARE 
nf  THF  MTN.  RADIUS  sn  THAT  THF  risf  tiff  is  LONG  ENOUGH.  SIGNP 

c 

IS  THE  PRQTCN-TRITCN  REACTION  CROSS  SECTION.  SIGND  IS  THE 
nFHTFRnN  REACT  \r-\     CROSS  SECTION-  RPA'R  IS  THF  PROTON  RANGE  IN 

c 

AIR  DATA.   THE  CROSS  SECTIONS  ARE  IN  8ARNS  (10.»«f-24)  CM.»*2). 
FNTARS  =  Fl PAT (  N  T  A  3.  S  ] 

A!\RAD  =  .015*2.54 
r.ATRiC  =  l-965*?-54 

CAT2  =  CATRAD»*2 

r. atm am  =  r. atp. ah  -  anrah 

CYLENG  =  8.5*2.54 
Y«AX  =  7-S*?.54 

PRESS  =  12.005*76. 

PNM  =  1.67248E-24 

T  V  M  =  5  .  n  ?  ?  n  K  F  -  ?  4 

DN  =  3.34728E-24 
F,\M  =  1  .A74RF-74 

HE3f  =  5.01976E-24 
Fl  Vni  T  =  1  „A0?07F-?? 

E  N  L I W  =  15.E6*ELV0LT 

FRANGF  =  FNI  TV  -  1  „  F<S*FL  VOLT 

C.NO  =  -3.27E6*ELVOLT 
C\'P  =  -7^4Fft*FI  vni  T 

TFACTR  =  PRESS»ALCGiCATRA0/ANRAD}/(2.*.721£6*V0LTS) 
TFMP  =  ?q^.l-ft 

DIFMIN  =  TIME/TFACTR 

RFAIIf?.?)   (^T^MOff  ).!  =  ltNTARS) 

READ(2,2)  (SIGNC(I)„I  =  i,NTAES> 
RFADf?.?)  (RPATRtl),!  =  l.NTAHS) 

2  FCRMAT(  8F10.0) 
41(11  =  0. 

A2(l)  =  0. 
a  a  ( 1  )    -    n  „ 

A7(l)  =  0. 

idlll  =  ?-»STfiNPI >  )     -    SIGNPI2) 

SPACNG  =  ENLIM/FNTABS 
SRFf.  IP  =  FNTARS/ENLIM 

RCHNKD  =  FLOAT  (NCHAi\S  J/ENLIM 
CAM  CLOCK  I  (TTM1  ) 

CALL  LENGTH(NTABSiENLIM,  ENER,  ERG  t  RP  A  IR  ,  EL  ,  ELM  AX  ,ELSP  » ELL»  EEL ) 
CALL  CLOCK  I  (TIM2) 

C 

CLOCKI  IS  A  LIBRARY  ROUTINE  TO  PROVICE  TIMING  INFORMATION. 
DC     160  NEUEN  =  l.NENS 

EN  =  ERANGE*IELCW  +  FLOAT ( NEUEN ) * ( EH  IGH  -  ELOW) /FLOAT INENS ) ) 
TNTFRP  =  INT(FNTABS*EN/ENLIM) 

PART  =  (EN  -  ERGMNTERP)  )*SRECIP 

STOP  =  SIGNPt INTERP)  +  PART* ( S I GNP { INTERP+ i )  -  S  IGN? ! INTER? )) 

r. 

SIGD  =  SIGNC{ INTERP)  +  PART* ( S IGND ( INTERP+1 )  -  S  I  GNC ! INTERP ) ) 
FN  IS  THF  NFUTRCN  ENERGY.   SIGP  IS  THE  P-T  RE  AC  T  I  ON ..  CRO  S  S  „„ 

c 
r 

SECTION.   PNEUTL.  IS  THE  NEUTRON  MOMENTUM  IN  Mrr*,~    LAS.  AND  PNUTCM 
fS  TN  THE  r„y„     SYSTFM.  VC.MSYS  IS  THE  VELOCITY  OF  THE  CM. 

c 

SYSTEM,   ETCM  IS  THE  TOTAL  ENERGY  IN  CM. 
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SIGTGT 
-CHOICE. 


1 


.WTFCTR 

PNEUTL 
VCKSVS 

PNUTCM 


=  SIGP 
-sl-S-IGP/ 

=  .6025 
.(22400- 

=  SQRT( 
.=  .PNEUT 
=  PNEUT 

-PNUTCM  ■:■ 


=  VCMSY 

PHE3CM-* 
ENCM  + 
UEN>-=~- 
EUEN)  = 
1:2X1    = 


S*273.16*SIGT0T*( 2  .*CYtENG )*CATPAN*3. 141592  65/ 


»ENJ 

...+..  HE3M)  . 

MSYS*ENM 

.;  c  \  y  ) 


&HE3M-L-- 


-LA- 


TCTAL(NEUEN)  =  0. 
-DC— 16-  I  =  1.+-NCHANS 

CHAM (I )  =  0. 
_C-CNT-L\U-F 


NEU 

\  F  l_! 


TRN  =  0 
IRK    =    NEU 


1£AL 


STATEMENT  17  BEGINS  THE  RANDOM  GENERATION  LOOP  (SEE  STATEMENT 
-LAO  BELOW-)..—  WS-_CAN_.CCNS.IDEA  .ALL  EVENTS  TO  OCCUR  TO  ThLE_RUGlHI. 

OF  THE  CENTER  OF  THE  CYLINDER  (IMAGINING  THE  CYLINDER  TO  LIE 
J-dnaLZ0NXAU.-Ya.-3  EC A  USE-  0  F_  ..S-Y.---EX-.Y.  I  h  IS  IS  SLMPLFR  BECA  USE 
NO  TRACK  CAN  REACH  AS-FAR  AS  THE  LEFT  END  SINCE  AT  15  MEV  OR 

i  F<;<;:   .■"  '      -?*r:'<^     A«-  SHP3TFR  ~H-'.M  ":-'A"_    Tj-F  END  HAii   F c  F  F  P  ~  ' 's 


THEN  DUE  CNLY  TO  THE  RIGHT  END.   TO  CHOOSE  THE  REACTION  •  S I TE 
JLAND.QM— :  N  THE— C-VLIAi:  E-i-U  PIIIMF  '.-.' F  P.HPPSF  HP.  THF  wr)RT7P,NTA1 
POSITION  MEASURED  TO    THE  RIGHT  FROM  THE  CENTER,  AND  R,  THE 
RADIAL  DISTANCE  F  aOJ^JLH£-£EAL7-ER— AXJ-S^. THF  REACT  IDN._S.:  IF.?. AJ_x 


AT 


NOT  UNIFORM  IN  THE  VOLUME  SINCE  THERE  IS  MORE  VOLUME  IN  R  TO 
o^no  pro  i  aqpf  ?  'u>v  fhr  svam  r.  s^np.f  '■"=  arf  p.hppstnp  thf 


VALUE  OF 
-CLLSL3TRUT 


R  UNIFORMLY,  WE 
1£__. SINCE  THF 


MUST  APPLY 

VOLUME  IN. 

WEIGHT-  MUST  3E 


HEIGHT 


CORRECT 


PROPORTIONAL  TO  R?  THE 
i'HEM  AN  EVENT  IS  CPU  H I  ELa 

THE  WEIGHT.   IN  ORDER 


TO 
>  y  A ' 

PROPORTIONAL  ~0  R. 
I  TO  THF  PRO?- 


BUT 

_S.?.£i 


»'E  AD 
■oj  c 


S'v-RAI 


TO  COMPARE  THE  RESPONSE 

>:^"-RP,y  -'-^P'-S  T7  v.r.'.'in  HA"-  ?.-FN' 


NECCESARY 
XiLQS-S—SEC 


TO    USE    A    Mi 

LIEN— M.A.8  I  F  S 


IGHT 

y. :  T  H 


ANYWAY 
NEUTRON 


TOTAL 


WTFCTR     IS 

■XAADIUS    0 


SECTION). 


EACTION 
E_J.:_iGB_I_-. 

(NUMBER    OF    HELIUM    ATOMS/CC )«( LENGTH    OF    CYLINDER)«P 
f   rvi  TNr<=s    trnRRPr--n   fpr    sp>  m   ANnnniciTfiTM    ( 

IMAGINE  A  COORDINATE  SYSTEM  AT  THE  REACTION  SITE 

t  <;     oao>.!  I  c;     -n    tk-     f^nCF     AMP     7-AX'S    A'PliXP     a,    RADTAL 


:  ■-■ 


LINE  FROM 

:ht<;  s  y s t 


cosine  c: 

P7C.K  THF 


THE  ANODE.   WE  CHOOSE  A  NEUTRON  DIRECTION  RELATIVE  TO 

FN*    ry    pwppstnp    A7  tmdthai A_L_LF,    PH!?    AT    RANJIOM    AND 

THE  POLAR  ANGLE,  CTH,  AT  RANDOM.   IT  IS  NECCESARY,  TO 


:OSINF  AT  RANI 


DISTRI3UT 


ION. 


the  c:m. 

NEUTRON 


SYSTEM  IS  IMAGINEO 
niRrCTIOXU 


TO 


TIONAL 

HAVE  X-AXIS 


HP  = 

R  = 


CYLENG* 

P.ATMAN-R 


RGEN(UGH) 
GFNiUGH)  + 


AN'RAD, 


R2 


R*R 

_WTFP.TR-_i 


TC" 

O'rl.' 


AL(NEUEN) 

=     6  -  ?  P,  ?,  - 


=  TOTAL (NEUEN) 


+  W" 


CCSPHI  =  COS (PHI) 
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SINPHI  =  SIMPHI) 


STH  =  SCRT.I.  -  CTHas2) 

T  F  (  R  G  F  ^tULGKL_ j^L.. CHOI  C£-__GJD__TI!_1S 


CHOICE  IS  THE  FRACTION  OF  EVENTS  THAI  ARE  PROTON-TRI  TCi\  EVENTS. 
_LF__A._NUKB£2...RAM3CWLY...CiiQS.£JM  .  BETWEEN..  C_.AJ_t_l._I  S_.GI ...  C___,_  WE 
PICK  THE  DEUTERCN  CASE.   IN  THIS  CASE  WE  MERELY  SUBSTITUTE  THE 
^■-=P.py  v**,*.   :-na  thf  prpitp,..  and  ts.tpn  masses  and  qno  for  the 
REACTION  ENERGY. 


P,v 


_C__ 


KASE  =  C 

p.  p  t ^  iq 


PN  =  DM 

TV  =  HM 


0  =  GND 
_ASF  =—2- 


19  PPCM  =  SQRT .2.*(ETCf  +  Q.*PM»TM/(PM  +  TM  >  ) 

PHTP  =  /i.7H^18S?l  --RPFNflir.H) 


CTHP  =  2.*RGEN (UGH)  -  1. 

STHP  =  Sj__TJ_.___.CTH.P._'2  )  ._ 

XPCM  =  PPCM-STHP-COS(PHIP) 

YPCM  =  RE._lJ_J.__.--.  h.  ( ______ 


z?cy  =  PPCM -CTHP 

XMA  =  X^P.:V*.STH  -  ZPf.M  _~  H. 

ZMB  =  ZPCM*STH  -:-  XPCM-  C7H 

xvr  =   xvA»r.nspm  -  y?c.v«si ,\  p.h  l 


YNB  =  XMA*SINPHI  +  YPCM-CCSPHI 

PPP.M  IS  THF  PROTON  Mr:-;  £„  I.U  _____  C  .„ X  M  g  -  yj_  .  Z K  B  j_A_R.E_IHE_X_jY.ti_. 


C  COMPONENTS  OF  PRCTCN  MOMENTUM  IN  A  COORDINATE  SYSTEM  THAT  IS 
_£ °A  .All  Fl  TP  THC  !  AR  SYSTEM.,  RU~  MOVING  WITH  THF  C  _  f  _   NOTE  THAT 

C  THE  COMPONENTS  OF  TRITON  MOMENTUM  IN  THIS  SYSTEM  MUST  BE  THE 
X __EJ_______L__J_j_Y__I3______B„ . 


VCSYSX 

VP.SYSY 

=    VCMSYS     _TH*CCSPHI 
=    VCMSY5.     STH-SINPHI 

VCSYSZ 
PI  ARXP 

=    VC    _YS*C7H 

=    PM-VCSYSX    ->    XMfi 

PL.A8YP 
PI  AH7P 

=     PM-VCSYSY    +    YMB 
=    PM-VCSYS7     +    ZNB 

PLABXT 
PI ARYT 

=    TM-VCSYSX    -    XM8 
=    TM*VCSYSY    -    YMfl 

PLABZT 
P!  XP?     = 

=    T. '-VCSYSZ    -    ZMB 
:     p;  APX?  -_? 

PLZP2  =  PLABZP**2 

XP7P  =  PLXP?  +  PL I P 2 


PLABP2  =  X?ZP  +  PLAEYP- -2 

PI  ABP  =  SORT!  Pi.ABP.2J 

PLXT2  =  PLABXT««2 

?'■  77?    =     ?■   ARZT*^? 


XPZT  =  PLXT2  +  PLZT2 

PI  ART?  =  XP7T  +  PI  A ____*_ 

PLABT  =  SCRT.PLABT2) 

FP  =  PLAB_2/„__s_i_ . 


T  =  PLABT2/(2.*TM) 

pp  Ayn  ft  .  f-.c  "'  _  P'  _70N  AN.  "...  .TON  ENERGIES  IN  'HE  LA!..   PLEV: 
AND  TLENG  ARE  7\-ic    TRACK  LENGTHS.   XP.YP.ZP  ARE  THE  X.Y.Z 
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COORDINATES  CP    Tl 

C.  ?  "  C- .  '• :  ''  ~  t  -  -  C  V 


E  TIP 


OF  THE 


PROTON  TRACK  IN  A 

,    v-£,--::s   ■u:_^:^•  r;  the. 


SYSTEM    WITH 

a  m  n  n  P      .■•-. '  D 


Z-AXIS     PARALLEL    TO    THE    RADIAL    LINE    FROW    ANODE     TO    REACTION    SI" 
-INTP    =     INT(-S&£C-IPgEg-J 

IISTT    =    INT(SRECIP*ET) 
■PJXR-T-P— J5— (-EP   -    EBGtl-N-T-P  U-itSRECIP. 

PARTI    =     iET    -    ERGUNTT))»SRECIP 

c^-mh   -    p  km-o   ",^c.vy,    ■■•rT?-->!  '•'ir?^;'(A';e»'!-fi!;?'T?-::'  vV 


__ 


TLENG    =    EL!INTT;:<; 
p^TIC'P    =     P|  FMH/P1 
RATICT    =    TLENG/PLA8T 
.XP— s    PLA&XP-gJLAT.IQP-  — 

YP    =    PLA3YP-RATIQP    -:- 

,zp  -  pl/-52?-.--:at::.:p  v 


TT*IEL( INTT+l,KASE+2)-EL( INTT.KASE+2) ) 


XPSC    = 
-2-2-SQ— s. 


R?2 


=    X 

—2.L 


XP*XP 

1^-12 

PSC    -:-    2PSG 
ASOXiR-AXXOZ 


PL 


ABY 


ToRATIOT  +  HP 

Xr-RATIQT  -:•  R 


XTSQ 
-Z-XSC- 


=  XT* 


XT 
ZT- 


RT2  = 

1  P  V  'J  P_ 

YTMHP  = 


TSC 
YP 


v  ZTSI 


-  HP 

'AY  -  UP 


IF(RP2 


.GT, 

XOJJ 


CAT2)  ALPHA  = 

,:j:~-io?) 


RMSGRTC CAT 2*X? Z P/R 2-P LXP 2) -PL A  8 Z?)  / 


(R 


72 
_LX 


.GT 

PZT 


.     CAT2) 
iRJLLLCU- 


A  =  R»(SCRT(CAT2*XPZT/R2-PLXT2)-PLABZT)/ 


C 

ALPHA  AND  BETA  ARE  THE  FRACTIONS  OF  PROTON  AND  TRITON  TRACK 

•  p^TH^  THAT  RFMATN  T  N  S  T  n  F  THF  f.v-  TNPFR  WA!  1  .   THF  PROGRAM  FROM 

C 

HERE  TO  STATEMENT  12C  DETERMINES  IF  THE  RISE  Tiy£  IS  3IG  ENOUGH 

Avn  tc  sn,  w-AT  fwfrgy  Tn  r.niiMT.  hfnfra:  i  y,  thf  trat.-.s  are 

C 
r 

SKEW  LINES  RELATIVE  TO  THE  ANODE  ANC  THE  MINIMUM  DISTAi\CE  MUST 
pF  FOMNP  TT  HFT  TH'F  RTSF  TTMF,  YM'MP.  TS  THF  Y  COORDINATF  flF  THE 

C 
r 

POINT  OF  MINIMUM  DISTANCE  IN  THE  LAB  SYSTEM  WITH  ORIGIN  AT  THE 

OCA<--r7r>N,   <;TTc> 

IFIHP  .GT.  YMAX)  GO  TO  60 
'r-ivo  .ip.  y y  a ::  i  an   rr    ? 4 

IF{RP2  -LE.  CAT2)  GO  TC  20 
yniT  =  Yp^hpsaj  pha 

IFtYOLT  .LE.  YWALL)  GO  TO  25 
?n  B5UAI  B  =  MRePI  1RVP  +  PIAR7P»YWAII  lee?  -;•  PI  XP  ?«YWA1  1  *»? )  /PI  Aft  YP»#? 

12  W  A 
IHER 


LP 


IS 

iiE_ 


THE  SQUARE 
PafTDN  TRA 


OF  T 
CK_GO 


HE  DISTANCE  FROM  ANODE  TO  THE  POINT 
PS  THRU  T!-F  I  TVF  SPACE-DEAD  SPACE  WAS. 


IF (PLABZP  -GE.  0.  )  GO  TC  22 
YV[N3  =  -Pi  AftYPsPI  AR7P«R/XP7P 

I  F  (  Y  M I  N  R 

0  V  T  \-  7  0     - 

.GT.  YWALL)  GO  TO  21 

PI  »OT;  3  ?  /  Y  o  7  D 

RMAX2?  = 
ar    rr  ?t 

AMAX1 (R2 »R2WALP) 

21  R,VIN2P  = 
c?v/\Y2P  = 

R2WALP 
R? 

GO  TC  23 

??   0  V  T  V  9  0      = 

o-> 

RivAX2P  = 

R2WALP 
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23  TIPP  =  (1.  -  YWALL/YPfhP) *PLE\C 


24  IF(RP2  .LF.  CA72)  GC  TO  31 
_25  -L£±2JLAP,7.?     .GF.  n  .  3  _£  CL_Tr  ?2_. 

RfI\2P  =  PLXP2*R2/XPZP 
GE_X0-_2-8 . 


27  RfIN2P  =  R2 
■>  o  o  v  a  y  '->  o  -  r  A  T  ? 


TIPP  =  (  1-  -  ALPHA) aPLENI 
TIPP  TS  THF  I  FN G.~B__C£_ 


!J±£_  Sl\C_Jdf_?.RtlT.CN_IRAC 


AI_JLS_CLlCL£iDJ 


OF  THE  CYLINDER. 
_jGjG 1 0  33 


31  IF(PLABZP  -GE.  0.)  GO  TC  33 

vy-"P     -     —  p|    '  ::yp--,pi    A  ?  7  P  -'R  /  X  ?  7  P 


IF(ARS( YMINR)     .GE. 
RyT  N7P     =     PI  X 


ABS;YPMH?) )    GO    TO    32 


RivAX2P    =    AMAXi  «R2,RP2) 

c,r.   ~  fi    ^4 


32    RI^IN2P    =    RP2 

o  v  A  y  ->  ?     -      p  7 


GC  7C  34 
3H  Ry7K?P  =  R7 


RXAX2P  =  R?2 
3A  TTP?  =  n. 


>9  IF(YT  .  LE.  YKAX;  GO  TC  44 

-"i'?7?  .i  -.  gat?)  r,r  tg  4n 


YCUT    =     YT;vHP*cETA 

IF(YCUT    .LE.    YisLAI.L)     GC    T.0_45_ 


40    R2WALT    =     ((R-PLA3YT    +    PL A8ZT*YWALL )«*2    ♦    PLXT2*YWALL»*2 ) /PLA3YT»*2 

IFtPLABZT     .GJi.     fl.)     GO    TC.    4? 

YCINR    =    -PLABYT*PLABZT-R/XPZT 
T  =  ;  v '/  T  \i  R     _  a  ~  -     Y  '■-■'  AM)     G  G    ~  r    41 

RPIN2T    =    PLXT2*R2/XPZT 

R.VAX2T     =     A.VAXMR?,  R7VMI  T) . 


GC  TC  43 
Al  RNTN7T  =  R7WAI T 

RyAX2T  =  R2 
r,r    tp    a,^ 

42  R,VI.\2T  =  R2 

RVAX7T  =  R7WA1  T 

43     TIPT    =     (I.    -    YWALL/YTiwHP)*TLEi\iG 

or.  TO  5  5__ 


44  IF<RT2  .LE.  CAT2)  GC  TC  51 

a=:  TFfQ'AR?T  .  Gf-'.  fl.  )  GO  Tf  4? 


RPIN2T  =  PLXT2«R2/XPZT 
nr  ~n  ar 


■47  R,VIN2T 
4  8  RVAX7T 

TIPT  =  (1.  -  6ETA)*TLENG 

r.  r    -  n  SJ5 


R2 

_CAT2_ 


51  IFIPLABZT  .GE.  0.)  GO  TC  53 
V^7\R  =  -Pi  AftYT*Pl.AR7T*R/"  ?ZI_  . 


IF (ABS (YMINR)     .GE.     ABS(YTMHP))     GO     70     52. 

»»fM?T    8     PI  XT7»R7/XP7T ! 

RNAX2T    =    AiVAX;.  (R2,RT2) 

r,r     -  r     S  A 


52    RPIi\2T 


RT2 
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RVAX27     =    R2 

ar.r.r.,.^: 


53  RMN2T    =    R2 

RXAX  2-T-~    JU2 .._ 

54  TIPT    =    0. 

&.5-RM  N2  -=  -  A-V  I Nl  1  R.V  I.\2P ,  R  .v  I N2T  ) 

R^AX2    =    ANAXKRf  AX2P,R^AX2T) 
IF(RKAX2.  -.    RMN2     .LL D IMMJ.U )._ GO-  Tu    1?.C 


IF  (TIPP    .EG-    0. )    GO    TC    56 
4L_=— INJ-U".IP-PV-ELSIUK.ASE*1  )  ) 


cTIPP  =  EELCIL,KASE+1 )    *     (TIP?  -  ELL  (  IL  ,  KASE  +  1 )  )  »  (  E£L'(  I  L+  1  ,K  ASE  +  I ) 

ECEPP  =  EP  -  ETIPP 
GC— 4-G— 5-7 


56  EDFPP  =  EP 

_57—  IF.t-T.-IPT-.EO-.  0.  J--GO  TC-53 .  .. 

IL  =  INT(TIPT/ELSP<KASE+2) ) 

E-T-I-?T.-=-SEL.CIL,X.ASE*2  J  *  .  {.TJLP.T.  .-.._£LLI  IL4KASE.±Z11*  (  FEL.Ui.-fcl  ,KA_SEfc2  i 

I  -  EEL(IL,KASE+2))/ELSPCKASE+2) 


GC  TG  59 

5£  EDEP-T— ~-ET ..  -  ..  

59  ECEPCS  =  EDEPP  +  EDEPT 
-C-     EOEP-OS  IS  TKE-TL0.T-A1 ENERGVL 


LE20  £II£C_.  !__"_£.  JL  LVE_  _PJ_£, 
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GC  TC  120 


v  k  a  y  i  r;  f  T  r 


T) 


IF!YT  .GE.  Y.VAX)  GC  TC  130 
-I  F  ( RX2— --LE CAXL2;.~-GC-  TG  65. 


61 

IF(YTNHP*BETA       -GE.    YWALL)    GO    TO    I3C 
TC'°?A3ZT    ,GP.    C.)     GO    T C    6 ? 

YNINR    =    -PLABYT«PLA8ZT«R/XPZT 

;F'Y''TVD    -^T,    vwa|  i  )    nn    yp    ffP 

RPIN2T    =    PLXT2*R2/XPZT 

GC    TC    63                                                                               .             .                . 

62  RPIN2T  =  UR«PLABYT  +■  PLABZT*YWAL-L  )  »*2  +  PL  XT2*YWALL**2 )  /PLA3YT*»2 

-63.  IF1CAT2  -  RMIN2T  .LT..  C I  F  M  ':  N  )  G  a._XL__£ 

TIPT  =  TLENG*(1.  -  BETA) 
l-'IPT  =  "'.  -v';-~  '  '..  -  ',"■', :.:  /YTMHP.i 

WTIPT  IS  THE  LENGTH  OF  TRITON  TRACK  THAT  EXTENDS  TO  THE  LEFT  OF 
THF  LIVE  SP_  CF-DFAD— S-P-ACF  HALL,.  _H£R_*_  T^F  RFACTTCN  SITE  IS  IN 

THE  DEAD  SPACE  BECAUSE  n?     IS  GREATER  THAN  YMAX. 
i!  =  TNTfT-pT/F!  S?(KASrS-?_ 


ETIPT  =  EEL(  IL»KASE+2  )  +  (TIPT  -  ELL  (  IL  ,  KA  SE  +  2  )  )*  <  EEL  (  I  L+ i  ,!<ASE+2  ) 

-  f  f  •  (  t  -■  .  ;<  a  s  f  -•-  ?  n  /  f  ?  s  ?  ■:  :<  /,  s  f  -.-  ?  1 

IL  =  INT  :;aTIPT/ELSP(KASE  +  2)  ) 

F '.-;~?T    =    FFI  f  T I   ,:_vSF->?)     -i-  (  VH1  £_..-    F!  i  (  I  L___S  £±2)  1*(c£LJL:  L  fcl  .KASE*?) 
-    EELUL.KASE+2) )/ELSP (KASE+2 ) 

-ECEPCS    =    E'.-JTPT    -    ETIPT 


GC    TC    120 

AS      Q7i.'flT      =      (   ■'  P  -PI    A^.v 


PI   A  n.  7  T  -r-  V '<■!  A.  '   !    )  '■>  '■>  7 


?!  X-?-YWA!  1  »»7)/PI  ABYT**? 


IF(PLABZT    .GE.    C.)    GO    TC    66 
v:;TM     =     -?l  ArlYT---P;  AP7T-JR/XP7T 


IFCYMINR    .LE.    YTMHP    .CR.    YMINR    .GE. 
RVM7T     =     ?!  XT7*R7/XP7T - 


Y'rtALL  )     GO    TO    66 


:  F  I  a  max:  (RT2 iR2WALT) 


RMIN2T    .LT.    Bi^IN)     GO    to    13C 


66 


3SCRT2    -    R2WALT!     .LT 


FMIN)     GO 


*    _         1  J  V> 
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67    WTIPT    =    TLENG«(1.    -    YWALL/YTMhP) 


EDEPOS    =    EEHIL.KASE  +  2) 

J -  £  r.  L  I  I  L  ,-KASi 


(WTIPT-ELLt IL,KASE+2) )  « (EEL ( IL+ 1 ,KASE+2 ) 

,2_LKASf.-;-2:i 


GC  TC  120 

70  IFtYT  .  Ill 


jm^ao_ 


IF(RP2  .LE.  CAT2)  GC  TC  75 

r  c  i  v  o  fj  u  p .-. '. '  p  h  a  .  r,  - .  v  w  /*.  i  :  )  r-  n  t  n  \  ~  n 


71     IF(PLABZP    .GE.    C.)    GO    7C    72 

YVTMR    =     -?:   ABYP--PI  AR7?*R/-,:P7? 


IF(YMINR     .GT.     Y'aALL)     GC    70    72 
R.VTN7?     =     Pi  XP7*3?/XP7P 


GC     7C    73 
■77     pvT\!?o     -     ffPr-oi   ^"P    -i     P1   A?7P---Y'-'A:  '    )o*?    <■     ?!  X°7*Y',-.'Al  !  **2 )  /PLABYP^-2 


73  ifccat2  -  rpi>\2?  .lt.   difmin)   go  to  13c 

:~°?  =   ?;,fmgj>jj ^.j^lphaj 


TPP  =  ?LE\G*;i.  -  YWALL/YPMHP) 

=  T  M3L!  T  T  P  p  /  F  i  S  P  LKASF-;-  i  )  J 


TPP  =  EEL( I L, KASE+1) 

-  F^'  (  '.•   ;"ASF 


•:■  (71??-  ell(  il, kase+1)  )*( eel ( il+ i » kase+1 

-  '  \/?\  s?;;<asfv;  ; 


IL    =     INT(WTIPP/ELSP(KASE+1) ) 

fvtpp    =    ffi   (Ti  .k^Fv  LJL_+ _  LU  TIPP    -  E  LL  LI  Lt&AS  E  +  1 )  )  »JLE  cJULLLt.  L*J<AS£+JL1 
1  -    EEL(  IL, KASE  +  1) ) /ELSP (KASE  +  1 ) 

^r--?"^   =   fwtpp  -   ft:?? 


GO    7C    120 
75     P:?MMO    =     'j^n;  jrv?    v.     r>  ■  a  g?  p  .-.yw  A:  :.  1  ■-•  ■-■  7     •>    PL  ".?  2*  YL'ALL  *  o  2  )  /  PL  A  3  YP**  2 

IF  {  PLABZ?    .GE.     0. )     GO    TC    76 
Y,VI.\R    =    -Pl.ARYP=JJLAa.Z£-g-R/XP.ZP. 

IFfYMINR    .LE.     YPMHP    ,CR.    YMINR    .GE.    YWALL)     GO    TO    76 
RV7MPP    =    pm:p?^r?/x?7p 


IFfANAXl (RP2,R2WALP) 
nr    I n  77 


RMIM2P  .LT.  OIFyiN)  GO  TO  13C 


76  IF(ABS(RP2  -  R2WAL?)  -LT.  CIFVIN)  GO  TO  130 


7  7    '.-;  ~  I  ?  ? 


.PLFNG«(1,    -   Y.ltLAL L/XPJ1HP.L. 


iL  =  INT(WTIPP/ELSP(KASE+1) ) 

:nFPn<;  =  FFl  f  T!  .XASF+1  )  +     '.  '.-.'.I.  IP?- ELL  {  I  L  .  KASE+  HI*  (  EEL  (  I  L+  1  .  KASE  +  1  ) 


1  -  EEL(IL,KASE+i))/ELSP(KASE+l) 


80    IFIRP2    -LE.    CA72)     GC    TC    ICO 

•r'RT?    .1  F.  .C.AI2  LJ3C    TC   9J1 


IF(YPPHP*ALPHA     .GE.    YWALL)     GO    TO    61 

rF^YTN'HP^RFTA .GF.     Y'/iALL  )     GO    t q_ij 

IF(?LA6ZT    -GE.    0.     !     GC    TO    81 

YVT  ,\u;     =-Pl   A?YT*P'   A  H7  T -a  R  /  X  P7T 


IFfYMINR    .GT.     YV.ALL)     GO    TC    81 
R^IMPT    =     Pi  XT7*R?/XP7T 


GC    TC    82 
8  1     RHTN?r    =     (  ( R*  PJL& BYT    +    PLA8 Z T_» Y  W A L.L  )»»2    +    PLXT2*YWALL«»2  )  /PLA3  VT-^-2 
82     IF1PLABZP    .GE.    C.)     GO    TC    83 

V  y  :  iv  R     -    -  P  ;  A  ^  Y  P  -■■  P '  A  .G  7  P  aR  /  X  ?  7  P 

IFlYflNR    -GT.     YWALL)     GC    TO    83 
3Jy_i-N?P    =    Pi  X£2«&2.ZXP.ZP. . 

GC     TC    84 
r^    Ry-M7P    =    (  (R*P;  ABYP    +  JPLABZPJtYWALLJ  *.*2._+_PLX? 2^YJ1ALJ.*_*_2.LZ.PJLA£.YPj?jl2_ 
84    RMN2    =    ANINKRf IN2T.RMIN2P) 

--(tat?   -   r^i\7    .;t.    r:-":v)    r.g   ~n    13C 


TIPT 


TLENG'd.     -    SETA) 


130- 


WTIPP  =  PLENG«(1.  -  YWALL/YPMHP) 
-14...  =..■■  I-MT  {  7  X-P-I/S4£P  I  KAS£*2  l-i 

ETIPT    =    EEL( IL.KASE+2)    +    (TIPT    -    ELLt  IL,KASE+2) )» (EEL ( IL+1     CASE+2) 
. =-£-£  .L4-IW{CAS£±24-U£  LSPCKASE+23 


=  INT(TIPP/ELSP{KASE*1) ) 

-)  D      —       C  CJ     '    T 


EluTPT    =    EELCIL.KASE+2)    +     (WTIPT    -ELL  I  IL  ,  :<ASE  +  2  )  )°  (  EEL  ( IL+  1  ,XASE+2  ) 
, =- E £ Ll  ZL,i<AS 2 ■:-2J-j  J SLSPt  KASS<-2) 


IL    =     INT(\rfTIPP/ELSP(KASE+l)) 

PJJ16J3      -      ~  Fl    <  T  I       •'  ic-:.;    '      ...       f  1-iT  T  OP     —  ' 


•;:,S3v 


f  Ll  ■--,  pX^SF^-1 


,KASE+1) )/ELSP(KASE+l) 

FVJTPT     —     '""T^     —     F  —  O- 


V '.'I A L '.   '     '■'-' 


"3    ft 


(PLA8ZP    -GE-    0-)    GO    TO    91 


■  o  -  o  •    ,-  r  7  o  .-<  r>  /  •  ■  o  -  p 


IF(YMINR    -G7.     YWALL)    GO    TO    91 

_c> ; ■■  1  y  7?    ~     °  L  ■'  P  2  ::-  ft  2  /  "  ^  Z  P 


92 


GC    TG    92 

_R* '  I ':  •'.  2  P    ~     t  '?.-"•?  L  A  BY  °    *N  ?  L  '  "  ?P -•"'•' '' A  r.  '•">"?     ■■•    P*.''P?-~V  '  A '  i  ^  ~-  7  ^  /  P '  A. ".  VP  •-  ■->  ? 

R2WALT    ={{R»PLABYT    -:■  PLA3ZT«YWALL )»»2    +    PLa72*YWALL**2)/PLABYT»«2 

-  c  '  p  •  ,•, 2  ~  ~     .  ^  j  3     Q  -  )    G  C  ~  0    c  A 


YMINR    =    -PLABYT*PLABZToR/XPZ 

TP./.VMTMfi         I   ^  .      " "'  MHO  \__TJ 


LiWJ 


L0_Z3- 


IF(YMINR    .GE.     YWALL)    GO   TO    94 


'JMZI- 


"2*32, 


GO    TO    95 

GC    TO    95 
a  a    P"TM2T    =    °?'-' ^  T 

95    R^IN2    =    AMIN1(RNIN2T,RMIN2P) 

IF£CAT2  -  Rr-i.\2  .lt.   difmzn") 

r-n 

m 

:•  ~n 

pt  =  tle\g-:l. 

OD      —       Dl    tllfi.'iM 


-   YWALL/YTHHP) 


tipp  =  ple;-;g-  ::.  -  alpha) 

-,     _    iMTH-iTTPT/F'  SP(kasp+7>  1 

EWTPT    =    EELJ  IL,XASE->2/    ♦     (WTIPT 

■2)  )  0  i  EEL  t  IL+  It  i<ASE->2  ) 

il  =   i\t:;:ti??/e^s?j:<;^se-:-:;  ) 

euTOD    -     -r  -  t   1  t  •    „  "  :■  c  p.:  •  •      .-.     ■•'.'-- -OD 

.,,!„,  PP.    m+l-KASF*.!! 

-    EEL{IL,KASE+1) )/ELSP(KASE+l) 


'IP?   =    EEL(IL,KASE*1)    *    (TIPP   - 


ELL; £L«KASc+1) 

•c-->"  • 


)*{EEL(IL+ltKASE< 


IDEPQS    =    EWTPT    •>    EWTPP   -    ETIPP 


100 


'2    oLE. 

•  •    ■_.  o  •■-  R  ' 


CAT 2)     GC    TO    110 

A    .G"n    v  •■■',!  '   )    r-Q_ 


F(PLA3ZT    -GE.    0.)    GO    TO    101 

'•jtm'i    -    _p-  .'-"/t.it  A"7~  ■•■■"' /'  1" LLT 


FtYMINR  0GT0  YWALL)  GO  TO  101 


GO  TO  102 


■  c? 

RfIN2T    =     ((RoPLABYT    *    PLA8ZT*YWALL )«»2    ->    PLXT2*YWALL*»2) /PLABYT»»2 
R7WAI  P    =     MRftP?    iiRYP     •■-     PI    iR7P*YW     II    l«?     -■■     PLXP?*YWALL*»2 ) /PLABYP*»2 

IFtPLABZP    .GEo     0.)     GO    7C     104 
VVTNR     =     -PI  ARYP-aP!  AH7P«R/XP7P 

IF(YMINR    .LE.    YPMHP)     GQ    TO    103 
TFfYMTNR     .r,F.     VWAI  I   )     CO.    7  ft     t  04 

RKIN2P    =    ?LX?2«R2/a?Z? 

r;  r    T  n    7  n  ^ 

103 

RPIN2P    =    R?2 

cr   to   :ns 

104 

RNIN2P    =    R2WALP 

Rf>T,M?     =     ACTi\WRMIiM7T.RNTNPP1 

IF(CAT2    -    RMI~N2    .LT.    DIFMZN)    GO    TO    ISO 
I..ITTPP    =     PI  FNft«  M  -     —    vi.  ••,■•    f y p m H °  ) 

WTIPT    =    TLENG«(1.    -    YWALL/YTMHP) 

-rpT     =     Tl    ->JG-:>  '  •  .,     -     RFTA1 

IL    =    INT(WTIPP/ELSP(KASE+1)  ) 

FNTPP    =     FFI  f  T!  ,t<ASF-;-T  5     -:-     '  i  -  7  P?    -  -\. "    I  U  ,  !</  SE  +  1  5  )  «•  (  -?.■    1  1 '  +  1  .  i<  A  S  F+I  / 

1                         -    EELUL,KASE+1)  )/ELSP(KASE  +  l) 

7'       =      T  VT  (WTT  PT/Fi    S  7>  I  :^  &S  --:-?*   *> 

EUTPT    =    EEL;  :L5;<ASE-:-2)    -:-    (WTIPT    -ELL ( IL.KASE+2) )»( EEL ( IL+1 tKASE+2) 
-                             -    FFI   <T!  .K/JSF  +  PM/Fi.SPlKASF  +  Pl 

IL    =    INT(TIPT/ELSP(KASE+2J) 

ct-pt    =    FF!  m..BKASF+71    +     (7IPT    -    ELL  {  IL  .XASE+2  )  ) «  ( EEL  ( IL+ 1 .  KASE+2  ) 

-    EELUL.KASE+2)  )/ELSP(KASE+2) 

FRFPn^    =    r''"?'3    +    F'-'T0"    —    FT  T  PT 

GC    70    120 


o  a  3  •    \  r  v  p 


PLA"?.?---V'--:^'-'-  V 


?:.  X??qYV.'A,.A<h>2)/PLA3Y?*:>2 


R2WALT    =     UR*PLABYT    +    PLABZT*YWALL )»*2    ■»    PLXT2*YWALL»»2) /PLABYT»»2 

"  E  ■  ?1  AR7P    .r^.    fln!     GO    70     ;  '? 

YNINR   =   -PLA3Y?*PLABZP«R/XPZP 

tp   '  y  ;  •  -  \]  Q      _  ■   p_      VPMH  9  ^      f7  0     "0      "71 


if: 


Y  M I N I 


YWALL3     GO 
>?--R?/Y>-7? 


c-r 


X2? 

to   : 


=    AMAXKRP2.R2WALP) 


'JLA 


N2? 
y-  p 


=    R?2 


1  1  ? 


70   : 

M2? 


=     R?' 


»*\r  rt 


X2? 

PI  AB 


=    R?2 

77   .(;-, 


0-)     Gf 


NR    = 


115 

lYT-oPLABZT*R/XPZT 

.    YT"H°)     GO    70    11 


YMINR    -GE.    YWALL)    GO    70 
dill 


PIXT?«R7/Xi 


X2T 

~n    i 


=    A,-;AXi(RT2,R2VJ. 

-6 


75 


114 


N2T 


=    RT2 

=    R?" 


_LL5_ 


bo 


TO    1 
'■'RT 


16 

=    R?'. 


=     K12 

LK2Pa.ftfIN?T) 

A!*AX1  CRHAX2P.RMAX2T] 

7  -  r,:r>i7  „  i  T„  Q  "  ^  7  -  ?M  ^Q 


1  ?  =  PI 


-  YWALL/YPMHP) 


-132- 


ewtpp  =  eel; z.oXase-: •:. )   *    :■::::         (  zl,xase+U)<»(Eelul+i,kase+i] 


■  S?CXASc  +  H- 


:nt<vjtip7/elsp(kase->2) ) 

- f-ei    f_Tj    .  (/  A  <:  g  J-  7  >     •:- 


-    EELUL»KASE+2  3  )/ELSP(KASE-5-2 ) 

ECEnCS  ~  "-'.'~ ° °  •'•  ^''"p" 


120 


INT CEDEPOSeRCHNWD) 

j.f_r_i  =  r.v.  _ .v  f." '  -:-  "~ 


-CHA- 


WTSM(NSUEN)  =  WTS.v{.\EUEN)  *  WT 

-2-Q.  0)  GO  T  0  14 Q 

GC  TO  549 

?  0  LT  S  T '.  .VE'JE.V  1  ~  '..""ST  ' v  -''"V '  •••  * 


WTLST(NEUEN)  =  W7LST<NEU£N)  +  W7 

c  ?  ^  c  ~  c  cr  .  ?  n  .  HI  h  °  7  H  "AH 

-X-!  £B — = — ',W-Z~£-0 

YY8  =  Yi-'S-O 


-Z-2-3-^- 
PLAPP 


■=0 


>£SV; '"  —  \-'CSYS"~C 


VCSYYY  =  VCSYSYeO 

V  T  '* v/  ?  I?    ~    '■ .' r  ^  v  c  7  •"•  ° 


1 1 R I T  E  { 3  »  5  4 

_£^-3— F-C^hat  'A2TLH 


}    XXB,YYB9ZZB»VCSYXX, VCSYYY, VCSYZZ,PLAPP,PLAT7,Y0UT 

LCS.7 i^^L^llZ^7'^^    =    (?; '7„7;;^);r;.?--.Hl    VCSY 


5481S  =  (2{F7-2,!H,) ,F7.2.9K)  ?^A3?  =,F7.2»8H  PLABT  =,F7.2,7H  YCU1  =  • 

-&4&2 S-Z^Ai 1 


549 0  =  1.E15 

oodv  —  oorw«! 


XPXM    = 

XPCM«0 

Z  P  Z  i-'    = 

ZPCM«0 

PLAYYP 

rs  i    ',  7  ■»  n 

=    PLABYPoO 

-      0!    AR70-in 

^Lma/>  1 

-             -     *  ,   V  f  — 

=    P  L  A  3  X  T  ■"  0 

PL-ZZT 
WRTTF ' i 

=    PLA3ZT--0 

^n)    -p.^„!<^<;c.rn<;pH-,c-^:pu-.rT!--rT:-p-^H"P-p??r-'.x?x:'?v?Y:-u 

1ZPZM,PL 

7  Y  P  r  V  P  j  7 

AXX?,PLAYYP,PLAZZP,PLAXXT,PLAYYT„PLAZZ7,EP,ET,PLENG,7LENG, 

P»  ITT  ..  VT  .  7T„  PP?  .  RT?„  A'   PH  A  ..  RFT.»i„  P  ?*■' A!  P  „R  ?' '  ■'  :  T»P,MIN?P-RMAX?P« 

3RMIN27, 

=;  =:  r>     p  r  P  '■'  i  T  ( 

RHAX27  jYM.INR,!  IPP»7I  PT9V17I?P, W7IPT  ,  Z3EPCS 

n?n     •-:?       r   kasf   r.sp   snp   r.TH   qtp   ?:-"?   ppt.mfi5  x?c;-'       ypcm 

550:       ZPCM 

PLA3X?    PLABYP    ?i.A3ZP    PLABX7    PLA3YT    PLAB2T    EPCE65     cT(Ec)     PL 
PN«/PA.?.P«i-?.  :?-•••  iX.F?.?t.F«5.5.           mF7-^?fiPPF7^.0P?F7.A/ 

55 


C3: 


i2H 


YP 

ii",ipo 


Z?     XT 

';  po  c  •'  MP" 


":'"?T  '/:■!-■ 


^P2 
•  pp 


RT2 

T  T  P" 


ALPHA 

1.JTT  DP 


i  C  .  «  r\.£ 


5505E 
"  ah  : 


'CS/6F 

:  :r:-!;T 


6.2?F7..2?r6<.2?2r6»3;>6F6.2jb>.~6o;>j 

p<\'  -■  y.  mit.-vs;  '  r,->:   ---n  -? 


LET  NUTRNS  3E  THE  NUMBER  OF  NEUTRONS 

n;-  nv-,f.ci'».i>^  r^\T'!-^v pp.  ^nw  amy  n" 
FRONT  CR  BACK)  IN  ONE  SECOND.   THEN 


ICSS  A 


RFACc 
THFR 


;E  NUMBER 


I  TCP  =  NCHANS  v  I 


-133- 


'  -'•  - 

ITOP    =    ITOP  -    : 

-  -  ,  r ,_,....  -  —  _. .     _  f;-f  _    n  ,  i    r;  -    — i    •  •  • 

IF { I  TOP    .EQ.     1!     GO    TO    142 

r.n     TO     1  A 1 

142 

"  A  ". 

WRITE (3. 143)     NEUEN 

prRMAT/77H     Nn     r, c     ppR     N'FIIFN    =»T->,! 

CHAN (I)     =     1. 

WT^M  t  MPI IPM  1      -      1  - 

1  •';  4 

r 

J    =    MOD  (NEUEN-I  »33     ■:-    1 

TTHP     7  "~     '":-'""     TI"l->     f.HANNF!      THAT     ,_.  •  c     ;...,,     [-("tllNTS. 

J    RA 

]ap<.    pRnM    . 

"0  3.   IV  IS  ONLY  TO  ALLOW  THE  RES-JLTS  FROM  3  CYCLES  ON  .NEUEN 

•n  -.:-  st-r-r  .-.;-:.-;  t:--\'  r:'.'--.'  .•.••_  at  rare  tr  cni-:s-RV'E  ?A?E^a 


:j(J)  =  NEUEN 


n  "■     U"      I     =     1  .  N r*rl  /!  M  s 

CCU.\TS:  Iff  J  5    =    CHAN  (I  J 

14-5 

CNORM(I«J)    =    CHAN(  I  )-"-1000./WTSM(NEUEi\l ) 
C-fiNTTNIJF 

r 

CPIT?     IS    PROPORTIONAL    TO    CHAN     AND    FIXE"     SO    THAT    THE    TOP 

CHANNEL 

c 

IS    1000. 
WRITF^-JAd)     NFUFN.  TNTFRPTFN?PAR7,  SIP"5,"  "*D^  IG~C~  -  CHC  I  C  -  -1 

ITFf.TR. 

1  PNEUTL,VCMSYS,PNUTCM»ENCH»PHE3CM, EHCMsETCM, I  TOP, CHAN (  I  TO?) , J 

sap  g,-^;/AV  •//:??:-:  w:=;;  =  -;  -NT"-*   -n  -> -.rt  sjr-?  rigo  s-r-qt  choice hie. 

560ICTR  PNEUTL  VCMSYS  PNUTCM  ENCM  PHE3CM  EHCM 

-~r,->  --'v/      ?'<=■  ---'.  ^  .  "-■•  .-=;.?-  ""p.,-,.  "-•  r  „?/-■•-:    — -~>   =  . -a  . -.-.;-?   r;-av  ' 'ATP ' 


5601 


E12.4,5H       J    =  ,135 

.■    ,.:■-„     v    „.■•■■-;„    r'EiiEN    „M".    v.- MS?    G" 


WRITE(3,146){{KJ(K)«L    -    1  s  3  ) ,  K    —    _  ~  J 5 

.aa    •r";."".l,-y       :      ~  .'  hf*.'  -, .  ?  ■'  -  -  h      r':;::i"(; 


JLZ^ 


1461  11H)       CNGRM(I,,I2, 1H> ) 3 

Tn-o-    =     --.no    ■••    -; 


IFIIT0PP5    oGTo    NCHANS)     IT0PP5   =    NCHANS 

nn    -  4fi    t    =    ■  ,'-■'-->-)'; 

ENERGY    =    ENLIM«FLOAT{  I)/{  1.E6'->cLV0LT*FL0A7  { NCHANS)  ) 
URTTFf*.1471     v  .cN?RGY.  C COUNTS (  I . K ) .CP ITP ( I ,K> .CNORM ( I .Kl . K 

. 

\ 

147 

1  A? 

FORMAT ( 14,F8.2»3(E15.5,2F11.3) J 

r.n\iTT  niif 

:<  J :   =   K  J  ( I ) 

K.I.I    =     KM. 1) 

J  AO 

WRITE(3»149)     (K»WTSM(K),K»LOST{K),K    =    :<J-;:<JJ3 
-r--'f-'//-f^^    WT<;M{..T?.^H)    =„F14.,fi,7H      LGSTC  9  1 2.  3H1     =,173) 

'/.'RITE  .'3,570  )     ('.<,l,'TLST(K),K    =    KJI,KJJ) 
!-iRTTF(^.S7--  1      (  K  .TIYTA!    ;-  •  ,«     =    K.r.K.J.J! 

570 
571 

F0RMAT(3(  19X  1 6HWTLST  (  ■>  12 1  3H3    =»E14.6)  ) 
FORMAT  (3<  !9X?AHT0TAL(  ,  T^.-JH)     =.F14-6n 

150 

IFCZFPLOT    .EQ.    0)     GO    TO    160 
Rr    *55     V    =     I, NCHANS 

IF{CHAN(Ii     .GT.    0.)    GO    TO    152 
:  r*Hi\j  i  -  •■     -    _occc 

]  <"? 

GO    TO    155 

r.wN  =   a:  rr-o  'Ghan;  t  ) y 

LCHN(I)     =     INT  (SIGN  :  10C0o--A3S(CnN  ;     ♦    . 5,CHN)) 

[■tir.HNin    .lf.    t--oooon   lchnii)   =  -<;<;<;<) 

155 

r 

CONTINUE 

ithm   -<-    innn   y    n  orartthmc  ra<;f    ?.n)   of   chank 

WRITE ( 14,1)    NCHANS,NTA8S»NENS» IFTES7» IFPLOT»NUTRNS» VOLTS»TI 

ME  i 

134- 


ELOW,EHIGH,£N 


,  NO  ANSI 


157  rCRyAT(I6I5) 
4-60--C-CNT4NUE 


CALL  CLOCXI (TIM3) 

:;:XIXZ  C3.» 5 75 J— 7- 1  :-U U-M2 «£IK3 —  -  

575  FCRMAK7H  TIM1  =,F13.2,7H  7IM2  =,F13.2,7H  71X3  =,.-13-2) 

^ZQ2 


-135- 


SIBFTC    LENGTH 

^  ftO."iiT1 


"APS,,  FVI    TM. 


AJjKj 


:i  ,ya: 


cp 


) 
D I  M£J\LSLCAL-£K£S.LAOjg)  ,  ER.GL60.01  ».RPA1FU  fcjQI 


.EL16HC  tAJ_. 


ENSIGN  ELNAXK) , ELSP( 4 ) , ELL ( 600, 4 } , EEL ( 600,4 ) 

. F  N  7 . .  aS_=_i=  LCA.IXM  A  7;  SJ 

ELVCLT  =  I. 602076-12 

R'-TIH^  = .  U-V-.29 


RHUKR  =  3-455 

?<:;  =  10. 

J^.CX _=  .  PH£*_RHQhJE  _*.  PKR*RHCKR_. 
WHE    =    PHE*RHOHE/DENCM 

'.-  K  O      -      I   -      -     '•  '  — 


D( 


31±U 


1 , NT  AGS 

JCQe_CJli'J>L,TES_ 
THE    LENGTH     IN    THE    A 

JZ5— I-UC:.L- 10  7  7^}  :     W.I 


PRCTCN,TRITCN,DEUTI 

:  p  m  -  f--  -  - ;  r.  a  t  f  t  w  f  v 


TRACK  LENGTH  VS  .. .  PAR.TLC.LE  .ENERGY.  ANC_ST0RES_ 
RSAY  EL.   THE  FORMULAS  ARE  FRCM  WANG.  P.73- 
Ih_.SLl£HT.  CC:^EC.7IC.\Sa.._J._=_I,Z,3,4  - .IMPLIES. 
RON, DEUTERCN,  RESPECTIVELY. 

-n*p;  vn:  -<-->:7A7,s  :■ , 


E  N  E  R  (  I  )  =  E 
.£&£-(  I  )  =  ENLUE1 
RPKR  =  RPAIR [ : 
RPHF  =  RPATRfT 
J3  =  INT(FLCA7 
.■?     =     -\!T(  -^  TAT 


aFI_DAILL).Z£NT  A  3S 

)*(1.89    -.25*ALCG10( E) )    -    . 
L5l«,82_+...04.3.*AL.0GL0.tE-L) 

ID/3.    •:•    .001) 

LI  )/?^L^.ClCL-J 


36 


PAR3    =    FLOAT (M 
par?    =    Fl  PAT  . ; ■• 

rpair3  =  rpair 
-R£AJ_R2_=_RPA lr 
rtkr  =   3.mrpa 

PTKP       -       ?,  .    ,;.  S  D   ■    ' 


UC(I,3)  )/3. 

QC.i.1,.2)  )/.2._    ___ 

(J3)  •:•  PAR3*{RPAIR(  J3  +  1)  -  R?AIR{J3)) 
LJ21_+_.P.AEL2*  LRPAIRU2*L)_Ti_RPAIiUJ2)  I. 
IR3M1.89  -  .25«ALOG10(E/3. ))  -  .36) 

P. '-I- :■!„;-•;  7  -.-  .0^7--- A'  P"'  .?■'  -/?  „  )  ? 


2.-MRPA 


RCKR 

JIC.RE. 

R7P  =  1./ (WHE/ 

ELLI  =  i  -/  ''/:-•:-/ 


1 .  /  (  W  H  E  / 

=  v-  7  ?  / 


IR2*< i.8< 

^2*_L.C2 

RPHE  +  Wl 

KJ.H£_v_..h. 

RCHE  -:•  kKR/RCKR) 

n  -  N  H  '•' 


>9  -.25*AI 
iKR/RPKR) 


nGi0!E/2. 5 )  -  .36) 
.aGLOCE/2^1) 


l:  :  ,2)  =  rtt/ 

:M  '  -3)  =  RTf!/ 
iLH  ,4)  =  EL(I 

■  CNTINLLE 


DENOM 

» 3  ) 


DC  7 


J  = 


,4 


?MTAP.S. 


!  L  S  P (J )  =  ELM 


iX{ J)/FNTABS 


EL(I  ,J)  =  0. 
F!  i  n_^J_L_.=..  Q^ 
EEL ( I  ,  J  )  =  0 . 

rr    fi  :  =  ■  .  g 


i  L  ( I  ,  J )  =  E  L  ( 1 

-1  ;t  -->)    FOR 


FORMULA  GIV 

A  rrMT'MlF 


at; i)/ io. 

_AJ?E_.1N7.ERP.0LATEC  FROM  EL  !  iC  ,  J  )__  BECAUSE  .THE 
ES  NEGATIVE  VALUES  IN  THIS  VERY  LGW  ENERGY  REGION. 


0, J) »FLC 
_L_=..L-_9 


:ntinue 


136- 


S  FORMAT (49H1 


E(MEV)  E*E6(ERGS)  L-P(CM]   L-T(CM)   L-0  (CM  )  ,10X  , 

„'^i:s:  l-7 :r;;:  e-;_-t-,  :- e.s:e?.gs:  l - c r c ; : : .  f  ;:.-'} 


S2£6(ESGS: ) 

E-LU-UJl-IS  TRACK  LENGTH  OF  gA.R_TICL-c.->L»£LLXNTAB.S,J  !  =£L  :  .MAES .,  JJ 

THE  DIFFERNECE  BETWEEN  EL  AND  Ei_L  IS  THAT  £LL  VALUES  VARY  . 
L-I-A£.AJU.-Y— £P~0;V— i_T-fl-M.A8S.»--UK.I4_£.  E_  i  I,JJ.__IS— THE  7RACX_L-EN£-Tl±_aE- 


"HS  J-TH  PARTICLE  CE  ENERGY  ERG(I).  WHERE  ERG  VARIES  LINEARLY 

■^i__  ■■■    •    "2  ■■ --  o  <:    "  :'_  '  '.*  '"-  ~  •:''-'^"  ""'-.',""  rr-^pr.vn1;  "^  c'  ■  - 

"HE    ELL-EEL    RELATICN    IS    THE     INVERSE    CE    THE    ERG-EL    RELATICN.     THE 
■X3CP    TC    STATSi-tEflT    IS.  -CCKP-UISS    AM£— R3INTS    THESE    ARRAXS. 


DC     15 


7\TABS 


ELLC  I  ,J)  =  ELSP<  J;  »FLCAT<  I  ) 


F(EL(ITRYjJ)  oGT.  E L L ( I » J ) )  GO  T C  i G 
■£4-X-T-R-Y— ^GE^-N-TABS)  GC  TG-JL2 


IF(EL(ITRY+1,J)  .GE, 

GC  TO  0 


ell{  :,o  ) )  gg  to  :: 


(ELLCtJ)  -  ELUTRY,J)  )/ < EL { ITRY+1, J  )  -  EL(I7RY,J)) 


-  ■z  v I  I  ••  .-.  i 


-2Aa.T-L/-FN7ABS 


G  C  T  C  13 

DAP  ~   -   I  : 


-     "=■    '  -TRY-1  ,A)) 


13 

EEL(ItJ)  =  ENLIM*{  FLOATt  ITRY-1  )  +  PARTJ/FNTABS 

CONTINUE 

E  E  L  ( I  ,  4  )  =  E  E  L  (  1 . 3  ) 

WRITE(3»14)  I tENERCI) »6RG(I)» tELt 1? J)»J  =  lf3)f(ELL(I»J)»EEL(ItJ)i 

14 

FCRMATC  I4,F9.2,6PF9.3,0P3F9.4,8Xf3t0PF12.4f 6PF13.4) ) 

RETURN 

This  report  was  prepared  as  an  account,  of  Government 
sponsored  work.   Neither  the  United  States,  nor  the  Com- 
mission, nor  any  person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation,  expressed  or 
implied,  with  respect  to  the  accuracy,  completeness, 
or  usefulness  of  the  information  contained  in  this 
report,  or  that  the  use  of  any  information,  appa- 
ratus, method,  or  process  disclosed  in  this  report 
may  not  infringe  privately  owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to  the  use  of, 
or  for  damages  resulting  from  the  use  of  any  infor- 
mation, apparatus,  method,  or  process  disclosed  in 
this  report . 

As  used  in  the  above,  "person  acting  on  behalf  of  the 
Commission"  includes  any  employee  or  contractor  of  the  Com- 
mission, or  employee  of  such  contractor,  to  the  extent  that 
such  employee  or  contractor  of  the  Commission,  or  employee 
of  such  contractor  prepares,  disseminates,  or  provides  access 
to,  any  information  pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with  such  contractor. 


